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Preliminary Remarks 


The content of this compilation comprises our early work on water diffusion and swell- 
ing in silica between 2011 and 2018. Newer results will be contributed in a second book- 


let. 


Water diffusing into silica surfaces gives rise for several effects on diffusion behaviour 
and mechanical properties. The most interesting observations from our point of view may 


be listed here. 


(D In absence of externally applied stresses: 


a) Density reduction with increasing water concentration [1]. 
b) Increasing concentration of molecular water at silica surfaces [2]. 
c) Decrease of measured water diffusivity with time [3]. 


d) Change of the shape of the water concentration with time [3]. 


These points are covered in the present booklet (Part I). 


(ID In presence of externally applied loadings: 

a) Increasing diffusivity under externally applied tensile stresses and a decrease in com- 
pression [4]. 

b) Decrease of water solubility in tension and an increase in compression [4]. 

c) Increased water concentration in tension at high temperatures [5]. 


d) Strongly increasing strength at high tensile stresses even under very low water va- 


pour pressure [6]. 


These issues are further elaborated in Part II. 


For an explanation of all the findings listed under (I) there is only volume swelling by 


hydrogen generation an appropriate tool. 


IX 


Preliminary Remarks 


In order to give a transparent description, two main assumptions will be made: 


(1) For an interpretation of the increase of diffusivity and water concentration under 
tensile loading, a further effect may play a role. Under tensile stresses, the nano-pores of 
the ring structure are possibly widened and diffusion and solubility may be enhanced. 
Such behaviour was observed in different materials and is interpreted in terms of the so- 
called free volume theory. This may also affect water diffusivity, generally described by 
the Doolittle equation [7], and solubility (see e.g. [8]). Due to the lack of results for sili- 


ca, free volume effects must be ignored in the following considerations. 


(2) Water molecules are believed to diffuse through silica glass in much the same way as 
the noble gases. Hence, the diffusion of water molecules is also expected to show little 
dependence on applied stress [5]. Following this reasoning, we believe that the molecular 
water in the glass is located in holes of molecular size in the silica glass network where 
the water does not contribute to the volume of the glass. Consequently, molecular water 
cannot contribute stress effects. Thus, the concentration of molecular water in the glass 
depends only on the vapour pressure in the external environment, but not on the stress in 


the glass. 


In the APPENDIX, analytical and Finite Element results (FE) on cracks and pores in 


glass are compiled. 
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1 Diffusion in the absence of external 
stresses 


1.1 Chemical reaction 
Water penetrated into silica reacts with the silica network according to 


=Si-O-Si= +H20 <+> =SiOH+HOSi= (1.1.1) 


with the concentration of the hydroxyl = [=SiOH] and that of the molecular water C= 


[H20]. The equilibrium constant of this reaction is at temperatures of 08<450°C 


ml (1.1.2) 
C 
In the range of 90°C<0<300°C the equilibrium constant can be expressed by [1.1] 
k, = 323eup{ 702i mel) (1.1.3) 
RT 


where 7=0+273K is the absolute temperature, and R the gas constant. The equilibrium 


constant is at temperatures 0>500°C 


2 
kas (1.1.4) 


1.2 Water diffusion 


Due to diffusion, concentrations of both molecular water, C, and hydroxyl water, S, de- 
crease with increasing distance, z, from the glass surface. If the reaction given by 


eq.(1.1.1) is in equilibrium, the diffusion process is governed by the diffusion differential 


1 Diffusion in the absence of external stresses 


equation: 


2 
OC _ pee (1.2.1) 
ot Oz? 
where Dis the effective diffusivity that takes into account the chemical reaction given by 
eq.(1.1.1) [1.2]. The diffusivity shows an Arrhenius type of temperature dependence. 


Solution of the diffusion equation requires an appropriate boundary condition, very often 


chosen as constant surface concentration of molecular water: 


C(z = 0,t) = C(0) = C, = constant (1.2.2) 


For a semi-infinite body, this assumption results in the following equation as a solution 


of the diffusion equation, że., eq.(1.2.1): 


c= cert u (1.2.3) 


2./Dt 


which is well established for liquid water as the environment [1.3]. 


In the case of water vapour as the environment, the behavior of water at the silica inter- 
face does not follow eq.(1.2.2). Experimental results show that the concentration of 
water is not a constant at the surface, but increases continuously with time, until it reach- 


es the solubility limit of water in the silica glass [1.4, 1.5]. 


Figure 1.la shows both the molecular water and hydroxy] water concentrations of silica 
glass at the surface (z= 0) for a temperature of 250°C and a vapour pressure of 
39 atm [1.6]. The fact that the equilibrium ratio between =SiOH and H20 is achieved in 
these experiments is supported by the fact that the ratio, S/C, slightly below 3 is reached 
after a relatively short exposure time, ~200 h, Fig. 1.1b. However, the fact that C(O) has 
not achieved a constant value in Fig. 1.1a suggests that the solubility limit of water in the 


glass has not been realized. 


The experimental results in Fig. 1.la [1.6] for ((z=0, d are a clear contradiction of the 
assumption of a constant surface value for the water concentration. One cannot assume 
that the surface concentration of water, S(z=0) or C{(z=0), stays constant under the con- 
ditions used for the experiment. There is, however, an incredible agreement with the 


thermal analogue of thermal shock behaviour under heat transfer boundary conditions, 


1.2 Water diffusion 


governed by a heat transfer coefficient. This calls for a surface condition 


oe LG 26) E (1.2.4) 
dz D 
where again G is the concentration of molecular water reached at z=0 for o. Follow- 
ing the suggestion by Doremus (Section 4.7 in [1.7]), the phenomenological parameter A 
in (1.2.4) may be interpreted as a reaction parameter for a slow surface reaction that 


limits the entrance of molecular water species. 


On the other hand, a simpler phenomenological description is possible by assuming that a 
barrier exists to the transport of water across the surface of the glass (for a physical de- 
scription of such a barrier see Section 1.3). The barrier gives rise to a mass transfer coef- 
ficient for diffusion, which slows the passage of water into the glass. Each of the assump- 


tions yields the same set of mathematical equations. 


C0 hydroxyl water S/2 _-- ® 
®t a) oo a| b 
S a 
(0) .. SIC 
Pa e 
1 
0.57 a ail 
A o | nE _ u --1----# 
a 8, 0 tga r 
O 0.47 4 
a / f E 
5 molecular water C 
= 0.3} ay 
= [ 
2 
0.2 
it 
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1 1 1 1 1 1) 1 


0 f f 1 N 
100 200 300 400 500 600 0 100 200 300 400 500 600 
t h t h) 
Fig. 1.1 a) Molecular water (C) and hydroxyl water content (5/2) at the surface of silica for 250°C at an 
equilibrium vapour pressure of 39atm, b) ratio S/C from [1.6]. 


As shown by Carslaw and Jaeger (Section 2.7 in [1.8]), the concentration profile, C(z) 


resulting from the boundary condition for a semi-infinite body is given by 


1 Diffusion in the absence of external stresses 


Zz h h? Zz t 
C(z, t)/C, =erfe ex z+—t lerfe +h (1.2.5) 
an is E D l Br 5 


At the surface, z=0: 


h? le 
= u Zell, (1.2.6) 
C(O, t)/C, =1 es} ey | 


For reasons of simplicity, we introduce a normalized dimensionless time t and normal- 


ized depth coordinate 6, defined by 


h’ ` Z (*) 


en ar: 


Equations (1.2.5) and (1.2.6) then read 


Sen = ete $| ~exp[¢vr + ein] $ + a (1.2.5a) 


0 


and C(0, r)/C, =1-exp[c]erfe[ vr]. (1.2.6a) 


The C-measurements of Fig. 1.la, carried out on specimens of the order of 1mm thick- 


ness, were fitted for the molecular water according to eq.(1.2.6). The result is 
Cy= 0.62 [0.523, 0.72] (mol H2O/litre SiO2), AN D=0.075/vh [0.046, 0.104] (D 


with the 90% Confidence Intervals given in brackets. This result is introduced into Fig. 


1.1a as the solid curve. 


The measurements of Fig. 1.2 were carried out by Davis and Tomozawa [1.4] and Waka- 
bayashi and Tomozawa [1.9] at a temperature of 350°C under 355 Torr water vapour 
pressure. The specimen widths, 2 W, ranged from 200 um to 1000 um. Since the maxi- 
mum penetration of water was less than 70 um (Fig. 4a in [1.4]), the solutions by 
eqs.(1.2.5) and (1.2.6) are applicable for these diffusion results. 


At 350°C, the equilibrium of the S/Cratio should be reached earlier than in the 250°C 
tests, because of the higher temperature. Consequently, S/C should be constant as sug- 


gested by the experiment carried out at 250°C and indicated in Fig. 1.2. The data points 


1.2 Water diffusion 


were obtained from IR-absorbance measurements and represent the IR-absorbance of the 
fundamental OH stretching band at wave number 3673 cm" per light path length through 
silica. In this case, the hydroxy] distribution and the molecular water distribution can be 


described by the same equation: 


S(z,t) _ C(z,0) | (1.2.7) 
S, C, 


where So is the saturation value of the hydroxyl concentration for >». 
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S(0) 
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Fig. 1.2 Hydroxyl surface concentration after hydration at 350°C in water vapour at 355mm Hg; measure- 
ments were made by Tomozawa et al.[1.4, 1.9], (ordinate scaled in terms of the IR absorption coef- 


ficient). 


By curve fitting eq. (1.2.6) to the data points in Fig. 1.2, the following parameters were 


determined: 
S= 18.9 cm” [8.3, 29.5], AIND=0.0315/Vh [0.0031, 0.060] (ID) 


The result is introduced in Fig. 1.2 as the continuous curve. 


1 Diffusion in the absence of external stresses 


For a temperature of 200°C, surface equilibrium between CO) and S(O) is not expected 
in the short exposure times of 5-20h reported by Helmich and Rauch [1.5]. The hydroxy] 
content at the surface, S(0), develops with time according to eq.(4.25) in reference [1.7]: 


dS(0 
un = 7C(0, t)- uS(0, t) (1.2.8) 


Re-arranging yields a linear differential equation of first order 


dS (0, t 

BED, aso, t) = yC(0, t) (1.2.8a) 
t 

The two coefficients u and y can be related to the equilibrium constant X for the SC 


ratio, the ratio of saturation concentrations 


So k (1.2.9) 


Ge ess i (1.2.10) 


It should be noted that the equilibrium constant X is stress-dependent. It is well known 
from the work of Le Chatelier [1.10] that chemical reactions that exhibit a change in 
volume will be sensitive to the ambient pressure. This also holds for mechanical stresses 
as will be outlined in Sections 5 and 6. When X is increased for instance by mechanical 
stress, the ratio A~y/u must increase, too. It has to be expected that the forward reaction 
rate x will be increased, the reverse reaction rate, u, will be reduced. Especially for ten- 


sile stresses in the GPa region, the stress effect on reaction rate can become strong. 


Since the concentration QO, A on the right-hand side of eq. (1.2.8a) is known from eq. 
(1.2.6) as the surface boundary condition, and the surface boundary condition S(O, Ò is 
known, too, namely S(O, Ò = & QO, Ò. The solution (1.2.8a) can be obtained by using the 
Mathematica Routine DSo/ve [1.12] resulting in 


S(0,t) _ 


> h + u(1 -expla t]erfce[ayt])- a exp[-ut](æ + Juerfil utp} (1.2.11) 
a +y 


0 


1.2 Water diffusion 


with the abbreviation a = A/VD and the imaginaty error function erfi. The total water 


content in molar units C, is then given by 


C 0, ) =C(0, t)+ŁS(0, t) > C,,=C,+S,/2 (1.2.12) 


In Fig. 1.3 the molecular, hydroxyl and total water concentrations at the surface are plot- 


ted versus normalized time. All concentrations are normalized to their saturation values. 
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Fig. 1.3 Surface values of the water species C, Sand the total water Cy as a function of time for an equilib- 
rium constant of k=2.1, valid for 200°C. 


Water concentration measurements were carried out on silica by Helmich and Rauch in 
[1.5]. The silica specimens were soaked in saturated water vapour at 200°C, 15.3 bar 
(1.53MPa). Fig. 1.4a shows the surface concentration of the water from H-measurements 
using the technique of Nuclear Reaction Analysis [1.5]. From data fitting and using 
So/Co=2.1 it results: 


Q=1.71x10” /cm’, a= 1.902/vh, u=0.130/h 000) 


From © we obtain: 5=2.1x@=3.42x10”"/cm’, Cwo=2.05 Q=3.51 x 10% /cm?. 


1 Diffusion in the absence of external stresses 
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Fig. 1.4 a) Total water concentration Cy at the surface at a temperature of 200°C in saturated water vapour 
by Helmich and Rauch [1.5] (data points), Curve: eq.(1.2.12) fitted to the data, hatched area: Re- 
gion of measurements by Zouine et al. [1.3] carried out in liquid water, b) hydroxyl and molecular 


water concentrations, c) ‚S/C-ratio vs. time. 


Since only six data points were available for the fitting analysis, and four parameters had 
to be determined, the results do not have high precision. Hence, confidence intervals 


were not determined. The two water species, C and $/2 are shown in Fig. 1.4b. Finally, 


1.2 Water diffusion 


the time-dependent S/Cratio is represented in Fig. 1.4c. From this diagram, we can ex- 


pect that S/Cequilibrium is reached after about 30 h. 


Surface concentration for tests in liquid water 

Water concentrations in silica glass exposed to liquid water are reported by Zouine et al. 
[1.3] for temperatures ranging from 25°C to 200°C. For several of the temperatures, the 
authors measured water content after different times. The ratio of longer to shorter 


times ranged from about 2 to 4. The related concentrations ratios, Chon/ Gong are plotted 


in Fig. 1.5 as a function of temperature. 
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Fig. 1.5 Ratio of the surface water concentrations for “short” and “long” soaking times from Zouine et al. 


[1.3]; dash-dotted line: mean value, dashed lines: +standard deviation. 


The mean value was found to be Caor/Cong= 1.015 with the standard deviation SD 
=0.052. The mean value is represented in Fig. 1.5 by the dash-dotted line, the region of 
+SD by the hatched area. Thus, there is no indication for an increasing ratio of surface 
water concentration with increased soaking time as there is for water vapour. From this 
result we can conclude that the results from [1.3] represent saturation values. Conse- 
quently, a constant surface concentration of water and the use of eq.(1.2.3) is justified for 


the analysis. 


1 Diffusion in the absence of external stresses 


It should be noted that not the high water pressure is responsible for the nearly time- 
independent water concentration for the tests in [3], since in the experiments by Oehler 
and Tomozawa [6], the saturation pressure was at 350°C clearly higher than in the tests 
by Zouine et al. [3]. Obviously the direct contact of silica and liquid water is responsible 


for a large value of A/VD. 


1.3 Mass transfer coefficients 


1.3.1 Parameter h/VD from literature data 


First computations of the mass transfer parameter were reported in [1.11]. These results 
are based on total water uptake and surface water concentration. They may be addressed 


first. 


The molecular water uptake per surface unit, mc, after time ¢ results by integrating the 


water profile of eq.(1.2.5) from z=0 to z>». With the help of Mathematica [1.12] we 


m. = [coaz =C, ane pi í “E pa (1.3.1) 


For large values of AV(#/Dp) the water uptake tends asymptotically to 


obtain 


m, CD p| (1.3.2) 


The time Vh at which the asymptote coincides with the Vtaxis is 
pe 033) 


Under the assumption that equilibrium between molecular and hydroxyl water is reached 
even for short soaking times, eqs.(1.3.1) and (1.3.2) are valid also for the hydroxyl S. 
Uptake measurements by Wakabayashi and Tomozawa [1.9] for the hydroxyl water S are 
plotted in Fig. 1.6 by the solid and open circles. 
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1.3 Mass transfer coefficients 


In the first evaluation given in [1.11], we fitted eq.(1.3.1) to all the data and obtained the 
curves in Fig. 1.6a. The fitted parameters AND are plotted in Fig. 1.7 as the open circles. 
An evaluation of the abscissa sections Vo in Fig. 1.6b via eq.(1.3.3) resulted in the trian- 


gles of Fig. 1.7. 


Absorbance Absorbance 


0.01 F 


300°C 


ran 1 1 
-25 5 75 10 125 15 
vto vt (h?) 


Fig. 1.6 a) Water uptake of silica by water vapour soaking at different temperatures and water vapour 


pressure of 355 Torr after Wakabayashi and Tomozawa [1.9], curves fitted to the data via 


eq.(1.3.1), b) determination of the asymptotes from results at large times using eq.(1.3.3). 


The straight lines in this diagram yield 


h Q (1.3.4) 
JD Aca 3 


with the fitting parameters: A=527 (1/h!?), Q=38.1 (kJ/mol) for the circles and A=306 
(1/h"®), Q=34.0 (kJ/mol) for the triangles. 


The surface concentration values for molecular water obtained from the measurement by 
Oehler and Tomozawa [1.6] (solid circle) are also in good agreement with the uptake 
data. The surface hydroxyl measurements (S) by Wakabayashi and Tomozawa [1.9] 


result in a clearly lower value (square). The 90% CI are given by the bars. 
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1 Diffusion in the absence of external stresses 


All data in Fig. 1.7 apart from the square reflect directly or indirectly the behaviour of the 
molecular water. This holds of course for the solid circle since this value was determined 


directly from the surface history of molecular water. 


8 (°C) 
550 500 450 400 350 300 250 
2 
h 
\D 
1 
(1m?) 
0.5 
0.2 
0.1 
0.057 
1 N N f N 
0.0012 0.0014 0.0016 0.0018 


UT (kK?) 


Fig. 1.7 Parameter A/ VD from [1.11] as a function of temperature for constant water vapour pressure of 355 
Torr; circles: evaluation of the data in Fig. 1.6a by fitting to eq.(1.3.1), triangles: results from the 
asymptotic behaviour via eq.(1.3.3). Solid circle represents data from Fig. 1.1a (Oehler and Tomo- 
zawa [1.6], saturation pressure). Square gives result from Fig. 1.2 (Wakabayashi and Tomozawa 
[1.9], Davis and Tomozawa [1.4], 355 Torr). 


On the other hand, the asymptotic evaluation of eqs.(1.3.1) and (1.3.3) is predominantly 
governed by the long-time data. The evaluation of Fig. 1.6a includes also some data 
measured at short times for which S/C-equilibrium might not have been reached. Since 
the majority of data points are very close to the asymptotes, the asymptotic behaviour 
prevails. Consequently, the shape parameters hINDfor S(as measured in [1.9]) should be 


approximately the same as for C. 
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1.3 Mass transfer coefficients 


1.3.2 Re-evaluation and replenishment of data from literature 


Considering again the data in Fig. 1.2a, large scatter in the measurements is noticeable 
and the deviations of the individual data points from the fitting curve are enormous. This 
holds especially for soaking times >400 h and is clearly caused by the high value at 900h, 


which appears extremely questionable. 


Therefore, we excluded the 900h-value and fitted eq.(1.2.6) to the remaining data. The fit 
result was found as 4/VD=0.11 [0.077,0.148] (INh). The related curve is introduced in 
Fig. 1.8a as the solid line. The dashed curve represents again the result from Fig. 1.2. 

In the paper by Wakabayashi and Tomozawa [1.9] additional data of surface concentra- 
tions are reported. These results were obtained by stepwise removal of the surface via 


etching. 


Figure 1.8b shows results for 550°C. These data were fitted with eq.(1.2.6) yielding AND 
=0.66 [0.56, 0.75] (INh). For the other temperatures, only two data points were available 
so that the application of a fitting procedure was not possible. Table 1.1 compiles the new 


results that are also introduced in Fig. 1.9 by the red symbols, from [1.13]. 


127 O 30 
a b) 
so) | ) | sc) 
|i a OPS SS SE Se O 
107 oi en 450°C 
(em)! Lee (em) 
8 s 
e 
550°C 
6 H 
I 350°C, 355 Torr 
4 
| sog’ 
2 (0) Davis and Tomozawa | aside 
land tomoan aO Or R © 200°C 
@ Wakabayashi and Tomozawa 
1 1 1 1 l I 
200 400 600 800 200 300 


t (h) t (h) 


Fig. 1.8 a) Re-evaluation of strongly scattering data by Davis and Tomozawa [1.4] neglecting the data 
point at 900h, b) additional data by Wakabayashi and Tomozawa [1.9], evaluation from [1.13]. 
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1 Diffusion in the absence of external stresses 


8 (°C) 


500 400 300 200 
T T T T 


1.2 1.4 1.6 1.8 2 2.2 
1000/T (K*) 


Fig. 1.9 Parameter A/ \Das a function of temperature from Fig. 1.7. Additionally plotted: parameters evalu- 
ated from Figs. 1.8a and 1.8b shown as red symbols. 


The scatter of heat transfer parameters from literature may be caused by the different 
vapour pressures (saturation pressure and 355Torr), and procedures (global water uptake 


and stepwise surface removal). 


In addition, the water-silica reaction is accompanied by the generation of swelling strains 
as will be outlined in detail in Section 2. Such swelling strains were observed in [1.1] for 
a GE-124 type-I glass by measuring the change of the curvature of silica disks after hot- 


water soaking. 


Temperature (°C) Basis data AND (1/ vh) Surface concentration 
350 [1.4, 1.9] 0.11 [0.077, 0.148] 9.4 (cm!) 
550 [1.9] 0.66 [0.56, 0.75] 20.2 (cm!) 


Table 1.1 Additional mass transfer coefficients from curve fitting via eq.(1.2.6). 
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1.3 Mass transfer coefficients 


1.3.3 A possible reason for a limited water entrance 


What we have demonstrated in the preceding Sections is that by assuming that silica 
glass surfaces present a barrier to the diffusion of water, we can duplicate the experi- 
mental diffusion results of Oehler and Tomozawa [1.6] and those of Helmich and Rauch 
[1.5]. We can also show that experimental diffusion results obtained in liquid water 
should be essentially different from those obtained in water vapour. The experiments of 
Zouine et al. [1.3] demonstrated the absence of a barrier when diffusion experiments are 
done in liquid water. These authors obtained the identical surface concentrations for 
specimens held for very different lengths of time at temperatures ranging from 23°C to 
200°C. A barrier is observed only when silica is exposed to water vapour. Immersion in 
liquid water differs significantly from immersion in the vapour, because silica dissolves 
in the liquid [1.14, 1.15, 1.16], but not in the vapour [1.14]. Hence in liquid water, the 
surface is constantly being regenerated and any barrier that might form is dissolved 


away. 


By comparison, specimens exposed to water vapour do not lose weight [1.14]. Not only 
that, their surfaces are still polished after exposure [1.6, 1.1] suggesting that whatever 
enters the surface during the exposure process, remains there when the specimen is re- 


moved from the moist environment. 


Recent investigation of silica surfaces show that silica is more reactive chemically at the 
surface than in the bulk, which may account for the formation of a diffusion barrier at the 
surface. These surfaces do not dissolve away when exposed to water vapour. In contrast 
to the SiO2-structure in the bulk of silica, under-coordinated Si and O atoms occur at the 
free surface because of missing neighboured atoms. These are more reactive so that when 
this surface is exposed to water vapour, water molecules react nearly instantaneously 
with the SiO2 surface molecules to form SiOH as has been shown by molecular dynamic 


computations (see e.g. Mahadevan and Garofalini [1.17]). 


The composition of a silica surface after contact with water molecules is illustrated in 
Fig. 1.10a. It is conceivable that such a layer can show a completely different diffusion 
behavior than the bulk material. Such a difference would be necessary to form a diffusion 


barrier. 


15 


1 Diffusion in the absence of external stresses 


When the silica surface is exposed to a water vapour environment, additional water mol- 
ecules are docked at the hydroxyl groups. Water molecules are trapped at the “free” ends 
of the surface molecules for a short while resulting in a dynamic equilibrium: water mol- 
ecules coming in from the environment are then captured by the surface, while water 
released from the surface is free to diffuse deeper into the bulk as has been studied by 


computer modelling of the diffusion process [1.18,1.19,1.20]. 


SiO2 bulk 


Surface layer panna N 
Ernai 


e N w Ff 


~ e b) 


Fig. 1.10 a) Ata free silica surface water reacts nearly instantaneously with under-coordinated silica species 
in a thin surface layer of about 0.7 nm thickness; per 2-coord O (and 3-coord Si) 1 H2O molecule 
reacts to SiOH (curves according to Mahadevan and Garofalini [1.17]), b) surface layer of water 


molecules on a silica surface (pink: H of SiOH or H20, grey: O of silica; blue: Si). 
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1.3 Mass transfer coefficients 


A surface layer formed by this process, Fig. 1.10b, was developed by Mahadevan and 
Garofalini [1.17]. In Fig. 1.11a the hydroxyl concentration is given in weight percents 
(right ordinate). Unfortunately, the swelling behaviour of the under-coordinated Si and O 


atoms of silica is not known. 


The measurements by Briickner, Shackelford and Shelby predominantly represent the 
behaviour of the bulk material with “correctly” coordinated oxygen and silicon. There- 
fore, we tentatively apply the dependency of volume swelling strain as a function of 


hydroxyl concentration, &,(5), given in [1.21] as an approximation. 


The hydrostatic swelling stresses can then be computed as shown in [1.21]. The result is 
represented by the left ordinate of Fig. 1.1la. Compressive hydrostatic stress reaches a 


maximum value of oswn = —2.8 GPa. 


108 F 
D/Do 


Fig. 1.11 Hydroxyl concentration (a) in weight percent, resulting in compressive hydrostatic swelling stress- 
es (upper curve). The lower curve (b) shows the strongly reduced water diffusivity D, normalized 


on the diffusivity D» for stress-free silica at room temperature. 
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1 Diffusion in the absence of external stresses 


Due to such high compressive stresses, the water diffusivity must be strongly reduced 
within this thin layer. When swelling stresses are present, the diffusivity is given by the 


following equation: 


AV, 
D=D,exp| o,,,— (1.3.5) 
0 il sw,h RT | 


with the effective volume A W= 16 cm?/mol as derived by Lezzi et al. [1.22]. 


Figure 1.12 schematically illustrates the consequences of the reduced diffusivity in the 
surface layer. In Fig. 1.12a the distribution of the diffusivity along the depth coordinate is 
plotted once more. Figure 1.12b shows the water profile for a fixed time. The water con- 
centration directly at the surface, Co, is assumed generally to be constant only depending 


on the partial pressure of water in the environment. 


Silica 


Fig. 1.12 a) Schematic representation of the diffusivity in the surface region, b) related water concentration 


profile, c) finite surface layer replaced by a surface with a finite mass transfer coefficient. 
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1.4 Diffusion in thin specimens 


From experimental water concentration measurements as for instance performed by Oeh- 
ler and Tomozawa [1.6], the distribution of Cis obtained. Extrapolation of the measured 
curves to 20 results in the dashed extension through the surface layer, Fig. 1.12b. The 
surface value of the extrapolated profile, C(O), (indicated in Fig. 1.12b by the circle) must 


of course be smaller than C and must increase with time as is found in the experiments. 


The suppressed water transport within 0<z<d can be taken into consideration by replac- 
ing the thin layer by the surface (d>0) and introducing a finite mass transfer coefficient, 


h, into the surface boundary condition (Fig. 1.12c). 


When silica is exposed to Liquid water, dissolution of the glass into the water has to be 
expected. Dissolution of silica in water has been reported by a number of authors and is 


summarized in [1.23]. 


From these data we can conclude that in liquid water the silica surface in Fig. 1.10b and 
the highly stressed layer in Fig. 1.11 are continuously removed by dissolution. By con- 
trast, in water vapour the transport of the SiOH from the surface is not possible. This is 
probably the reason for the difference in diffusion behaviour in liquid water and in water 


vapour environments. 


1.4 Diffusion in thin specimens 


1.4.1 Analytical solutions from literature 


In the experimental and theoretical considerations in Section 1.2, diffusion in thick speci- 
mens was discussed. The diffusion zones were negligible compared to the thickness of 
the bulk material, /Dxt <<W, so that the limit case of a half-infinite body is sufficient- 
ly fulfilled. 


For thin specimens of width 2 W, the zone thickness may be comparable with the speci- 
men thickness. Then the concentration distribution as a function of time reads according 
to [1.8] 


C(z,t)=C,1 5 Te er a = | _ Weh/p (141) 
par Wh'(Wh'+1) + æ; W 
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1 Diffusion in the absence of external stresses 


with the time-dependent surface concentration (z=0) 


2 2Wh' ,D 
C(t) =C,]} 1 ex a t (1.4.2) 
‘© Í 2, Wh'(Wh's1) +a 1 "W? ) 


where G is the saturation value and om, n=1,2,... are the positive roots of 
atana =Wh'=Wh/D (1.4.3) 


In eqs.(1.4.1-1.4.3) the parameter Wh/D is equivalent to the Biot number B (sometimes 


also written Bi) for thermal problems: 


wha “hp (1.4.4) 
D 
Then the equations (1.4.1-1.4.3) read 
čchi x 2Bcos[a,(z/W - Disceta) csp =i? 2 (1.4.1a) 
al B(B+1)+ a; B 
C (=G 1 25 5 exp] a a (1.4.2a) 
1 B(B+1)+a; B 

where now On, 2= 1,2, ... are the positive roots of 

atana=B (1.4.3a) 


The water concentration profiles as a function of time are shown for B=1 in Fig. 1.13a in 
normalized representation. After a normalized time of t=0.03, the diffusion fronts from 


the opposite surfaces start to interfere noticeably. 


The water concentration at the surface, C(0), is given in Fig. 1.13b for several Biot num- 
bers. At first sight, these results may appear somewhat astonishing since for increasing 
Biot numbers the surface concentrations at a fixed time t decrease. For constant D and 
W, this says that Cp decreases with increasing mass transfer coefficient A. In this context, 
it has to be noticed that the normalized time t also includes the parameter A via tcf’, so 
that for short normalized times, the curves become independent of A. This is visible from 
Fig. 1.13b. 
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1.4 Diffusion in thin specimens 
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Fig. 1.13 a) Concentration profiles for Dx W° = 1/B’= 1 and different normalized times t, b) normalized 
surface concentrations C/C) vs normalized time for different Biot numbers, c) results with modi- 


fied abscissa scaling t/B’ for small Biot numbers. 
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1 Diffusion in the absence of external stresses 


Figure 1.13c shows results of C(O) for small values of B with the abscissa t/B instead of 


t. These curves now show the properties to be expected intuitively. 


Since the water uptake mic is proportional to the average water concentration over the 


cross-section, it holds simply 


u (1.4.5) 
with the average water concentration C : 
_ 2W 
C(O = |C, ddz (1.4.6) 
0 


and the saturation uptake zmo. 


The molecular water uptake as a function of normalized time t and Biot number is given 
in Fig. 1.14a and Fig. 1.14b for Vr and Vt/B chosen as the abscissa. As a characteristic 
feature of the effect of a limited water transfer coefficient A, the curves do not start line- 
arly. This would only occur under the assumption that the surface concentration of mo- 


lecular water would be constant during the water-soaking test, i.e. for C(0,=Co=const. 
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1.4 Diffusion in thin specimens 
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Fig. 1.14 a) Water uptake mic/mm as a function of normalized time t and Biot number B, b) the same results 


with different abscissa, c) measurements by Davis and Tomozawa [1.4]. 


In real tests (see Fig. 1.14c), the surface concentrations increase for short times oVtand 
the layer thickness goes with V£ too. This makes that the uptake is proportional to time. 
In a plot with Vtas the abscissa, the uptake must start with a quadratic dependency. This 


effect is also visible from literature results by Davis and Tomozawa [1.4]. 


1.4.2 Results for high temperatures >500°C 


The equilibrium constant of the reaction (1.1.1) is at high temperatures, 7>500°C, repre- 


sented by the ratio 


2 
ko St (1.4.7) 


Results of Ş as a function of swelling stresses are given in Figs. 1.15a and 1.15b. The 
influence of Biot number is shown in Fig. 1.15c. In principle the same behaviour is visi- 
ble as in the case of lower temperatures. In Figure 1.15 the asymptotically reached max- 


imum hydroxyl content is denoted as 5 given by 


Ss S4/ ko Go (1.4.8) 


23 


1 Diffusion in the absence of external stresses 
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Fig. 1.15 Hydroxyl concentration at the surface, a) and b) effect of swelling stress, c) effect of the Biot num- 


ber. 
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1.5 Water uptake 


1.5 Water uptake 


The uptake of molecular water, mic, is obtained from eqs.(1.4.1a) with C (t) from (1.4.6) 


RR <2 2B? aT | 
m- =C(t)x2W =m, | 1 ` (1.5.1) 
eH OO ol = aaa On ) 


where the asymptotically reached value of molecular water uptake is denoted as 


Me) = 2WC, (1.5.2) 


Equation (1.5.1) holds in the absence of stresses. In this case the hydroxy] content is 
m, (t) = k, ome (t) (1.5.3) 


In presence of swelling stresses the equilibrium constant is 


0,AV 
k=k, ex] Er l (1.5.4) 


with the hydrostatic stress in first order 


Bee (1.5.5) 

91-v) ~ 
An approximate analytical solution can be obtained following the general procedure 
usual in perturbation theory. If we consider the equilibrium constant kı as the disturbance 
parameter, perturbation theory suggests to solve the problem in the absence of stresses, 
i.e. for ki=kio =const. and to insert the stress-dependent disturbance parameter into this 
solution. The analytical solution of the undisturbed problem in terms of molecular 


water is QZ, ð given by eq.(1.4.1-1.4.4). 


The first-order solution for the hydroxyl water uptake under swelling stresses then reads 


m,() =k,, | CC, De 2 une de (1.5.6) 


0 


with the abbreviation A defined by eq.(3.1.14). 
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1 Diffusion in the absence of external stresses 


The effect of the disappearing stresses due to the exponential term C(z,t)—C,—0 for 


>>» results in a reduction of the hydroxy] uptake as is shown in [1.24]. 


References 


1.1 S. M. Wiederhorn, F. Yi, D. LaVan, T. Fett, M.J. Hoffmann, Volume Expansion 
caused by Water Penetration into Silica Glass, J. Am. Ceram. Soc. 98 (2015), 78-87. 


1.2 Doremus, R.H., Diffusion of water in silica glass, J. Mater. Res., 10 (1995), 2379- 
2389. 


1.3 A. Zouine, O. Dersch, G. Walter and F. Rauch, “Diffusivity and solubility of water in 
silica glass in the temperature range 23-200°C,” Phys. Chem. Glass: Eur. J. Glass Sci and 
Tech. Pt. B, 48 [2] (2007), 85-91. 


1.4 Davis, K.M., Tomozawa, M., Water diffusion into silica glass: structural changes in 
silica glass and their effect on water solubility and diffusivity, J. Non-Cryst. Sol. 185 
(1995), 203-220. 


1.5 M. Helmich and F. Rauch, “On the mechanism of diffusion of water in silica glass,” 
Glastech. Ber. 66 [8] (1993), 195-200. 


1.6 Oehler, A., Tomozawa, M., Water diffusion into silica glass at a low temperature 
under high water vapor pressure, J. Non-Cryst. Sol. 347 (2004) 211-219. 


1.7 R.H. Doremus, Diffusion of Reactive Molecules in Solids and Melts, Wiley, 2002, 
NewYork. 


1.8 Carslaw, H.S., Jaeger, J.C. Conduction of heat in solids, 2nd ed. 1959, Oxford Press, 
London. 


1.9 Wakabayashi, H., Tomozawa, M., Diffusion of water into silica glass at low tempera- 
ture, J. Am. Ceram. Soc. 72 (1989), 1850-55. 


1.10 H. Le Chatelier, C.R. Acad. Sci. Paris 99(1884), 786 


1.11 S.M. Wiederhorn, G. Rizzi, S. Wagner, M.J. Hoffmann, T. Fett, Diffusion of water 
in silica glass in the absence of stresses, J. Am. Ceram. Soc. 100 (2017), 3895-3902. 


1.12 Mathematica, Wolfram Research, Champaign, USA. 


26 


1.13 G. Schell, T. Fett, Mass transfer coefficients for water at silica surfaces - Additional 
results, Scientific Working Paper SWP 73, 2017, ISSN: 2194-1629, Karlsruhe, KIT. 


1.14 H. Li and M. Tomozawa, Mechanical strength increase of abraded silica glass by 
high pressure water vapor treatment, J. Non-Cryst. Sol, 168 287-291 (1994). 


1.15 G. Perera, R.H. Doremus, J. Am. Ceram. Soc. 74 (1991), 1269-74. 


1.16 P.M. Dove, N. Han, A.F. Wallace, J.J. De Yoreo, Kinetics of amorphous silica dis- 
solution and the paradox of the silica polymorphs, PNAS, 105, [29] (2008), 9903-9908. 


1.17 T.S. Mahadevan, S.H. Garofalini, Dissociative chemisorption of water onto silica 
surfaces and formation of hydronium ions, J. Phys. Chem. C 112(2008), 1507-1515. 


1.18 Feuston, B.P., Garofalini, S.H., Water induced relaxation of the vitreous silica sur- 
face, J. Appl. Phys. 68(1990), 4830-4836. 


1.19 Leed, E.A., Pantano, C.G., Computer modelling of water adsorption on silica and 
silicate glass fracture surfaces, J. Non-Cryst. Sol. 325(2003), 48-60. 


1.20 Ma, Y., Foster, A.S., Nieminen, R.M., Reactions and clustering of water with silica 
surface, J. Chem. Phys. 122 144709 (2005). 


1.21 S. M. Wiederhorn, G. Rizzi, S. Wagner, M. J. Hoffmann, and T. Fett, Stress- 
Enhanced Swelling of Silica: Effect on Strength, J. Am. Ceram. Soc. 99 [9] 2956-2963 
(2016) 


1.22 P.J. Lezzi, Q.R. Xiao, M. Tomozawa, T.A. Blanchet, C.R. Kurkjian, Strength in- 
crease of silica glass fibers by surface stress relaxation: A new mechanical strengthening 
method, J. of Non-Crystalline Solids 379 (2013) 95-106. 


1.23. N.P. Bansal and R.H. Doremus, Handbook of Glass Properties, Academic Press, 
Inc. Orlando, FL (1986), pp. 25. 


1.24 T. Fett, G. Schell, Water diffusion in thin silica sheets, Effect of swelling stress on 
diffusivity, Scientific Working Papers 72, 2017, ISSN: 2194-1629, Karlsruhe, KIT. 


27 


2 Swelling in the absence of external 
stresses 


2.1 Swelling and swelling stresses 


2.1.1 Qualitative indications in literature 


Clear evidence has been reported in the literature for volume swelling, &, measured di- 
rectly via dimension changes and indirectly via stress generation within the silica glass. 
The swelling behavior has been determined by length [2.1] and curvature measurements 
on glass bars [2.2] and by X-ray diffraction strain measurements [2.3] on silica surfaces. 
Gorbacheva and Zaoints [2.1] soaked fused silica prisms (25x5x5 mm) in water at 80 °C 
for 20 months. They found that the measured length had increased by 0.17% during this 
treatment. Since only a thin surface layer of water penetrated glass, ~ 5 um, could have 
been generated during the water storage, the suppressed linear strain in the surface must 


be clearly larger than 0.17%. 


2.1.2 Swelling strains from density measurements 


In the past, the authors showed in a couple of papers the principle effects of volume 
swelling by hydroxyl generation in silica on mechanics and fracture mechanics properties 
[2.4, 2.5, 2.6]. 


A swelling effect in water-containing silica at high temperatures was early reported by 
Briickner [2.7, 2.8], Shackelford [2.9] and Shelby [2.10]. These authors showed that the 


density decreased by reaction (1.1.1) for a fictive temperature of 77=1100°C according to 


Ap _ go m q m (2.1.1) 
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2 Swelling in the absence of external stresses 


where Cy is the weight fraction of water. The literature data for a fictive temperature of 


T=1100°C are shown in Fig. 2.1 by the open symbols. The individual data for Brückner 
and Shelby were taken from the compilation by Shelby [2.10]. The data of Shackelford 
for fictive temperatures 7i=1000°C, 7i=1100°C and 7i=1200°C were directly taken from 
the original paper [2.9] transforming the “water” content (wt% OH) in water content 


(wt% H20) using the factor (18/17)/2 where 18/17 reflects the different molar weights of 


OH and H20. 
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Effect of water concentration on density of vitreous silica, results by Shelby [2.10], Brückner [2.7, 


2.8], and Shackelford [2.9]. The dash-dotted line represents the value given in (2.1.2). The open 
symbols are for a fictive temperature of 77=1100°C, the black symbols for 7i=1000°C and the red 
ones for 7;=1200°C. The circles represent results by Shelby [2.10] with the water produced in sili- 


ca saturated with hydrogen and then gamma-irradiated. 


By a least-squares fit of the data points for 7i=1100°C and 1000° exclusively, we obtain 


x = 9.836 [0.739, 0.933] 
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(2.1.2) 


2.1 Swelling and swelling stresses 


with the 95%-Confidence Interval in brackets. The dependency given by this value is 
indicated in Fig. 2.1 as the dash-dotted line. In our former papers [2.11, 2.4-2.6] we used 


the value of y=0.84 as was obtained simply as the average of the 3 values given in [2.10]. 


From the definition of the density as the quotient of mass Mr Mauasstm, and volume V, 


p=m/V, it follows for the volume swelling strain €v 


_AV _Am_ Ap (2.1.3) 


vV-” >e, 
P V m pP 


Since only the water content m2 can change during water soaking, the mass change is 


Am=m,, =C,m (2.1.4) 


for Mw<<sMelass. Consequently, the volume swelling strain ev simply results from 
eqs.(2.1.1-2.1.4) as 


é,=(1+ x)C, (2.1.5a) 


In the following, let us consider the case of swelling behavior for both the =SiOH and the 
molecular water that is present at low temperatures <500°C. The total swelling strain ev 


is then composed of the individual contributions of the hydroxyl, £, and the molecular 


water, £c, according to a molar “rule of mixture”: 


Pr ap (2.1.5b) 
C+S/2 C+S/2 


So far, no reliable information on volume change is available for the strain ec by the 
molecular water species. We believe that the molecular water in the glass is located in 
holes of molecular size in the silica glass network [2.12, 2.11], where the water may not 
contribute to the volume change, i.e. we assume in the following considerations that 
&c=0. On the other hand, it is well known that molecular water at high temperatures of 
1100°C is nearly absent. Under these circumstances, the volume swelling in [2.7-2.10] is 
due to hydroxyl generation by the reaction (1.1.1), i.e. &,=es. This may be confirmed and 
if necessary improved by future experiments of Molecular Dynamic (MD) computations. 


For this topic see also the Outlook-Section. 
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2.1.3 Transformation into hydroxyl concentrations and 


swelling strains 


In general Cy has to be converted to mole fraction of water, Xpo , in the glass via 


M 
Xino = TSG = (60/18)C,, (2.1.6) 


w 


Here we used the molar mass for water my =18g/mol and as an approximation Metass= 
60g/mol for the molecular weight of the glass. For every water molecule that reacts with 


the glass, two SiOH are formed, therefore: 


Kar = 2 Xp (2.1.7) 
Combining (2.1.6) and (2.1.7) yields 
Nine 2=C, (2.1.8) 


By assumption, the volume expansion is just due to the =SiOH formation. The hydroxy] 


mass concentration Sis 


Ms; 17 
S = Kon <= Xsion 60 (2.1.9) 


glass 
Consequently, the volume swelling strain in terms of the water concentration Cw is 
e, =(1+ y )C, =1.836xC,, [1.739xC,,,1.933xC,,] (2.1.10) 


Here it should be emphasized once more that this dependency holds for sma// water con- 


centrations. This relation reads in terms of hydroxy] mass concentration $ [2.13]: 


S 
&= 7+ 7) =exS (2.1.11) 
Then 


x =0.97 [0.92,1.02] (2.1.12) 
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In (2.1.12) the numbers in brackets represent the 95% confidence interval. The value of 
x=0.97 is recommended to be used for “best” prediction, the lower boundary value of the 
CI, «=0.92, for “conservative” predictions. Our coefficient used so far [2.11], «=0.92, is 


located at the lower boundary of the 95%-CI. 


The data by Shackelford for 7:=1000°C represented by the black solid symbols fit very 
well to the data for 7i=1100°C. Therefore, we included these data in the analysis with the 


result of 
«x =0.965 [0.916,1.014] (2.1.13) 


that is practically identical with the value of (2.1.12). Including all data for the three 


fictive temperatures gives 
K=0.92 [0.86,0.984] (2.1.13a) 


To our actual knowledge we suggest use of k=0.97. Having in mind the limited accuracy 
of water concentrations (in the Shackelford analysis +75 ppm H20) the equation for wa- 
ter concentrations in terms of mass units seems to be sufficiently accurate also for molar 


units because the mole masses deviate not strongly, i.e. 


m, _ 18 g/mol _ | 058=1 (2.1.14) 
Mson 17 g/mol 


2.1.4 Results by Shelby on hydrogen infiltrated and gamma- 


irradiated silica 


Shelby [2.10] generated water in the glass on a further way. He saturated the glass with 
hydrogen H2 under gas pressures up to 70 kPa. The specimens were then gamma- 
irradiated to doses of 1 to 15 Grads. In this way, SiOH was directly produced in the bulk 
material. The resulting density reduction by this procedure is shown in Fig. 2.1 by the 
solid circles. These results do strongly deviate by roughly a factor of two from all the 
other data. The regression analysis according to [2.14] for these data gives 

—=-y7C x =1.735 [1.60,1.87] (2.1.15) 


w? 
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For the computation of the volume strain it has to be taken into account that the resulting 
water has not been introduced from the environment, but was generated in situ with the 


oxygens taken from the SiOz network. 


Unfortunately it is not clearly visible in [2.10] whether the initial density was measured 
before or after saturation by hydrogen. Therefore, we will perform two different 


evaluations as performed in [2.15]: 


Assumption 1: The initial density was measured before hydrogen infiltration. Conse- 


quently, the mass change is 2/18 Cy since 2 hydrogen atoms are needed for each water 


molecule. 
Am=1C,m (2.1.16) 
In this case, the volume strain results as 
&,=6+r)C, (2.2.17) 
and according to (2.1.11) we obtain 
e =BG) =s (2.1.18) 


Introducing the result of eq.(2.1.15) yields 


x =0.977 [0.906,1.049] (2.1.19) 


Assumption 2: The initial density was measured after the hydrogen infiltration. In this 


case it is simply Am=0 


er (2.1.20) 
From this it yields 
K =0.92 [0.847,0.99] (2.1.21) 


Within the 95%-confidence intervals, also the results of (2.1.19) and (2.1.21) are in rea- 


sonable agreement with the results from heat-treated specimens. 
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2.1 Swelling and swelling stresses 


Finally, the volume increase ev is plotted in Fig. 2.2a as a function of the water concen- 
tration Cw and in Fig. 2.2b versus the hydroxyl concentration .S for all measurements in 
Fig. 2.1. The symbols in Fig. 2.2 are the same as given in Fig. 2.1. The straight line, 
introduced in Fig. 2.2b, represents the coefficient k=0.97 according to (2.1.13). The good 


agreement with all the data is evident. 
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Fig. 2.2 a) Volume expansion strain vs. water concentration, b) expansion vs. hydroxyl concentration S, 
(symbols as in Fig. 2.1), line: eqs.(2.1.20) and (2.1.21), from [2.14, 2.15]. 


2.1.5 Swelling stresses under isotropic swelling 


The effect of the water-silica reaction is the generation of swelling strains. Swelling 
stresses are a consequence of the mechanical boundary conditions. A volume element in 
a thick plate that undergoes swelling cannot freely expand. If the diffusion zone is small 
compared to the component dimensions, expansion is completely prevented in the plane 


of the surface and can only take place normal to the surface plane. 


At a free surface, the stress state is plane stress and, consequently, also stresses caused by 


swelling are equi-biaxial (0x=0) 
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RER EE (2.1.22) 
y : 3(1-v) 


In (2.1.22) Eis Young’s modulus and v is Poisson’s ratio. Consequently, the hydrostatic 


stress term reads 


2E, E 
O sw,h =3 (Owy Oow) zu E (2.1.23) 
i ' i 9(1-v) 


Introducing the proportionality (2.1.11) in (2.1.22) gives 


saaa E E TOS (2.1.24) 
í 9(1-v) 


2.2 Swelling from measurement of disk curvature 


2.2.1 Evaluation of curvature measurements 


Swelling stresses were derived in [2.11] by measuring the change of the curvature of 
silica disks after hot-water soaking. At the heart of this study is an experiment in which 
water is diffused into one side of thin vitreous silica disks. This technique has been used 
earlier to determine the volume change caused by the penetration of water into silica 
coatings on silicon wafers [2.2]. The consequent volume expansion at the glass surface 
results in a bending of the disks, which are no longer flat after exposure, but concave- 
shaped [2.16]. By measuring the radius of curvature at the surfaces of the curved disks 
and by understanding the mechanism of water migration, the surface stress can be calcu- 


lated. The technique can also be used to calculate the diffusivity of water in silica glass. 


Two different experimental procedures were applied in [2.11] for tests at 90°C and 200°C. 


90°C: Silica disks of W=250um thickness were coated on one side by e-beam evaporation 
with a 3 nm bond-coat of titanium, followed by 10 nm of gold. Then the initial curvature, 
Ro, was measured using the wafer-curvature, laser-scanning technique. Disks were ex- 


posed to water at 90 °C for periods ranging from two to six weeks. One specimen series 
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was immersed in water, a second test series suspended over the water — the uncoated side 


facing the water, the coated side facing the relatively dry furnace air. 


A plate is considered with all dimensions very much larger than the water diffusion zone. 


The distribution of the swelling stresses is given by 


fe = overt =) (2.2.1) 
2b 
with b=./Dt (2.2.2) 


The bending moment caused by a thin swelling layer is then given as 


(2.2.3) 


0 


M, =o, ferfe Z dz=Wo, 2 
2 2b Vn 
RE Ae 


2b 


lee 


The measurable change of the curvature by swelling is (see e.g. Timoshenko and 


Goodier, Theory of Elasticity, page 290) 


pe ey (2.2.4) 
R EW? 

In [2.11] the diffusion parameter b was computed with the diffusivity given by Zouine et 

al. [2.17]. Using the surface water concentration reported by Zouine et al. [2.17] together 

with the S/Cratio from IR-measurements in [2.11], the surface swelling stress oo was 


obtained via 


Ge KE o (2.2.5) 
3(1-v) 


The bending moments computed via eq.(2.2.3) are compared in Fig. 2.3 with the experi- 
mental results obtained from (2.2.4). Figure 2.3 shows the measured and computed mo- 


ments plotted versus square root of time. 


The prediction made with the surface concentration from [2.17] is entered as the triangle. 


In a second attempt the total water uptake mcw was used as measured by IR as well as the 
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S/C ratio also obtained from IR. As shown in [2.11] it holds for the generated bending 


moment expressed in terms of the hydroxyl concentration S 


-pW E .m, (2.2.6) 
2 3(-v) 


b= 


where zns is the uptake of the hydroxy] water from one side. 
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Fig. 2.3 Measured bending moments (via disk curvature) versus time (circles) for 90°C. The predicted 


moment obtained from /R data is given by the square. The triangle is determined via the water 


concentration data from Zouine et al. [2.17]. 


The result from this equation is plotted in Fig. 2.3 as the square. The prediction using the 
surface water concentration by Zouine et al. [2.17] (triangle) is in very good agreement 
with the curvature measurements. Even though the square lies just outside of the confi- 
dence limits, these limits are for the mean value of Ms and not the individual data points. 


A certain fraction of the data should lie outside of the confidence band. The two predic- 
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tions indicate support for the assumption that only the =SiOH species contribute to the 


volume expansion during soaking. 


200°C: At this temperature it was not possible to prevent water entrance sufficiently by a 
gold layer. Therefore the bare silica disk was heat-treated on both sides and then stepwise 
etched from one side in a buffered HF solution. After each etching step the change of 
curvature was determined resulting in the released bending moment versus total removed 


surface layer as plotted in Fig. 2.4a for 2h and 20h soaking times at 196°C. 


In order to avoid the complicated equations in integral form used in [2.11], the procedure 


may be explained here alternatively in differential form. 


The bending moment per unit length caused by a thin surface layer of stresses o in a plate 


of thickness Wis given by 
M, = foed (2.2.7) 
0 


If the surface is removed by an amount of d, the remaining moment after a total surface 


removal (single etching depths accumulated) is 
M,(d) = food (2.2.8) 
d 
and the change of moment 


AM,(d) = z food (2.2.9) 


The left-hand side of (2.2.9) is known from measurements of curvature by using 


AM, (d) = Ew’ f u > (2.2.10) 
120-v)|R R, 


where Ris the actual and =o the initial curvature. 
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Fig. 2.4 a) Change of the bending moment during stepwise etching one surface, b) swelling stresses. 


Equation (2.2.9) is an integral equation with respect to the unknown stress distribution 
o(Z). Since the integrand does not depend on d explicitly, its solution can simply be ob- 


tained by taking the derivative with respect to the removed layer thickness d 


ol = 204M, (2.2.11) 
zd W ôd 


This equation holds for the case that the removed layer-thickness is small compared with 
the thickness of the disk, i.e. for d<< W. The measured change of the bending moment is 
plotted in Fig. 2.4a. In the procedure described in [2.11] the swelling stresses were first 
represented by appropriate erfc-type functions with unknown parameters and then inte- 
grated according to eq.(2.2.7). The resulting bending moments were then fitted to the 
measured data resulting in the best set of parameters for the stress distribution o(2). In the 
procedure based on eq.(2.2.11) the measured data were transformed into stresses and 
these were then “smoothed” by fitting with erfc-type functions. For taking the derivative, 
the measured data were interpolated by cubic splines using the Mathematica subroutine 
Interpolation [2.18]. Then the splines were differentiated with the result shown in Fig. 
2.4b. 
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Fig. 2.5 a) Comparison of measurements (circles) and predictions (squares) based on measured IR-data and 
data reported by Zouine et al. [2.17], and Helmich and Rauch [2.19]. The dashed lines represent 
the 95% error limits of measurements; b) Ratio of measured and predicted bending moments, bars 
represent +1 SD, circles: predicted with «=0.97, eq.(2.1.12), (1) based on IR-measurements by 
Wiederhorn et al. [2.11], (2) based on data by Zouine et al. [2.17]. 


Fig. 2.4b gives the stress distribution obtained from the curves of Fig. 2.4a. Finally, Fig. 
2.5a represents the total bending moment released by removing the swelling layer com- 
pletely. This value is directly visible for the case of 2h soaking time in Fig. 2.4a. In cases 
where a plateau value is not yet visible in Fig. 2.4a, the value can be computed by 


eq.(2.2.7) using the erfc-type functions from the fitting procedure as the integrand. 


All predicted data in Figs. 2.3 and 2.5a (squares and triangle) are within the 95% error 
limits supporting again the assumption that only the =SiOH contributes to the volume 
expansion. This can be seen from Fig. 2.5b. In this figure, the ratio of the measured to the 
predicted bending moments is plotted. The symbols represent the predictions made in 
[2.11], the bars represent +1 Standard Deviation. The results of Fig. 2.5b are compiled 
also in Table 2.1. (1) gives the bending moment for water-soaked silica disks at 90°C by 
Wiederhorn et al. [2.11] together with the prediction based on the water content obtained 


via IR-measurement. The numbers in brackets represent +1 standard deviation (S.D.). 
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In (2) the same experimental data are compared with the predictions 


centrations at the surface as were reported by Zouine et al. [2.17]. 


using water con- 


(1) 


(2) 


(3) 


Fig. 2.6 


Moment (Nm/m) 
-Mh 90°C) 


-M (90°C) [2.11] 


-M (200°C) [2.11] 


predicted 
0.0047 


0.0040 


0.0145 


measured 
0.00413 
[0.00363, 0.00463] 


0.0173 
[0.0128, 0.0217] 


ratio 
0.878 
[0.772, 0.985] 
1.032 
[0.908, 1.16] 
1.19 
[0.882, 1.496] 


Table 2.1 Comparison of predicted and measured bending moments from [2.11]. 


-Osw 


(MPa) 


40 


30 


20 


10 


201°C 
e 


> 
ok 


10 12 
vt (ht?) 


Maximum swelling stresses, results from [2.11] for temperatures in the range of 188°C<0<216°C. 


Finally, Fig. 2.6 summarizes the surface swelling stresses from [2.11] plotted as a func- 


tion of heat-treatment time. In context with the data scatter, it should be mentioned that 


the individual data points represent slightly different soaking temperatures as indicated at 


the individual data points. 
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2.2.2 Conclusions from curvature measurements 


In literature a different interpretation of the mechanical consequences of the silica-water 
reaction (1.1.1) is actually discussed by Tomozawa and co-workers [2.20], [2.21]. In 
these papers it is suggested that the mechanical response on water entrance would be an 
instantaneous and complete relaxation of stresses pre-existing at a silica surface. This 
hypothetical mechanism is claimed to occur even at 200°C. This makes it possible to 
have a look on such a mechanism with regard to the bending tests at this temperature 
[2.11]. 


Our conclusions from the point of the results on disks are: 


° Before the tests, the silica disks were etched in 5% volume fraction HF in order to 
remove 5-10um of glass, which might have contained a diffusion layer of water at 
the surface or might have residual stresses due to polishing. Consequently, stress re- 
laxation could not be the reason for curvature since no stresses to relax were pre- 


sent. 


° We have to conclude that the mechanical effect of the reaction (1.1.1) is an active 
process that generates strains and depending on the mechanical boundary conditions 
stresses. It is not a passive effect like stress relaxation that needs pre-existing stress- 
es. 

e — On the other hand, Fig. 2.5b clearly shows that the predictions and measurements of 
curvature (or bending moments) are within 0.85-1.20 of the ideal agreement (dash- 
dotted line =1). We have to state that there is no free space for any further stress- 
generating effect. This confirms the assumption that the water species causing 


swelling strains is exclusively the hydroxyl Seven at rather low temperatures. 


In contrast to a suggested relaxation effect, it has furthermore to be emphasized that the 
swelling mechanism is already proved by the density reduction as the consequence of an 


increase of the water content [2.8-2.10]. 


Moreover that relaxation is not present at 200°C or at least negligible. If relaxation would 
be present at 200°C, the generated stresses in the disks causing the observed curvature 
must have relaxed instantaneously and of course not any curvature could result as ob- 


served. 
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2.2.3 Equilibrium constant 


The equilibrium constant ķı according to eq.(1.1.2) is plotted in Fig. 2.7a. The solid cir- 
cles represent the data obtained in [2.11] and the open ones are results from literature 
published by Zouine [2.17], Davis and Tomozawa [2.22, 2.23], and Oehler and Tomo- 
zawa [2.24]. The data follow an Arrhenius curve, śe., a straight line in Fig. 2.7a that is 
expressed by eq.(1.1.3). The activation energy of the least squares fit line is Q=10.8#1.1 
kJ/mol. The parabolic curves on each side of the straight line are the 95% confidence 


limits of the fit. 


2.3 Observation on etching rates 


In the etching procedure in hydrofluoric acid/sulfuric acid solutions at room temperature; 
2 N NaOH at 80 °C, we found as an interesting detail that the etching rate decreased with 
increasing swelling stresses. Results are shown in Fig. 2.7b plotted versus the swelling 
stress for disks that were soaked for 2 h at about 200°C under saturation water vapour 
pressure of 1550 kPa. The open circles show results at 216°C, the solid circles stand for a 
test at 196°C. 


From Fig. 2.7b we found that the etching rate decreased with increasing compressive 


swelling stresses. A fit over all data according to 


=a, Expla, y] (2.3.1) 


(dô is the removed layer thickness during the time interval dd resulted in the parameters 
ao=0.1179 um/min [0.1154, 0.1203], 2=0.0032 MPa’! [0.0024, 0.0041] 


with the 90% Confidence Intervals in brackets. Since the hydrostatic swelling stress at 
the surface is 2/3rd of the equi-biaxial stress, a representation via an Arrhenius relation 


gives the stress dependence of the etching rate by 


ze a, exp Fonk (2.3.2) 
dt RT 
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with an activation volume V, 
V=11.7 cm?/mol [8.8, 15] 


again with 90% Confidence Interval. 
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Fig. 2.7 a) Plot of S/Cas a function of 1/T. The solid circles are from [2.11] via IR measurements on disks 
exposed to saturated steam at 200 °C and 90 °C. The open circles are IR measurements made by 
other authors, summarized in Zouine er a/.[2.17], b) etching rate as a function of the average stress 


in the etching intervals for silica disks hot water soaked at about 200°C for 20 h. 


In contrast to the Charles/Hillig theory Tomozawa and co-workers [2.25, 2.26] report 
negative activation volumes. This means that under compressive stresses the corrosion 
rates would be larger than under tension. This is not necessarily in contradiction with the 
results in Fig. 2.7b. It may be possible that the etching rate in completely stress-free spec- 
imens containing high hydroxyl concentrations would strongly decrease and overcom- 
pensate a possible inverse stress effect. Soaking tests after very long times would result 
in a homogeneously distributed S concentration over the whole specimen with disappear- 
ing swelling stresses. The only conclusion from the etching rate measurements is so far 


that the higher the S-concentration is, the smaller are the etching rates. 
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2.4 Volume swelling and shrinking in other glasses 


In the previous Sections, only silica was considered. In this Section, we will look on 
other glasses, too. Density measurements on several glass compositions were reported by 
Scholze [2.27]. The materials were: (B) B>O;-glass, (K) K2O-SiO2-glass (20-80 mol-%), 
(Ca) Na20-CaO-SiO2-glass (16-10-74 mol-%), (Al) LizO-Al203-SiO2-glass (20-5-75 mol- 
%), and (Si) silica. 


The relative density changes due to water content Cw are given in Fig. 2.8a. Whereas for 
silica a non-linear behaviour Ap( Gy) is reported in [2.27], for all other glasses straight- 
line behaviour was found. It has to be taken into account that the data scatter for silica is 
rather large as the data reported by Shelby [2.10] show which are introduced in red. The 
straight line represents the average of data for as-received water in silica measured by 
several authors and quoted by Shelby [2.10]. The circles are results by Briickner [2.8] 


measured by water removal (also quoted in [2.10]). 


Since the studies resulting in the red line ([2.7-2.10]) showed reasonable agreement, we 
have used this curve for our computations. For a discussion of the trustworthiness of the 


different results for silica, see Shelby [2.10]. 


Glass Composition Po Ap Ev/ Cy, 
(mol-%) (wt-%) (g/cm?) PCy eq.(2.1.5a) 
B B2O03-glass - 1.83 6.3 —5.3 
K K:0-SiO>-glass 28.2-71.8 2.38 1.41 0.41 
(20-80) 
Ca N2>0-Ca0-SiO>-glass 16.6-9.4-74.1 2.49 0.46 0.54 
(16-10-74) 
Al Liz0-AlO3-SiO2-glass 10.7-9.1-80.2 2.33 0.37 0.63 
(20-5-75) 
Si SiO2 - 2.2 —0.84 1.84 


Table 2.2 Swelling and shrinking strains derived from density curves given by Scholze [2.27] using 
eq.(2.1.5a). 
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Whereas silica shows a density decrease with increasing water content, the densities of 
the other glasses increased with increasing Cy. The relative density changes Ap/po were 
transformed into volume strains £v according to eq.(2.1.3). The results are plotted in Fig. 
2.8b and compiled in Table 2.2. While the B2O3-glass shows strong volume shrinking 
due to water uptake, silica shows clear volume expansion. The effects on the other glass- 
es are rather small. Materials of Table 2.2 containing Ca and Al also swell. K shrinks 


slightly. 
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Fig. 2.8 a) Density change of glasses as a function of water content. Black curves given by Scholze [2.27], 
red line: average of data for as-received water in silica quoted by Shelby [2.10], circles: data from 
Brückner for water removal, also quoted in [2.10]. B: B203-glass, K: K2O-SiO2-glass (20-80 mol- 
%), Ca: NaxO-CaO-SiO2-glass (16-10-74 mol-%), Al: LizO-Al203-SiO2-glass (20-5-75 mol-%); b) 


volume strains computed from Fig. 2.8a via eq.(2.1.3). 


The swelling effect is evident also from quantum mechanical computations by Zhu et al. 
[2.28]. These authors made quantum mechanics computations on a SiO2-nanorod of 36 
SiO2 atoms to which they added a single water molecule and determined the effect of the 


hydroxyl] generation on stress-strain behaviour. 
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From their computations, a linear swelling stress of Osw.z=—1.44 GPa was found by the 
hydroxyl generation. From the axial modulus of 189 GPa it follows a /inear swelling 
strain of Eswz = 0.8%. The numerical results must of course depend on the special SiO2- 
arrangement (axially orientated rod) and the mechanical boundary conditions of the rod, 


nevertheless two points become evident: 
a) The material must swell by the silica-water reaction. 


b) The swelling stresses can become very large (in this special case —1.44 GPa). 
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3 Effect of stress-dependent diffusivity 


3.1 Computation of diffusion profiles for 
temperatures <500°C 


The diffusivity is a function of stress, commonly expressed by the hydrostatic stress 
component, On, given by eq.(2.1.23). The diffusivity for the case of stress-enhanced dif- 


fusion is given by the following equation [3.1] 


AV. 
D= D, exp| o, — (3.1.1) 
o il h `g 


where D» denotes the value of the diffusivity in the absence of a stress. Tis the absolute 
temperature in X; AV, is the activation volume for stress-enhanced diffusion and R is 


the universal gas constant. 


When an externally applied stress O.ppı and swelling stresses Osw are present simultane- 


ously, the total hydrostatic stress is 
oO; = Oh appl + Cn sw (3.1.2) 


The swelling stress depends linearly on the water concentration, on « C. The saturation 


value of Onw for C=Q is in the following considerations denoted as ono. 


Due to the swelling stresses, the diffusivity is a function of water concentration. Conse- 
quently, diffusion of water is now governed by the partial differential equation for the 


uniaxial diffusion 


ôC ô po% (3.1.3) 
ot oz Oz 


Equation (3.1.3) was numerically solved using the Mathematica procedure NDSolve 
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3 Effect of stress-dependent diffusivity 


[3.2]. Figure 3.1a shows water profiles as a function of time with the actual water con- 
centrations normalized on the value at the surface C(O). The solid curves represent the 
solution in the absence of any stress effects on the diffusivity. With increasing 
time, the widths of the diffusion profiles increase. The bold dash-dotted curve in Fig. 


3.1a is the limit case for saturation conditions Q0) = © [3.3]: 


C=C, ot (3.1.4) 

aD: 
In Fig. 3.1b the depth Zı,2 at which the water concentration decreased to the half of 
the surface value is plotted as a function of time. Whereas the diffusion zones in 
the absence of swelling stresses increase with time, they are clearly reduced with increas- 
ing time in the case of swelling included. This effect was already visible from Fig. 


3.1a. 


For reasons of simplicity, we introduce a normalized dimensionless time t and normal- 


ized depth coordinate 6, defined by 


h? 


= a 2 (3.1.5) 
T D, tg [Bat 
Equations (1.2.5) and (1.2.6) then read 
Gest) = ert £| = exp[6Vr + lert $ + J| (3.1.6) 
C, 2 2 
and C(0, r)/C, =1-exp[c]erfe[ vr]. (3.1.7) 


In the absence of swelling, analytical solutions for the limit-cases t=0 and t—>œ can be 
derived as was outlined in detail in [3.4, 3.5]. For the ratio C({€,t)/C(0,t) it results from 
eqs.(3.1.6) and (3.1.7) 


an tA oer 


C(0,7) 1- exp[terfe[ Vz ] 


(3.1.8) 
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3.1 Computation of diffusion profiles for temperatures <500° 


0.6 


0.2 


Fig. 3.1 a) Diffusion profiles as a function of time without and with consideration of swelling stress effect 


on diffusivity, b) diffusion depth C12 versus time, arrows indicate limit values for x. 


At very short times, we obtain by a series expansion with respect to qt that by setting t—>0 


reads 


CE,0) exp] E va cent $] (3.1.9) 
Cc(0,0) 4 2 2 


At very long times >, only the first term on the right-hand side of eq.(3.1.8) remains 
finite. Consequently, we obtain the well-known solution for constant surface concentra- 


tion: 


CG, o) = ert $| (3.1.10) 


These limit cases are represented in Fig. 3.1a by the dash-dotted and the blue bold curve. 
The depths at which these limit distributions decrease to C{C,t)/C(0,t)=1/2 are 


Ca =0.6995 for 70 
61/2 = 0.9538 for To 


(3.1.11) 
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3 Effect of stress-dependent diffusivity 


In presence of swelling, the solution for to is derived in [3.4, 3.5] and reads 


AV. 
GT pie exp new ert £| (3.1.12) 
GC Tp AV, RT 2 


The red curve in Fig. 3.1a shows this solution for onoA Vw/R7=-3. The depth C12 at which 
AE) AO) decreases to C(C)/C(0)=C(C)/Q=1/2, is 


1 
= ert!) co,-—__-___ (3.1.13) 
Gin = 2er |. natal 


where erf ! is the inverse error function and the abbreviation X stands for 


a FHV, (3.1.14) 
RT 


Numerical values for dr are compiled in Table 1 in [3.5]. The limit values from (3.1.13) 


are introduced in Fig. 3.1b by the arrows. 


3.2 Semi-analytical approach 


The numerical effort for the time-dependent solution of the diffusion equation, eq.(3.1.3), 
in presence of swelling and stress-dependent diffusivity is hardly manageable. An ap- 
proximate analytical solution can be obtained following the general procedure usual in 
perturbation theory. If we consider the diffusivity D(C) as the disturbance parameter, 
perturbation theory suggests to solve the problem for the case of D(0)=D» and to insert 
the disturbance parameter into this solution. The disturbed solution is then obtained by 
introducing D(C) instead of D» into the undisturbed solution. Since the analytical 
solution of the undisturbed problem, here denoted as the first-order solution C‘') is 


given by 


CY Z h h? 7 t 
= erfc ex z + —t Jerfci +h | (3.2.1) 
Co eon iD D, l 55. a 
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3.2 Semi-analytical approach 


the analytical second-order solution C® simply reads 


(2) 


= erfc u 
C, Ee (3.2.2) 


h? 


h Z t 
= ex E Z+ DIG”) de, Die +h | 


with the diffusivity, eq.(3.1.1), taken at the local- and time-dependent concentration 
C=O", 


A comparison of the perturbation solutions and the numerical results from eq.(3.1.3) is 
given in Fig. 3.2. Whereas for onoA V«/RT= —2 (a) the iteration solution C® is in com- 
plete agreement with the numerical solution, first small deviations are visible for 
OnoA Vw/ RT=—3, Fig. 3.2b. It is therefore suggested to use the perturbation solution if the 


stress term is not too strong. 


ir by 
CPC, 
OnoAVw 
0.87 RT 
f t=1 
0.6 
og C® and 
numerical solution 
0.2r 
a) 
L 1 1 
0.5 1 1.5 


Fig. 3.2 a), b) Semi-analytical solutions C’” compared with numerical results for x} /D=1. 
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3 Effect of stress-dependent diffusivity 


Higher iteration solutions are possible by introducing D(C(n)) on the right-hand side, 


resulting in Art), etc. 


cm 7 
= erfc 


Cy 2J D(C™)t (8.2.3) 


ex a Z+ h? t lerfc Z +h : 
Dem" DE”) 2 Dem): VDE") 


Figure 3.3a shows the effect of the maximum swelling stress ono on the shape of the con- 
centration profiles. Figure 3.3b represents the time-dependent surface concentration 
normalized on the saturation concentration. The higher the compressive saturation swell- 


ing stress is the earlier is saturation reached. 


C/C(0) C(0)/Co 
i 1 


0.8 0.8 


0.6 0.6 
0.4 


0.4 


0.2 0.2 


Fig. 3.3 a) Effect of the saturation stress ono, on the shape of the profiles for >», b) surface concentrations 


vs. time. 
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3.3 Comparison with experimental results from literature 


3.3 Comparison with experimental results from 
literature 


Profiles of time-dependent residual water content via stepwise surface removal were 
measured by Davis and Tomozawa [3.6] in water vapour of 355mm Hg pressure at 350°C. 
The results, Fig. 3.4, showed reduced diffusion depths with increasing soaking time. This 
behavior obviously agrees amazingly with the trend shown by the dashed curves in Fig. 


3.1a. It can be interpreted as further evidence for the occurrence of swelling stresses. 


0.5 
time (h) 


Residual Hydroxy! Content/Total Uptake 


2 
Z/Vt (umh) 


Fig. 3.4 Residual water content during stepwise surface removal Diffusion profiles by Davis and Tomozawa 


[3.6]. 


Effective diffusivity from diffusion profiles 

Approximate methods for the determination of the diffusivity from measurements of 
water profile and water uptake were proposed by Tomozawa and co-workers [3.6, 
3.7]. Basis of their evaluation is the use of an equation that is only valid for a constant 


surface concentration, although the concentrations clearly increased with time: 
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3 Effect of stress-dependent diffusivity 


C =C(0)erfe] ——— (3.3.1) 
ADi 


One possibility for the computation of an “effective diffusivity” is via the measured depth 


Zır at which the concentration decreased to 50% of the surface concentration. 


l Tn oAVw 
10712] 4 = OH uptake Des/Do 
Der Im @ OH profile 
eff | A H20 uptake 
s & v H20 profile 
(cm?/s)| 
[ v 
a) 

10° 1 1 1 ı ı 


100 200 300 400 500 
t h 


Fig. 3.5 a) Measurements of effective diffusivities by Oehler and Tomozawa [3.7] at 250°C and 39 bar 
water vapour pressure, b) effective diffusivity computed from the curves for the diffusion zone 


half-width 21/2 of Fig. 3.1b. 


3.4 Influence of externally applied stresses 


Beside the swelling stresses also externally applied stresses as present for instance in 
bending tests must affect the diffusion profiles. The diffusivity in the presence of a hydro- 


static swelling stress term on and an externally applied hydrostatic stress On.appı iS 
AV, 
D= D, ex (Osna + Ci a) e| (3.4.1) 


For bending tests, the stresses are sufficiently constant within the range of the diffusion 
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3.4 Influence of externally applied stresses 


layers. In such tests, the upper and the lower surface of a bending bar show tensile and 


compressive stresses. 


Ir 1r 
[ ===" =0 | 0AVw _ 
RT RT ~ 
0.8 F 0.8 F 


0.6 
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0.2 
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Gr an E E E E E Ye Ey 
0.5 1 1.5 2 
G 


Fig. 3.6 Diffusion profiles under externally applied stresses; effect of swelling for soaking time ~1. 
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3 Effect of stress-dependent diffusivity 


Concentration profiles computed with eqs.(3.1.1) and (3.1.3) are plotted in Fig. 3.6 and 
Fig. 3.7 for varied hydrostatic swelling terms onoA Vw/RT and varied time. It is obvious 
that in all cases surface concentrations are smaller under externally applied tensile stresses 
than in the absence of such stresses. Compressive stresses cause an increase of hydroxyl 
water at the surface. The enlarged surface concentrations for compression decrease with a 
steeper profile than under tension. In deeper regions the concentrations under tension 


exceed those for the stress-free and the compression case. 


Since the swelling stresses ono depend on the saturation concentration which is pro- 
portional to the vapour pressure and since the activation volume A Vw in eq.(3.1.1) is not 


known sufficiently, the parameter onoA Vs/RT was tentatively varied in Fig. 3.6. 


Comparison of Fig. 3.6a to Fig. 3.6c shows clearly that the reason of the difference in 
surface concentrations is not caused by the swelling stresses since the effect occurs even 
in the absence of swelling, (Fig. 3.6a). With increasing swelling the difference at the 
surface is only slightly increased. From Fig. 3.6a we therefore have to conclude that the 
effect of the externally applied stresses on surface concentrations is the increase of water 


with time, i.e. the parameter / in eqs.(1.2.4-1.2.6). 


In Fig. 3.7 the time is varied for constant swelling On oA Vw/R7=-2. The results indicate 
that for short times the concentration differences at the surfaces are larger than for longer 
times. For t — œ the differences must of course disappear completely because then the 


water concentration reaches saturation independent on the applied stress. 


The stress effect on the surface concentration can simply be summarized: 
Under tensile stresses the diffusivity is increased and diffusion is enhanced. This must 


result in 


e an increased zone thickness compared to tests without externally applied stresses. 


e _ Due to the higher diffusivity under tensile stresses, the water diffusion into the bulk 
material is also increased. Since the water entrance from the environment into the 
surface is limited by the finite value of A, the concentration of water at the surface 


must decrease. 
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3.4 Influence of externally applied stresses 


Under compressive stresses the diffusivity is reduced. This results in a reduced penetration 


depth and an increased surface concentration. 


1r 17 
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Fig. 3.7 Diffusion profiles under externally applied stresses; effect of normalized soaking time, t, for 


OnoA W/RT=2. 
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3 Effect of stress-dependent diffusivity 
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Fig. 3.8 Concentration profiles by Nogami and Tomozawa [3.8] as a function of the hydrostatic stress, a) 
bending bars for silica at 192°C (water vapour pressure: 12.3 bars), b) externally pressurized spec- 


imens at @=350°C (water vapour pressure: 25 bars), C: in absorption coefficients (1/mm). 


Figure 3.8a shows experimental results on bending bars by Nogami and Tomozawa [3.8] 
at 192°C carried out after 300h treatment in water atmosphere at 12 bars water vapour 
pressure. Measurements were carried out on the tensile and compression side under load 


and in an additional test also on unstressed bars. 


In Fig. 3.8b similar results at 350°C are shown for 5h soaking (water atmosphere at 25 
bars). In this case, Nogami and Tomozawa [3.8] applied hydrostatic pressures directly. 
The results of Fig. 3.8 completely confirm the computed curves in Fig. 3.6 and Fig. 3.7. 
The effective activation volume was experimentally determined in [3.8] as A Ver =170 


cm’/mol for 192°C and 72 cm?/mol at 350°C. 


For the interpretation of the results in Fig. 3.5, the soaking temperature of 250°C is con- 
sidered. Fig. 3.9a shows hydroxyl concentration S-data in mass units under water satura- 
tion pressure as could be computed from measurements by Zouine et al.[3.3] using the 
equilibrium equation from Wiederhorn et al. [3.9]. The straight line introduced in this 


plot reads 
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3.4 Influence of externally applied stresses 


S.. = 0.000265 exp(0.0143 0) (3.4.2) 


sat — 


with the temperature 9 in °C. An extrapolation to 0=250°C results in the saturation con- 
centration of S=0.94 wt%. Under the assumption that even at moderate temperatures only 
hydroxyl generation contributes to the volume change, the corresponding hydrostatic 


swelling stress term at a free surface is 


~_2KE S=-18.7GPaxS (3.4.3) 


Ow h Ar 
i 9(1-v) 
i.e. Ono=-177 MPa. Consequently, it holds for 250°C: 


On AV, 


=—0.0407 AV, (3.4.4) 
1r P 
L / LF Ol h,0AVw 
S i _— 
0.5 F 
(wt%) | 
0.2 F 
AVy=72 cm?/mol 
0.1 f 
0.05 F 
lı “O 1 N 1 N N 
0 50 100 150 200 250 
8 (°C) T 
Fig. 3.9 a) Concentration of hydroxyl water at silica surfaces under saturation pressure by Zouine et al. 


[3.3], b) results of Fig. 3.5b together with the curve related to the effective volume A V=72 


cm?/mol (red curve) as was derived by Nogami and Tomozawa [3.8]. 
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3 Effect of stress-dependent diffusivity 


In Fig. 3.9b the results of Fig. 3.5b are shown once more together with the parameter 
chosen for the results obtained by Nogami and Tomozawa [3.8], namely, A Vw =72 
cm?/mol (the value of AV, =170 cm?/mol seems to be extremely large). For the parameter 
(3.4.4) it results: onoA Vw/RT=—2.93. The related D/D) dependency is shown in Fig. 3.9b 


by the red curve. 


Finally, Fig. 3.10 illustrates the ratio Det/D) from Fig. 3.9 as a function of the swelling 
parameter onoA Vs/RT by the circles. The dash-dotted straight line represents eq.(3.1.1). 
The red arrows indicate the activation volumes suggested by Nogami and Tomozawa 
[3.8]. For strongly negative swelling parameters the integration procedure for the numer- 
ical solution of the differential equation (3.1.3) becomes rather unstable. A rough estima- 
tion for the very large value of 170 cm*/mol needs an extrapolation of computed results. 
Such an extrapolation is tentatively shown by the dashed extension of the computed data 


curve. 
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Fig. 3.10 Diffusivity ratio of Deti/D) from Fig. 3.9b for normalized time t=3 versus the swelling parameter 
OnoA Vw/RT (symbols) compared with the diffusivity ratio D/D» according to eq.(3.1.1) (dash-dotted 


line); the dashed line is tentatively introduced for extrapolations. 
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3.4 Influence of externally applied stresses 


A modified representation of the results in Fig. 3.6 and 3.7 was given in [3.10] as plotted 
in Fig. 3.11 for the lower temperature region <500°C, where simply SocC: In Fig. 3.11 
the hydroxyl concentration at the surface, S, is normalized on the asymptotically reached 
value for disappearing external load, S®, as is for instance fulfilled in the neutral axis of 
a bent bar. The hydroxyl content is under compression higher than under tension, Fig. 
3.11a. If reaction (1.1.1) shows stress-enhancement as will be discussed in Section 5, the 
S-concentration in tension is by a factor Sœ exp[constx o/ RT] increased and in compres- 


sion decreased with the result given in Fig. 3.11b. 


For long times when the water concentration approaches saturation, the S-content in 
tension exceeds that in compression due to the stress-enhanced equilibrium constant. The 


influence of stresses changes at the circle in Fig. 3.11b. 
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Fig. 3.11 a) Hydroxyl concentrations at a silica surface as a function of soaking time ¢based on the stress- 
affected diffusivity, exclusively; b) for stress-enhanced equilibrium of reaction (1.1.1). S® is the 
saturation concentration in the absence of external stresses. Diffusivity parameter 


|oz\/3x Vw/RT=0.5, swelling parameter Osw.n.o Vw/RT=-—3. Curves plotted for temperatures <500°C. 


Hydrogen concentrations close to silica surfaces were reported by Agarwal et al. [3.11]. 


These authors used the nuclear reaction analysis (NRA) on water profiles in specimens 
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3 Effect of stress-dependent diffusivity 


that were hot-water soaked at 250°C under 355 Torr vapour pressure and a bending stress 
of 60 MPa. These results are represented in Fig. 3.12a. The water content under pressure 
(compressive side of a bent specimen) is always higher than under tension. This is in 


complete agreement with the data by Nogami and Tomozawa [3.8]. 


Since no surface concentrations are given in [3.11], we plotted in Fig. 3.12b the data at a 
depth of 0.12 um below the surface versus soaking time. The slopes roughly fulfil the 
surface predictions from Fig. 3.11 for short times, namely C/G œ Vt No indication for 


saturation of the surface concentration is detectable. 
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Fig. 3.12 a) Water concentration below the surface at 250°C and 355Torr vapour pressure under bending 


stresses of 60MPa, b) data at a depth of 0.12 um below the surface plotted vs. soaking time. 


Figure 3.13 shows the concentration of hydrogen atoms, H, in flexure specimens of silica 
glass, exposed to water vapour of 355 Torr at 650°C. The maximum bending stress in the 
specimens was +45 MPa. Under stationary conditions, an increase of the concentration at 


the tension side and a decrease at the compressive surface were found. 
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3.4 Influence of externally applied stresses 


These data confirm clearly the effect of higher water concentration under compression 
for short time and higher water concentration in tension for longer times. They are also in 
agreement with curves in Fig. 3.11a. The only difference is caused by taking the square- 


root for the ordinate since SoC” for temperatures >500°C. 
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Fig. 3.13 Hydrogen concentration below the surface at 650°C and 355Torr vapour pressure under bending 


stresses of 45MPa. 


Conclusions 
There are several important findings from the computation of diffusion profiles for water 


vapour soaked silica. The computations show that: 


. in the absence of swelling stresses the diffusion profile must tend for long time 


against the well-known description by the complementary error function; 


° in presence of swelling stresses in the diffusion layer, the profiles must become 


steeper with increasing time; 


° the effective diffusivity, Der, defined by the distance 2,2 =; Dyt at which the 


water concentration decreased to 50% of the surface value, must increase slightly if 
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no swelling stresses would occur. In presence of swelling, it should decrease strong- 


ly. 


The agreement of computations and the experimental results by Davis and Tomozawa 


[3.6] and by Oehler and Tomozawa [3.7] clearly show that swelling of silica is present. 


From the computations, also the results by Nogami and Tomozawa [3.8] and Agarwal et 


al. [3.11] are easily understandable. 


3.5 Diffusion in thin specimens 


3.5.1 Computation of swelling stresses 


The swelling stresses given by eqs.(2.1.22, 2.1.23) are valid for the case of thin swelling 


layers negligible compared with the specimen widths 2 W. Now let us consider thin spec- 


imens. In this case, the negative hydrostatic swelling stresses in the surface region where 


the diffusion zones develop must decrease after longer time. This is due to the fact, that 


all stresses must disappear when a constant water concentration is reached after long 


times. Then no restrictions on swelling are possible and the specimen can freely expand. 


The diffusivity trivially must tend to the value Dp valid for on=0. 


The condition of mechanical equilibrium requires for the total stress distribution over the 


thickness 2 W that 
2W 
Í odz =0 
0 


This condition is fulfilled by the effective swelling stress © eff 


z l E = 
O off (z) = 0,,(Z) Ow = 3 l-v (E, é,) 


as can simply be validated by inserting (3.5.2) into (3.5.1). 
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(3.5.1) 


(3.5.2) 


3.5 Diffusion in thin specimens 
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Fig. 3.14 a) Deviations of concentrations from their average over the width. Red curves: before maximum is 
reached, blue curves: after the maximum (the ordinates are proportional to the swelling stresses), b) 
concentrations for B=10 dependency of surface values with time, c) difference between surface and 


average concentration. 


Figure 3.14 shows the influence of the time on the stress distributions (Figs. 3.14a, 
3.14b) and the stresses at the surface (Fig. 3.14c). The stress values are scaled with the 


swelling stress Oswo that would occur at the surface of a semi-infinite body for t+. In 
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3 Effect of stress-dependent diffusivity 


case of the low-temperature reaction with equilibrium constant Ai it holds for equi-biaxial 


swelling stress 


o = IE ckC, (3.5.3) 
sw,0 
3 1-v 


As a consequence of eqs.(3.5.2) and (3.5.3), it can be written with C defined by 
eq.(1.4.6) 


Oy (Zt) C(z,t)—C(t) (3.5.4) 
o Co 


sw,0 


For stress-enhanced swelling the equilibrium constant kı depends on the swelling stress. 


Then it must hold 


Op (Zt)  S(z,t)- S(t) (3.5.5) 
oO ki oCo 


sw,0 
where io is the equilibrium constant in the absence of any stress. 


Figure 3.14a represents the deviations of the local water concentrations from the average 
value for a Biot number of B=1, ie. the right hand side of eq.(3.5.4). The red curves 
show the results for times before the maximum surface value and the blue curves illus- 
trate the distributions after the maximum. Figure 3.14c gives the surface values as a func- 


tion of normalized time t for different B. 


3.5.2 Expectations from first-order analysis 


Swelling of silica has several consequences on material behaviour. Two effects on diffu- 
sion behaviour may be addressed in this section. Since the analytical solution by the 
eqs.(1.4.1a) to (1.4.3a) does not include stress effects, these equations allow only to per- 


form first-order estimations. 


In [3.10] we could show that the surface water concentration increases at the surface 


under compressive stresses and decreases under tensile stresses. These findings are a 
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3.5 Diffusion in thin specimens 


simple consequence of the /imited mass-transfer coefficient [3.10] or s/ow-reaction 


boundary condition [3.12] at the silica surface. 


0.8 


0.7 
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Short time limit __------ 
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oe = 
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Long-term limit 


Fig. 3.15 Effect of diffusion boundary conditions at the silica surface and stress-affected diffusivity, 


eq.(3.4.1), on the normalized surface concentration of water, G(t)/C. Dashed curve: Limit case for 
the half-space under swelling at the surface; Solid curve: Limit case for diffusion in the absence of 


swelling stresses; Coloured curves: Possible limit-case interpolations. 


Two limit cases may be considered: 


(1) For very short times, ¢ the thickness of the water diffusion zone, b=(D xd!” is small 


(2) 


compared to the specimen thickness, b<< W. Consequently, diffusion behaves like 
in a semi-infinite body. This limit case was outlined in [3.10]. The surface concen- 


tration G for onoA Vv /R T=-3 is shown in Fig. 3.15 as the dashed curve. 


For very long times, the specimen is completely soaked by water, resulting in the 
concentration Qz, JQ. Any swelling stress must disappear in the specimen as a 


consequence of mechanical equilibrium, eqs.(3.5.2, 3.5.4). The solution for the half- 
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3 Effect of stress-dependent diffusivity 


space, computed with a disappearing stress-enhancement of onoA Vw /R T=0, is rep- 


resented in Fig. 3.15 by the solid limit curve. 


The real diffusion behaviour must start on the dashed curve for short times and tend to 
the solid curve for long times. This is schematically indicated by the coloured thin 
curves. Depending on the material parameters this change can occur very slowly (blue 


curve) or more abruptly (red curve). 


From Fig. 3.15 it can be expected that for long times for which the layer thickness band 
the specimen width Wbecome comparable, a maximum in the surface concentration C() 


occurs. 


3.5.3 Numerical results for stress-enhanced diffusivity 


Equation (3.1.3) was numerically solved using the Mathematica procedure NDSolve 


[3.2]. Two boundary conditions had to be satisfied simultaneously, namely 


Na (3.5.6) 
dz D`’ 

age acy. aw (3.5.7) 
dz D 


In the following considerations, it is assumed that the coefficient A might be independent 


of stress. The swelling stress related to a certain concentration Cis 


le 6, le) (3.5.8) 
where the quantity ono is the hydrostatic swelling stress value reached for C=Q. Figure 
3.16 shows the surface values of the effective swelling stresses defined by eq.(3.5.2) as a 
function of time t, Biot number B and hydrostatic swelling stresses onoA Vw/RT=-3 
[3.13]. The upper limit curve with B>% represents the case of the half-space, since 
Bo W. In this case, any diffusion layer is small compared to the specimen width and an 


unloading effect with increasing layer thickness is impossible. 
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3.5 Diffusion in thin specimens 


Fig. 3.16 Effective swelling stresses at the surface. 
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3 Effect of stress-dependent diffusivity 


On, oAVw 
RT 


Fig. 3.17 Water concentration at the surface (effect of swelling stress, Biot number B, and time T). 


Figure 3.17 represents the surface concentration as a function of time t and Biot number 
B. An interesting effect is visible for stronger swelling as illustrated in Figs. 3.17a-3.17c 


for several Biot numbers. 


These figures show a relative concentration maximum, followed by a temporary decrease 
of C(O). For large times the concentration © is asymptotically approached. 
It should be mentioned that for stronger swelling, |onoA Vso/R7]>4.5, the Mathematica 


Routine NDSolve did no longer converge. 
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4 Surface concentrations and swelling 
stresses 


4.1 Hydroxyl concentration under water saturation 
pressure 


4.1.1 Results from literature 


Hydroxyl concentrations at silica surfaces below 500°C are available from investigations 
by Tomozawa and co-workers [4.1,4.2] for water vapour pressure of 355 Torr and for 


water saturation pressure by Oehler and Tomozawa [4.3] and Zouine et al. [4.4]. 
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Fig. 4.1 Results by Zouine et al. [4.4], a) saturation pressure with symbols indicating the test conditions, b) 


solubility of water at silica surfaces under saturation pressure. 
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The results from Zouine et al. [4.4] are shown in Fig. 4.1. The saturation pressures p, 
under which the tests were carried out, are represented in Fig. 4.la. For these measure- 
ments the nuclear reaction analysis (NRA) was applied that subsumes molecular and 
hydroxyl water. Whereas in [4.1-4.3] the S-concentration was measured via the IR- 
method, the data by Zouine et al. [4.4] represent the total water content Cy at the surface. 
The results are shown in Fig. 4.1b. From these measurements, also the water species $ 


and Ccan be obtained. 


4.1.2 Computation of molecular and hydroxyl water 


concentrations 


In molar units, the total water concentration is given by 


C,, =C +45 =C(1+4k) (4.1.1) 


where the quantity kis the equilibrium constant describing the ratio of A=.S/C. The exper- 
imental results on equilibrium ratios from literature were expressed in [4.5] for the tem- 


perature range of 90°C< 7<350°C by the empirical relation 
pes => (4.1.2) 
C RT 


(A=32.3 and Q=10.75 kJ/mol). Equations (4.1.1) and (4.1.2) result in 


Gia Gs, (4.1.3) 

1+4k 

C 

s=- 7 (4.1.4) 

1 1 

pan + Pair 

2 k 

and in mass units considering the different mole masses of water and hydroxyl 

Se 17 C, (4.1.4a) 


ae D 
—+— 
2 k 
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4.1 Hydroxyl concentration under water saturation pressure 


Dividing Cy, C, and Sby the density of silica yields the water species in parts per mass of 
glass, mostly in ppm or wt-%. For the temperature range of 20°C<0<200°C, the straight- 


line dependence of saturation pressure, Fig. 4.1a, can be approximated as 


B, 
Psa = B, exf- =| (4.1.5) 


(Bi=5.79x107 kPa, B>=41.28 kJ/mol). Under the common assumption that the concentra- 
tion of molecular water species at the surface is proportional to the pressure of the water 
vapour, Cp, the combination with the two 1/7-depenencies for k and p makes clear that 
the total water concentration at the surface, eq.(4.1.1) not necessarily shows a straight 
line in an Arrhenius plot. This behaviour is visible from the dashed curve in Fig. 4.1b. 
The concentrations of molecular and hydroxyl water were computed from the measured 
total water concentration by using eq.(4.1.3) and (4.1.4a) with the result plotted in Fig. 


4.2 versus 1/7: The dashed curves are guidelines for the eyes. 
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Fig. 4.2 Surface concentration of the molecular (blue) and hydroxyl (red) water species obtained via 
eqs.(4.1.3), (4.1.4a). 
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4.1.3 Semi-logarithmic description 


Figure 4.3a shows the total water concentration in a semi-logarithmic plot with linear 
temperature scaling. The measurements can approximately be described by the simple 


expression 


C,, = 0.000780 exp(0.00868 9) for 23°C<0<200°C (4.1.6) 


with the temperature 0 in °C. 
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Fig. 4.3 a) Total water from Fig. 4.1b and hydroxy] data of Fig. 4.2 plotted with linear abscissa scaling, b) 
hydroxyl water concentration, curve: eq.(4.1.7) obtained by combining eqs.(4.1.6) and (4.1.2). 


Figure 4.3b illustrates the S-data in a similar plot with a slightly curved averaging curve 


given by 


_ ız 0.000780 exp(0.00868 0) (4.1.7) 


S 
8  14texp(Q/RT) 


with the parameters A and Qin eq.(4.1.2). For simple computations in the reduced tem- 
perature range of 70°C<0<200°C, tolerating a few percentage deviations from the curve 


in Fig. 4.3b, a straight-line representation is suggested that reads 
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4.1 Hydroxyl concentration under water saturation pressure 


S = 0.000265 exp(0.0143 0) for 70°C<0<200°C (4.1.7a) 


This approximation is introduced in Fig. 4.4a as the red line. For extrapolations, 
eq.(4.1.7a) should not be used. Already at 8=250°C, the hydroxyl concentration is over- 
estimated by about 30%. 
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Fig. 4.4 a) Hydroxyl data of Fig. 4.3b together with the simple straight-line approximation by eq.(4.1.7a), b) 


Hydrostatic swelling stresses under saturation water pressure. 


The swelling stresses at the glass surfaces are given by eqs.(2.1.11) and (2.1.22) as 


O swz T Oswy TTE 2 S =-28 GPa xS (4.1.8) 
! y 3(l=v) 


Introducing the hydroxyl concentrations from Fig. 4.2 results in the swelling stresses 
plotted in Fig. 4.4b. For the computation of the swelling stresses, the necessary hydroxy] 
concentration can be determined via eq.(4.1.4a) with the equilibrium constant according 
to eq.(4.1.2). Introducing the hydroxyl concentrations from eq.(4.1.7) results in the swell- 


ing stresses 


_ _ 13.75MPa exp(0.00868 6) (4.1.9) 
a 44 Lexp(Q/RT) 
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4.2 Hydroxyl concentration under constant pressure 


In order to eliminate the effect of increasing surface concentrations with increasing satu- 
ration pressure, the concentrations were normalized on the pressure of 355 Torr (see also 
[4.4]). In contrast to Fig. 4.2, the surface concentrations decrease with increasing temper- 
ature. This makes clear that the increasing curves of Fig. 4.1b and Fig. 4.2 reflect strong- 


ly the increasing saturation pressure. 
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Fig. 4.5 Concentrations of molecular and hydroxy] water from Fig. 4.2, normalized on a constant water 


vapour pressure of p=355 Torr. 


The data in Fig. 4.5, describing concentrations of molecular water and hydroxyl in the 
absence of externally applied stresses, can be represented by straight lines for /7> 
0.0022/K: 


Con © exp] 18.9442> beni met (4.2.1) 


and 
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4.2 Hydroxyl concentration under constant pressure 


S(O)asstox = exp - 15:36 ZZ (42.2) 


RT 
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5 Stress-enhanced equilibrium 


5.1 Thermodynamic description 


At temperatures T> 500°C, the equilibrium constant of the reaction 


=Si-O-Si= +H20 © =SiOH+HOSi= (5.1.1) 
is 
2 
ER (5.1.2) 
C 


According to Le Chatelier [5.1], the equation governing the equilibrium constant is 


ôlnk __ AV 
Op RT 


, AV =2V,-Ve- (5.1.3) 


where p is pressure, AV is the change in the partial molar volume of the reaction, V, 
and V, the partial molar volumes of hydroxyl and molecular water, respectively. R is the 


universal gas constant and Tis the temperature in °K. At high temperatures, some of the 
water will react with the silica network according to eq.(5.1.1) to form hydroxyl groups, 
S= [=SiOH], in the silica network. The concentration of the unreacted water is C= [H20]. 
We believe that the molecular water in the glass is located in holes of molecular size of 
the silica glass network [5.2], where the water does not contribute to the volume of the 
glass. As a consequence, molecular water cannot contribute to the stress effect. Thus, the 
concentration of molecular water in the glass depends only on the vapour pressure in the 


external environment, but not on the stress in the glass. 


This same conclusion was reached earlier by Agarwal et a/. [5.3] using slightly different 
reasoning. They noted that noble gases such as helium diffuse through silica glass with 
little sensitivity to the applied stress [5.4]. Water molecules are believed to diffuse 
through silica glass in much the same way as the noble gases. Hence, the diffusion of 


water molecules is also expected to show little dependence on applied stress. Following 
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this reasoning, all of the sensitivity of the hydrogen concentration to the applied stress in 


the experiments by Argawal et al. [5.3] resides in the hydroxyl concentration. 


This is equivalent to V, =0 in eq.(5.1.3). Consequently, molecular water cannot contrib- 
ute to the stress effect as had been shown by deformations on water-soaked disks [5.5]. 
Thus, the concentration of molecular water in the glass depends only on the vapour pres- 
sure in the external environment, but not on the stress in the glass. By replacing the hy- 


drostatic pressure p by the hydrostatic stress on in a solid 


o, =4(0, 0) (5.1.4) 


we obtain with the hydroxyl concentration So for on=0 


o,V; 
S=S his ; mn (5.1.5) 
vel RT ) On P 


At temperatures T< 450°C, hydroxyls are immobile. Only those hydroxyls can react in 
the reverse reaction step, which are directly neighboring. Consequently, the reaction is of 
first order. In this case the equilibrium constant for the reaction is 


k= (5.1.6) 


ala 


The result for the hydroxyl concentration is therefore of the same form as given by 


eq.(5.1.5) possibly with a different reaction volume V in the different temperature re- 
gions of <500°C and >500°C. 


5.2 Experimental evidence for stress-enhancement 
5.2.1 Water concentration at crack surfaces of subcritically 
grown cracks 


An interesting experiment was performed by Tomozawa et al. [5.6]. Fracture surfaces 
were formed by the passage of a crack through Double Cleavage Drilled Compression 


specimens (DCDC) [5.7, 5.8], while the specimen was exposed to water. The experi- 
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ments demonstrated that the surface concentration of water was significantly higher near 
the fracture surface than would be the case for silica glass merely exposed to water for a 


time equal to that required for the passage of the crack through the specimen. 


The study [5.6] used nuclear reaction analysis (NRA) to determine the concentration pro- 
file of water near the fracture surface of the glass. This technique provides a direct chemi- 


cal analysis of the concentration of H as a function of distance from the fracture surface. 
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Fig. 5.1 Hydrogen concentration in the surface region of a growing crack via NRA-measurements by To- 
mozawa et al. [5.6] (circles: from Fig. 4, squares: from Fig. 3 in [5.6]); a) linear ordinate: curve 
given by eq.(5.2.1), red bar indicates 90%-Confidence Interval of the surface concentration, b) log- 


arithmic ordinate. 


In Fig. 5.1 the results from [5.6] are plotted versus depth (circles). We included also the 
high values from Fig. 3 in [5.6] (indicated in Fig. 5.1 by the squares). It should be noted 
that the measured data are affected by the finite breadth of the instrument curve. A de- 
tailed analysis that includes the effect of the rather large breath of the instrumental is giv- 
en by Wiederhorn et al. [5.9]. Here we are interested only in the surface values for which 


the effect of the instrument curve is of minor importance. 
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In order to allow an approximate extrapolation of the measured data to the surface, we 


fitted the data cu in Fig. 5.1a arbitrarily by a complementary error function 


Cy =Cyp ere Z| +0.005%10 /cm’ (5.2.1) 


Background 


taking into account the background level of 0.005 10”/cm?. The coefficient cno is the 
surface concentration of cu. By fitting eq.(5.2.1) to the measurements, the best set of 


parameters was found as 
cx0=1.30 [1.055, 1.55] 10”/cm?, b=4.95 [3.28, 6.61] nm (90% CD. 


with the 90%-Confidence Intervals in brackets. The red bar in Fig. 5.la indicates the 
90%-CI of the surface concentration. The depth bat which the concentration decreased to 
=50% of the surface value in Fig. 5.1a is about 5 nm. 

Using Avogadro’s number, N=6.02x107? molecules/mol gives for the water concentra- 


tion at the surface in molar units: 
C (0) = T =0.0108 [0.0095, 0.0121] (mol/cm’) (5.2.2) 


Defining the mass concentration as the mass of water per mass of glass, the surface value 
is with p=2.2 g/cm? for the density of silica glass and a molecular mass of water of 


my20=18 g/mol 


C, (0) =- 


m 
= Su 0 _ 8 817.7,9.8] (wt%) (5.2.3) 
m p 


2N 


glass 


For high water concentrations it is of advantage to use the total mass m=my+Melass in the 


dominator 


C,,(0) =" — =8.1[7.1, 9.0] (wt%) (5.2.4) 


glass w 


From the result by Zouine et al. [5.10] for unstressed silica surfaces in water, cao=1.4 


10” H-atoms, the water concentration in the absence of stress, Cv, reads 
C,,. = 0.000116 (mol/cm’) = 0.095 (wt %) (5.2.5) 


The value of œ is shown by the blue line in Fig. 5.1b and after including the background 


by the dash dotted line. From these results, an increase of the water concentration by a 
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factor of C.(0)/Cv093 has to be concluded. The right ordinate in Fig. 5.1b represents the 


total water concentration. 


5.2.2 Computation of the equilibrium constant 


In molar units, the total water concentration is given by 


C,, =C +45 =C(1+4k) (5.2.6) 


where the quantity X is the equilibrium constant given by eq.(5.1.6). For the case of satu- 
ration, the surface concentration of molecular water, C(0) is generally assumed to depend 


only on the water vapour pressure in the environment, C(0)=C=const. This leads to 


C, =C, +48 =C,(1+4k) (5.2.7) 


In the absence of stresses, the equilibrium constant may be denoted as ko. In this case, the 


total water concentration can be written 


Cpo =Cy + 589 = Col + Fk) (5.2.8) 
again with subscript “0” for disappearing stress. 


The ratio of total water concentrations under stress and in the absence of stresses is 


C, _ l+ak (5.2.9) 
C 1++k, 


w,0 
Since the left-hand side is known from experimental data, the stress-enhanced equilibri- 


um constant can be computed as 


k= (24h) 29 (5.2.10) 


w,0 


Experimental results on equilibrium ratios from literature were expressed in [5.5] for the 
temperature range of 90°C<@<350°C by the empirical relation (1.1.3). 
For room temperature, we obtain by extrapolation to 23°C (7=296°K) 


k,(23°C) = 0.41 
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5 Stress-enhanced equilibrium 


The experimental data for Cw/ Gw, in eq.(5.2.10) yield 


k=222 > <= 540 (5.2.11) 


0 


an enormous increase of the equilibrium ratio by the stress field near the crack tip. The 
result clearly illustrates the strong effect of tensile stresses. The equilibrium of the reac- 
tion (5.1.1) is strongly shifted to the right side. It is evident that the strong increase of the 
equilibrium constant must be caused by a positive reaction volume AV >0 in (5.1.5) as 


was derived in [5.11]. From eqs.(5.1.5) and (5.1.6) it follows simply 


S_k (5.2.12) 


The hydroxyl concentrations from eqs.(5.2.3) and (5.2.5) are 


S(0) = ZC, (0) (4 +4) =16.5[14.4,18.4] (wt%) (5.2.13) 


and So = ŽC, o (E+ L) = 0.0368 (wt%) (5.2.14) 


The factor 17/18 is the ratio of the mole masses of hydroxyl (17 g/mol) and molecular 
water (18 g/mol). In [5.12] we derived from the spherical pore model a maximum possi- 
ble hydroxyl concentration of Smax=15.7 [13.3, 19.2] (wt%). The 90%-confidence inter- 
vals of S(O) and Smax totally overlap confirming the expectaion that the hydroxyl concen- 


tration at the crack tip is the maximal possible one. 


5.2.3 Bending tests by Agarwal et al. 


In this section, we use the experiments of Agarwal et al. [5.3] to derive information on 
stress-enhanced swelling of silica glass. Figure 5.2 from Agarwal et al. [5.3] shows the 
concentration of hydrogen atoms, H, in thin flexure specimens of silica glass, exposed to 
water vapour at 350 Torr and 650°C. The maximum bending stress in the specimens was 
+45 MPa. We evaluated the data in Fig. 5.2 for the times, £=10 h and 20 h. As indicated 
by the work of Davis and Tomozawa [5.13], water and hence hydrogen at 650°C is pre- 


sent primarily in the form of =SiOH groups. 
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5.2 Experimental evidence for stress-enhancement 


The concentration of H in atoms per cm? of glass in Fig. 5.2 was obtained by the tech- 
nique of Nuclear Reaction Analysis (NRA) [5.3, 5.14]. The OH concentration, S, in gm 
of OH per gm of silica glass, was calculated from the concentration of hydrogen atoms, 


H, using the following relation: 


3 
S(gmOH / gm silica glass) = H cm’ x17g,, mol ; (5.2.15) 


6.02 x10” /mol: 2.2g „„/cm 
a 3 0.8 
] S 
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Fig. 5.2 Hydroxyl water concentration below the surface of bent silica glass specimens: water vapour pressure 
355 Torr; temperature 650°C; exposure time 4 h, 10 h and 20 h, taken from Agarwal et al. [5.3], shaded areas: 
+1 SD. 


The 10 data points in the tensile region at +45 MPa showed a mean value of $** = 0.638 
wt% with a standard deviation of SD = 0.022 wt%. The standard error (SE) of the mean 
is SE = 0.0071 wt%, or S** = (0.638+0.0071) wt%. 


The 11 data points in the compressive side showed a mean value of S“ = 0.539 wt%, 
SD = 0.0173 wt%, and SE = 0.00523 wt%, hence, S* = (0.539+0.00523) wt%. 
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5 Stress-enhanced equilibrium 


Using the error propagation law gives for the ratio of the two concentrations 


S*°  0.638+0.0071 _ 0.638+1.1% 


= =1.184+1.5% (5.2.16) 
s®  0.539+0.00523 0.539 +0.97% 


The hydroxyl concentration at the neutral position of the cross-section is determined as 
the geometric mean of S'* and S*, i.e., (S5 x S®)"?. Since Agarwal et al. [5.3] mention 
that the hydrogen concentration of the unstressed samples fell between the concentrations 
of samples under tension and compression, we will identify this mean concentration with 


the concentration in the absence of an applied stress: 
S(0) = 0.586 wt% . (5.2.17) 


The concentration of H was measured by Nuclear Reaction Analysis. Normally, H is 
divided between H2O and =SiOH, but at 650°C almost all of the H is present in the form 
of =SiOH. Hence, the measured H concentration can be considered equal to the concen- 


tration of =SiOH. 


In [5.3] the concentration of =SiOH was measured additionally by IR. These measure- 
ments showed 5"/,S‘ = 1.19 that is in good agreement with .S*/S = 1.184 from the NRA 
measurements. Within experimental scatter, the results were the same. The absolute value 
of So however depended substantially on which technique was used for the measurement. 
The NRA measurements were a factor of three times larger than the equivalent IR meas- 


urements. The reason for this discrepancy is not yet clear. 


The experiments in [5.3] were carried out under uniaxial stresses (45 MPa on the tension 
and -45 MPa on the compression side). For an evaluation of these results, a relation may 


be used similar to eq.(5.1.5) that describes S(oappı), where Oappı = +45 MPa: 


S-S, exp Zeita | (5.2.18) 


A Vet is an effective activation volume that will of course deviate from V, in the 


eq.(5.1.5). The ratio of the concentrations in tension and compression is consequently 


5” exp Fa Ver (5.2.19) 
= RT 
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5.2 Experimental evidence for stress-enhancement 


From the bending experiment with 5"/S’=1.184, eq.(5.2.16), we obtain for the case of 
stress-enhanced swelling and 7>450°C 


Vg = 14.4cm?/mol +10% (5.2.20) 
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6 Swelling anisotropy 


6.1 Motivation for anisotropic swelling 


Strongly different stress states can result in the same hydrostatic stress term on. An ex- 
ample for this fact is illustrated in Fig. 6.1 where the hydrostatic stress disappears com- 
pletely, on=0. This is trivial when all stress components disappear, Fig. 6.1a. In the case 
of Fig. 6.1b where the transversal stresses are half of the longitudinal stress with inverse 
sign, it is again on=0 according to Gn=(0x+0y+0,)/3. 

The same holds for torsion loading, oy=-o:z, Fig. 6.1b’. It is hardly imaginable that high 
tensile stresses in bond direction should have no effect on bond fracture. It has to be 
expected that the reaction in tensile direction should produce swelling even for the hy- 


drostatic stress state on=0. 


Serer rrr rrr rrerer 


Fig. 6.1 Three stress states acting on a Si-O-bond exhibiting the same disappearing hydrostatic stress on=0: 
a) isotropic swelling is expected, for b) and b’) swelling in z-direction should dominate, (case b’ 
describes torsion loading), c) water molecule reacting with a stressed bond oriented in z-direction, 


producing two hydroxyls, 2... 
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6 Swelling anisotropy 


Figure 6.2 shows the load-displacement curves for bond breaking by reaction (1.1.1). 
Perpendicular to the loading direction, Fig. 6.2a, the necessary energy for splitting the 


bond is Uthat must be delivered by thermal fluctuations. 


The portion M (subscript x means the bonds oriented in x-direction in which no stress 


acts) of opened bonds results from the Boltzmann equation as 


N, exp- (6.1.1) 


The same equation holds for the y-direction 


U 
N exp| — — (6.1.2) 
ye | a 


The potential barrier Ubelongs to the case when no stresses are present or when the bond 


is perpendicular to the stress, Fig. 6.2a. 


Fig. 6.2 Energy conditions for bond breaking. 
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6.1 Motivation for anisotropic swelling 


Under action of a force Fin bond-direction, Fig. 6.2b, the potential barrier is reduced to 


U-FSı. Consequently, the occurrence of opened bonds in z-direction is increased 


where Vis an activation volume. In the following considerations, it will be shown that 


the strongest anisotropy effect has to be expected for the torsion loading case. 


For the computations on anisotropic swelling let us follow the derivations made in [6.1]. 
The stress effect on the equilibrium constant for reaction (1.1.1) can be derived from 
Eq.(6) of the article by de Boer [6.2]. Following the procedure in the review article by 


Hamann [6.3], we obtain the following equation for the equilibrium constant: 
RT Ink = Homo — 2 Hosion + 2AW + 2AU , (6.1.4) 


where pomo and po,sion are the chemical potentials for molecular water and hydroxyls, 
respectively in the unstressed state; A Wrepresents the work per mol of SiOH done by the 
mechanical stresses; and AU is the change of the elastic strain energy per mol due to a 
reduction of material stiffness as a consequence of the water reaction. 

The work term A Win Eq.(6.1.4) is for loading in 3 directions ox, Oy, Oz 


AW/V, = 0,69), + o,E0 +0,E% (6.1.5) 


z“ sw,z y“ sw,y x" sw,x 


with the Cartesian components of swelling strains &w,i, the glass volume, Vo, convenient- 
ly chosen as the molar volume of dry glass, and the superscripts ‘” standing for 1 mole of 
S. Since the total swelling volume is in general the sum of the components, it holds for 
the individual strain components &w,” and the volume swelling strain &‘ for 1 mole of 


Sin the volume Vo 


SWZ 


VD p0 a) a) 
bp SER ee re (6.1.6) 
The change of strain energy is for uniaxial loading by O= Oappi 


AU /V,=-18}, AE (6.1.7) 


27 appl 
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6 Swelling anisotropy 


where AFis the change in the Young’s modulus as a consequence of the material damage 


by hydroxyl] generation. 


The reduction in stiffness results in a change of the quantity of elastically stored energy, 
which in turn, contributes to the driving force of the reaction [6.4]. Since AZX0, the con- 
tribution AU has the same sign as AW. In the evaluation in [6.1], the energy term, AU, 
was neglected assuming that AZAS is small due to the low hydroxyl concentrations in 
our calculations. Hence, we use the undamaged Young’s modulus in all equations having 


in mind that the results of our computations are lower limits. 


Equations (6.1.4) and (6.1.5) now yield 


a) 0) 0) 
k u k exp 2V (OE wz FOE wy + OLE sy) 4 (6.1.8) 
0 0 RT 
and, consequently 
a) (a) () 
S =S,expl V (Osu t FyE sy + Ex) | (6.1.9) 
0 0 RT 


where ko includes all of the terms that are independent of stress (uo,sion, Ho,n2o and In[C)). 


The molar volume for glass is Vo=27.27 cm?/mol. For one mole of hydroxyl with the 
mass Mou=17g/mole in the volume Vo we obtain the hydroxyl concentration in weight 


units (mass S/mass SiO2) 


go __Mou__ 17 g/mole _ 4 59, (6.1.10) 
M zus 60 g/mole 
with the related volume strain according to eq.(2.1.11) 
e® =S” =0.274 (6.1.11) 
Then it holds for Vey 
Vie = 27.27 x 0.274 = 7.5 cm’/mole (6.1.12) 
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6.2 Pure tensile loading 


6.2 Pure tensile loading 


Under purely uniaxial tensile stresses 62 = Oappı >0, Ox= Oy =0, it results 


a) 
On 
S=S,exp V, le (6.2.1) 
0 0 RT | 
or rewritten 
> fe) AY) Toy ED 6.2.2 
q v SW,Z a] SW,Z 3 
S = S, exp} V, ai FU = S, exp RT FU x7.5cm’/mol (6.2.2) 


by introducing eq.(6.1.12) in (6.2.1). 


The experimentally found dependency S=f(Gappi) in Fig. 5.2 must agree with the theoreti- 


cal relation eq.(6.2.2), i.e. it must hold for uniaxial loading Oappi=oz: 


a) 


AV, 
P appt eff _ Z apo! esta x 7.5cm*/mol 
RT RT e? 
(6.2.3) 
eu, _ 14.4cm*/mol 
ma „ 1,9 
& 7.5 cm”/mol 
Consequently, we obtain: 
Esw/&=1.92. (6.2.4) 
and for the other components [6.1]: 
E Ei, 
Ey Et ED En EI = = 0.46 (6.2.5) 
E, E, 
Finally, it holds for the uniaxial load 
5 =S,exp) Zoe Var (6.2.6) 
RT 
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6 Swelling anisotropy 


6.3 Swelling stresses 


Since in water diffusion layers the swelling strains €sw,i are mostly present in thin surface 
layers, free expansion is not possible due to the bulk material that remains free of water. 
Due to this restriction, the swelling strains must cause swelling stresses Osw,. and Osw,y in 
the surface layer. These swelling stresses are responsible for the fact that even under 


purely uniaxial loading a biaxial stress state prevails near the surface. 


Then, the total stress in the surface layer is 


0, = O appl,z + O swz $ O, = O appl, y + Ow y (6.3.1) 
For biaxial stresses we get from (6.1.9) 
a) 0) 
O,€5,, tOy€ 
S = S, exp} V, (Enz ARE AR (6.3.2) 
RT 


In our special application an applied stress in y-direction is not present, Gappi,y = 0, so that 


eq.(6.3.2) reads after renaming Oappi,z = Oappı and use of eq.(6.2.5) 


-046 Var 
S= So exp (om H O swz H 1.92 m RT | (6.3.3) 


Hooke’s law gives the total strains in a thin surface layer as the sum of the swelling and 


the elastic strains 


€ total,z zz E swz + gO = vo...) (6.3.4a) 
1 
E total,y = E my H E (Oswy VO oy 7) (6.3.4b) 


With &orasz = Etay =0, equations (6.3.4a) and (6.3.4b) are simultaneous equations with 


two unknowns: swz and Oswy. The solution of this system of equations is 


epee ete we (6.3.5a) 
v 
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6.3 Swelling stresses 


oO = eu Em + VE wz) (6.3.5b) 
v 


Introducing eqs.(6.3.5a, 6.3.5b) and (6.2.5) results in the exponent of eq.(6.3.3) 


E é 
Opa to ayy =e | + (5 Av) - 201-20) (6.3.6) 
Z g 7 Al-v)leu. E 


v 


so that eq.(6.3.2) can be written as an implicit equation with respect to $ 
AV. 
S=S, ox (a + BE (6.3.7) 


The abbreviation A reads by using the material data Z=72GPa, v=0.17, and «=0.97 


v 


psn E E E E EET (6.3.8) 
4v) Ew € 


The explicit form of eq.(6.3.7) is: 


RT 


AV, AV, 
S PLog eff B Sy EXP| © ppi Sei is (6.3.7a) 
AP RT RT 


where PLog stands for the product-logarithm function (see e.g. [6.5]): 


The axial swelling stress results from eqs.(6.3.5a) and (6.1.6) as 


Ow 57K -= aı-" |+% |xS=-7,S, n,=133 GPa (6.3.9) 
l-v ZN 2 


with &=€sw,/Ev as an anisotropy coefficient. 
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7 Strength of silica fibers 


7.1 Tensile strength of fibers 


Soaking of silica under tensile loading improves strength. To the knowledge of the autors, 
this increase was qualitatively found first by Proctor et al. (Section 3.3 in [7.1]). In a 
recent study by Lezzi et al. [7.2], uniaxial tensile stresses, 1 GPa to 3 GPa, were applied 
to silica fibers, which were then exposed to water vapour, 6 Torr vapour pressure. After 1 
min of exposure, followed by cooling to room temperature, fiber strengths were measured 
in normal lab air at 20°C. Figure 7.la shows the strengths as a function of soaking 
temperature for different tensile stresses Gap. Perpendicular bars show + 1 standard 


deviation. 


of | b) 


(GPa), 


200°C 
1 L 1 L 1 ! ae oer eee ee E A 
100 200 300 400 500 600 1 2 
T co) Oappı (GPa) 
Fig. 7.1 a) Strength measurements on silica fibers thermally treated for 60 s in nitrogen atmosphere with a 


partial water pressure of 6 Torr under different tensile stresses; b) strength at 200°C vs. applied 


stress during thermal treatment (results from Lezzi et al. [7.2]). 


103 


7 Strength of silica fibers 


In the evaluations of strength we concentrate on 200°C and 500°C since at these 


temperatures the best data base is reported. 


In Fig. 7.1b the 200°C-data are plotted in some more detail. For 500°C, Fig. 7.2a shows 


the increase of critical stresses, i.e. of strength, Ao,, 


AG, =O, (Fag) — F-(0) (7.1.1) 


as a function of the applied stress where o.(0) is the strength of unstressed fibers. The 


circles indicate the median values and the bars represent + one standard deviation. 
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Fig. 7.2 a) Strength results on silica fibers at 500°C, Red data: corrected strength; b) Effect of 2 test 
specimens cracked during soaking at 3GPa. These specimens were ignored in strength evaluation 
by Lezzi et al. [7.2]. 


In our opinion the measured strength data at 3 GPa are slightly too high [7.3]. The 
authors mention that at the highest load some fibers (about 10%) fractured already during 
soaking. Under constant test conditions, these fibers, which were inevitably excluded 
from the later strength tests, must have had the largest defects. This corresponds to the 


well-known problem of a proof-test preceding the strength measurements. In principle, 
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7.2 Strength predictions including swelling 


one has to evaluate a strength distribution truncated at the lower strength values [7.4]. 
The way of evaluation has of course an effect on the strength of the surviving specimens 
and may be illustrated in Fig. 7.2b for a Weibull modulus of m=8. The black curve 
represents the Weibull distribution for the case that no specimen would have failed 
during the soaking procedure. For the case that two specimens might have failed during 
soaking, the failure probability has to be computed on the basis of 17 specimens and not 
on the basis of 15 specimens. The true characteristic strength oo is 1.5% smaller than the 
value reported in [7.2] and the median value by 3% smaller. The change in Ao, is then 
slightly reduced by 4.2%xAo, = -95 MPa. Figure 7.2a shows the results after re- 
evaluation as the red data point together with the original data. The slight shift in Gap» is 


for an improved distinction. 


7.2 Strength predictions including swelling 


7.2.1 Effect of global swelling stresses 


The swelling stresses in Section 3 were obtained under condition of the absence of exter- 
nal loading. So we could always assume isotropic swelling, ie. strains were assumed to 
be identical in all directions. As this might not be the case for some recent high tensile 
stress experiments carried out on silica glass [7.2], a reconsideration of the equilibrium 
conditions is necessary according to Section 6.3 for a valid determination of the effect of 
applied stress on the equilibrium constant for reaction (1.1.1). In this section information 
on stress-enhanced swelling of silica glass is used to discuss the experiments of Lezzi et 


al. [7.2] in terms of water penetration and the consequent swelling of silica glass. 


In the high-stress soaking tests on silica fibers by Lezzi et al. [7.2], the externally applied 
loads were uniaxial. Therefore, we concentrate here also on the uniaxial description of 
the effect of stress on =SiOH concentration in the glass by applying o, exclusively. As 
outlined in Section 6.3, the swelling stress in axial direction is according to eq.(6.3.9) and 
[7.3] 

O =—133xS (GPa) (7.2.1) 


sw,z — 


The strength increase with respect to the strength of unstressed fibers is 
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7 Strength of silica fibers 


Ao = —Ow.z (Capp) F Ow. (0) (7.2.2) 


c 


where Osw,z(0) is the swelling stress in the absence of an externally applied stress, O.ppi=0. 


Figure 7.3 shows the hydroxyl concentration at a water vapour pressure of 355 Torr as a 
function of temperature, reported by Davis and Tomozawa [7.5]. The circles in Fig. 7.3 
indicate data from IR-analyses. The red triangle is from the IR-data in reference [7.6], 
whereas the blue triangle is from the NRA analysis in reference [7.8]. The blue square is 
an extrapolated result, based on the slope of the Davis and Tomozawa data [7.5] and the 


position of the blue triangle. 


Assuming that the ratio Syr/Sirx3 is constant from 650°C to 500°C, we can use the 
general temperature trend shown in Fig. 7.3 to scale the value of S from 650°C to 500°C. 
From the result S(0) = 0.586 wt% at 650°C, we obtain from Fig. 7.3 for 500°C: 
S(0)=0.586x1.75=1.026 wt%. This concentration results in the concentration of 
S0,500°C,355Tor= 1.28 wt%. Finally, the concentrations, So, were transformed from 355 Torr 


to 6 Torr water vapour pressure (as was used by Lezzi et al. [7.2]) via 


oe Oi (7.2.3) 


=S 
6 Torr 355 Torr 355 


=> So,s00°c, 6Torr = 0.167 wt.% 


Table 7.1 compiles the S-concentrations derived. 


Temperature S0) at 355Torr (0) from Fig. 7.3 (355 Torr) Sy at 6 Torr 
(wt%) (wt%) 
NRA IR NRA IR NRA IR 
650°C 0.586 0.19 s H 
500°C 1.026 0.333 0.167 | 0.047 


Table 7.1 Hydroxyl concentrations in wt-% used for the strength predictions; S(0) at 650°C represents the 
hydroxy] content in the neutral axis of the bending tests by Agarwal et al. [7.6]. These results were 


used in Fig. 7.3. 
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7.2 Strength predictions including swelling 
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Fig. 7.3 The hydroxyl concentration Sin silica, reported by Davis and Tomozawa [7.5] at a water vapour 
pressure of p=355 Torr (circles); the red triangle indicates the result of the IR-evaluation in [7.6]. 


The blue triangle represents the concentration obtained by an NRA analysis in [7.6]. 


7.2.2 Application on strength 


Strength at 500°C: In a recent study by Lezzi et al. [7.2], uniaxial tensile stresses, 1 GPa 
to 3 GPa, were applied to silica fibers; the fibers were then exposed to water vapour, 6 
Torr vapour pressure at 500°C. After 1 min of exposure, followed by cooling to room 
temperature, fiber strengths were measured in normal lab air at 20°C. No measurements 
of hydroxyl water by the NRA-method are available for that experiment for comparison 
with our computations. However, Agarwal et al. [7.6] measured the hydroxy] distribution 
at 650°C by both NRA and IR-analysis. Therefore, a comparison of the two measure- 
ments is possible at that temperature. From the IR measurements in Fig. 7b of their paper, 


we get the following surface values from the bent specimen: 
S= 20.3 Abs/cm, S* = 17 Abs/cm, S(0) = 18.6 Abs/cm 


(1 Abs unit/cm = 0.1 wt% OH). The ratio of the concentrations at the two surfaces, S*/.S” 
= 1.19, is essentially the same as found for the NRA-results presented above, S'/S = 
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1.184. By contrast, the value of S(O) for the unstressed specimen differs from the NRA 
method by a factor of about three. For the IR measurement S(0) = 0.19 wt%, whereas for 
the NRA measurement S(0) = 0.586 wt%. 
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Fig. 7.4 a) Strength increase at 500°C vs. applied stress during thermal treatment for anisotropic swelling; 
Circles: results from Lezzi et al., bars: + 1 SD, curves: Effect of global swelling stresses -Osw 
computed with the product of a Vos=14.4 cm*/mole, a=1.92 and So from NR-analysis (solid curve) 
and from IR-data (dashed curve), b) sensitivity of strength prediction on the anisotropy coefficient a. 


For the red data point, see Section 7.1. 


Figure 7.4 shows the strength increase measured by Lezzi et al. [7.2] as a function of the 
applied stress. The circles indicate the median values and the perpendicular bars +1 
standard deviation. As already mentioned, the measured strength data at 3 GPa given in 
[7.2] are slightly too high. A correction of this data point is given in Section 7.1, resulting 


in the red circle. 


The swelling stresses are plotted in Fig. 7.4a by the solid curve for the Smeasurement 
via NRA and the dashed curve for the IR measurement; both values were calculated from 
Agarwal et al. [7.7] and are indicated in Fig. 7.3 by the triangles. In Fig 7.4b the strength 
predictions via eqs.(7.1.1) and (7.1.2) are shown for different values of the anisotropy 


coefficient a. It can be stated that the influence of æ is not strong. 
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7.2 Strength predictions including swelling 


Strength at 200°C: Figure 7.1b shows the strengths as a function of the applied stress 
during soaking at 200°C. For the prediction of strength a modification of the derivation in 
[7.3] is necessary. In the low-temperature region of 0<500°C, the equilibrium constant of 


reaction (1.1.1) is 


ee S (7.2.4) 
C 
The procedure in the article by Hamann [7.8] results in the equilibrium constant: 
RT Ink, = fo mo — Hosion tAW + AU » (7.2.5) 
The equilibrium constant under stress is 
k, =k, on V, Orfa t = + Ena | (7.2.6) 


Combining eqs.(7.2.6) and (7.2.4) yields the stress-enhanced hydroxyl concentration 


a) a) a) 
(O, Ew. + Oyewy T O xE wy x (7.2.7) 
RT 


S=S, os 


This is the same relation as was derived in [7.3] for the high-temperature region of 
02500°C. This makes it possible to take over the main results from [7.3].The hydroxyl 
concentration data in mass units under water saturation pressure are approximated in 


Section 4 as 


S,. = 0.000265 exp(0.0143 0) (7.2.8) 


sat — 


with the temperature O in °C. For 200°C we obtain Sjx=0.5 wt%. This concentration is of 
course not the surface value in a stress-free state since this concentration causes swelling 
stresses and is therefore the concentration under a compressive swelling stress. In the 
absence of externally applied uniaxial loading, the swelling is isotropic. This implies 
(Section 6.3) 


o=1/3 and Bi3=-56.1 GPa 


The hydroxyl concentration So results from 
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ly g® 

So = S sat exp - bi S sat Voß, = 0.6 wt% (7.2.9) 
RT 

Since saturation pressure at 200°C is 1554 kPa, the concentration for p=6 Torr is 


So, 6 Torr = 3.1 ppm 


Introducing this value in eq.(7.2.7) gives the curve ‚S(O.ppı) Shown in Fig. 7.5a by the right 


ordinate. 
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Fig. 7.5 a) Effect of soaking under stress for 6 Torr/200°C: left ordinate: strength increase, right ordinate: 


related hydroxyl concentration, b) prediction of strength increase (curve). 


The axial swelling stress was obtained from eq.(7.1.1) and the strength increase from 
eq.(7.2.2). For a stress of 1 GPa there is not any strength increase visible and even for 2 
GPa the increase of tensile strength is only #250 MPa. For a stress of 3 GPa we can ex- 
pect a strength increase of about 0.9 GPa and a hydroxyl concentration of S*0.7 wt-%. 
Figure 7.5b finally compares the measurements by Lezzi et al. [7.2] from Fig. 7.la 
together with the predictions. The measurements are again reported by the solid symbols, 
the predictions by the open ones. Figure 7.6 shows the strength predictions for the 
different applied tensile stresses together with the measured data as a function of soaking 


temperature. 
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7.2 Strength predictions including swelling 


In the preceding computations only the effect of global swelling in the diffusion layer 
was taken into account. For the alternate strength description by crack-like defects, it is of 
importance to compare the size of the surface cracks, a, with the thickness of the water- 
affected zone, b, given by eq.(2.2.2). From Zouine et al. [7.7] we obtain D,=8x 10°"? cm?/s. 
Consequently: by 70 nm. 


predicted for Oappi=3 GPa 
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Fig. 7.6 Predictions (open symbols) compared with the strength measurements of Fig. 7.1a. 
The crack depths can be concluded from strength results. Assuming the surface cracks as 
edge cracks, the crack depth is given by the fracture condition A= Kj, and 
K.=0,Fıra (7.2.10) 


with the geometric function /=1.122. For kj.=0.8 MPaxm [7.9] and o=4.6 GPa one 
obtains from (7.2.10): 


a=7.6 nm 


Since a<<b, cracks are fully embedded in the water-affected zone. 
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7.3 Effect of local swelling at surface cracks 


During the soaking treatment the crack-tip region is surrounded by a swelling zone. The 


shielding stress intensity factor caused by such a zone is given by [7.10] 


Keen (73.1) 
v 


where @ is the height of the swelling zone. This shielding stress intensity factor causes an 


increase of inert strength Ao, 


Ks. en EE o, (7.3.2) 


The zone of an arrested crack is heart-shaped resulting in the coefficient y=0. After high- 
temperature soaking, these swelling zones are maintained at room temperature where the 


strength tests are performed. 


In the strength tests, the crack must extend and escape of the initial swelling zone. Con- 
sequently, the coeffient yin eq.(7.3.1) must become dependent on the crack extension Aa 
in the strength test. By application of the McMeeking and Evans [7.10] procedure as well 
as by FE-analysis, y was computed in [7.11]. The red curve in Fig. 7.7 represents 
analytical results for a heart-shaped swelling zone we got by application of the 
McMeeking and Evans [7.10] procedure (see also Appendix Al). A maximum of y is 
reached at Aa/@=8/V27=1.54 with wx0.3. In the case of a grown crack with extended 
swelling zone (black curve) it results yx0.33. In order to evaluate eq.(7.3.2), the zone 


height œ has to be known. 


For room temperature Lechenault et al. [7.12] used a neutron reflection technique to 
measure the penetration of heavy water (deuterium oxide) into the silica glass. The work 
of Lechenault er al. [7.12] predicts the water profiles for crack velocities of = 10° m/s and 
=4x10° m/s. An average zone height of &=6.5 nm may be concluded for room tem- 
perature. This value is a lower limit for high-temperature applications since at 500°C the 


diffusivity is increased and increased zone heights have to be expected. 


The negative intrinsic stress intensity factor causes an increased externally applied stress 


intensity factor for crack extension. Consequently, the strength must increase. 
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7.3 Effect of local swelling at surface cracks 


Due to the water content in the crack-tip region, the Young’s module is reduced. Conse- 


quently, eq.(7.3.1) must overestimate the shielding effect and (7.3.2) the strength increase. 
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Fig. 7.7 a) Swelling zone for an arrested crack, b) grown crack for weak swelling, c) shielding stress 
intensity factor for crack extension through the initially swelling zone, red curve: crack passing the 
initial zone developing for an arrested crack, obtained by application of the procedure by 


McMeeking and Evans [7.10]; black: grown crack after additional extension Aa. 
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8 Prediction of residual deformation 


8.1 Swelling stress 


In Section 7 and [8.1] we predicted the strength increase of silica fibers by swelling and 
compared it with measurements by Lezzi et al. [8.2]. Good agreement was found by 
using the results from Agarwal et al. [8.3] obtained with the Nuclear Resonance Analysis 
(NRA). In this Section, we compare the swelling effect on residual bending displace- 


ments with deformation measurements by Tomozawa et al. [8.4]. 


In Section 6.3, we derived an equation for the hydroxyl concentration $ as a function of 


the applied stress, Sapp, eq.(6.3.7a) 


S RT PLog ABS, oe. Z) (8.1.1) 
AB RT el RT 


with %=14.4 cm?/mol and the parameter ß=-135 GPa. The swelling stresses in axial 


direction are given by eq.(6.3.9) as 
Oy, =—-133xS (GPa) (8.1.2) 


In this section, the bending moments and curvatures caused by these swelling stresses 


shall be computed and compared with measurements from [8.4]. 


Knowledge of the parameter A in the compression zone of a bent fiber is hardly critical 
for computations. In compression it principally holds 0< S<6% i.e. allowing only a small 
region of variation. On the other hand, the stress-affected zone in compression is clearly 
smaller than in tension due to the stress-affected diffusivity. Consequently, the value of 
the parameter A in the pressure range may therefore be approximated to be equal to the 


tensile case. 
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8 Prediction of residual deformation 


8.2 Effect of swelling zones on bending deformation 


Volume expansion under tensile loading must have consequences on the deformation 
behaviour, similar to the strength behaviour discussed in [8.1]. Under bending load, the 
effect of stress on the SiO./water reaction (1.1.1) is stronger in the tensile region than in 
the compression region. The enhanced swelling under tension and suppressed swelling 
under compression must result in an internal bending moment acting against the applied 


one. 


In order to give a transparent analysis, let us consider a thin bent fiber of radius x, Fig. 


8.la, loaded by a constant displacement between the fiber ends, Fig. 8.1b, and constant 


force F, Fig. 8.1c. This fiber is water vapour soaked under the externally applied outer 


fiber bending stress op. 


As long as b<<n is fulfilled, the bending stress in the diffusion zones is sufficiently 


constant. The applied stress at the surface as a function of the angle ọ is given by 


O appi = Op COS O (8.2.1) 


appl 


The thickness b of the swelling layer is governed by diffusivity and time via 


b= Dt (8.2.2) 


The diffusivity to be used in eq.(8.2.2) is a function of stresses, commonly expressed by 
the hydrostatic stress component op. When D(0)=D(Capp=0) denotes its value in the 
absence of an externally applied stress, the diffusivity for the case of stress-enhanced 
diffusion reads [8.5] 


AV, 
D=D() exp Ha a (8.2.3) 
The bending moment by swelling stresses can be written for << 
2a 
M,, =r [0,(P)b(p)cospdp (8.2.4) 
0 
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8.2 Effect of swelling zones on bending deformation 


When the swelling zone becomes comparable with the fiber radius, e.g. b>m/5, the 


moment has to be computed via 
2m To 


M,, = f [0,9 (r eosp)r dr do (8.2.5) 


g=0 r=0 


After unloading the internal swelling moment M,, causes a remaining curvature of 


un. Te, (8.2.6) 


load constant 
deformation constant 


b) c) 


Fig. 8.1 a) Fiber cross-section with swelling zone under bending, b) bending under constant displacements 


atthe ends, c) constant load. 


Equation (8.2.4) gives the bending moment by swelling for the case of very thin swelling 
zones with b/m—0. For thicker zones, comparable with the fiber radius 0, the solution 
becomes more complicated since then also the equilibrium conditions have to be fulfilled. 
Because of the Bernoulli hypothesis that plane cross sections must remain plane, the fiber 


curvature caused by swelling results in a linear strain distribution over the cross section: 
€,=C,X+C,, X=rcos@ (8.2.7) 


(for x see Fig. 8.1a). 
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8 Prediction of residual deformation 


The integrals of swelling stresses over the cross section of the fiber allow consideration 
of different types of tests. The requirements of disappearing normal force of swelling 
stresses and disappearing total moment result in two conditions from which the unknown 
coefficients cı and œ can be determined by integrating over the whole cross-section area 
A 


[Go,,+qx+0,)dA'=0 (8.2.8) 
(A) 

[Gon +0x+0,)xdA'=0 (8.2.9) 
A 


For constant load, the equations (8.2.8) and (8.2.9) have to be solved simultaneously. For 
this purpose we used the Mathematica [8.6] subroutine FindRoot. Boundary conditions 
for constant displacement are realized by setting cı=0. Consequently, eq.(8.1.1) now 


reads 


s RT PLog ABS, os (Gg + Ec, + Ec,rcosp)A (8.2.10) 
AB RT RT 


8.3 Predictions on the basis of diffusivity affected by 
applied stresses 


8.3.1 Results at 550°C 


Figure 8.2a shows the hydroxyl concentration at a water vapour pressure of 355 Torr as a 
function of temperature as reported by Davis and Tomozawa [8.7]. The data from IR- 
analyses are indicated by the circles. The ordinate shows the hydroxyl concentration 
normalized on the value at 650°C. At 550°C, we can conclude that the ratio $550°c/S¢so°c iS 
about 1.4. The values of 5S) expected from the NR-Analysis [8.3,8.1] are compiled in 
Table 8.1 for temperatures >550°C. The last column of Table 8.1 gives the diffusivities 
used for predictions according to eq.(8.2.3). The effect of the different loading conditions 
is shown in Fig. 8.3a. Since in 2-point-bending tests the bending moment changes along 
the fiber length axis and due to the fact that the water reaction only took place in a 


central region (indicated by the arrows in Fig. 8.1c), no well-defined mechanical bounda- 
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ry conditions 


curvature. 


are given. As limit cases, we consider constant moment and constant 


NRA NRA DO) 
Temperature | S/Sssosc | S(O) (wt%) | So (wt%) | (cm?/s) 
650°C 1 0.586 0.652 
550°C 1.4 0.821 0.971 1x10" 
625°C 1.07 0.63 0.71 3x10" 
700°C 0.87 0.51 0.56 7x10"! 


Table 8.1 Hydroxyl concentrations from NRA-measurements used for the deformation predictions according 


to [8.3,8.1,8.7]. 
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Fig. 8.2 a) Hydroxyl surface concentration S'in silica, normalized on the concentration at reference 


temperature 650°C, Ssso°c for a water vapour pressure of p=355 Torr (circles) (Davis and Tomozawa 


[8.7]); 


b) Diffusivity from [8.7] (circles); Squares: D(0) used for predictions with eq.(8.2.3). 


Computations were performed for a bending stress of 493 MPa at 355 Torr vapour 


pressure and a diffusivity of D(0)=1x10°' m?/s as indicated by the white square in Fig. 
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8.2b. Under constant load (cı#0, c>#0) the dash-dotted line resulted and for constant 
displacement (cı=0, #0), the solid curve was obtained (see Fig. 8.3). The circles 
represent the experimental data from [8.4]. The effect of the applied bending stress is 
visible from Fig. 8.3b. The results agree sufficiently with the experimental data from [8.4] 
for all applied bending stresses. The perpendicular lines indicate an upper limit for the 
time range in which the swelling layer thickness is smaller than 4=10um. It is self- 
evident that the residual curvatures cannot increase unlimited with time. For long times 
the diffused water has homogeneously soaked the full cross section, i.e. So=const. The 


saturation moment results then from eq.(8.2.5). 
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Fig. 8.3 Residual curvature of water vapour soaked silica fibers (125um diameter, 550°C), computed for 
A Vv=18 cm’/mol. Circles: results for constant displacement tests by Tomozawa et al. [8.4], a) 
Effect of different boundary conditions: Solid curve: constant displacement (Fig. 8.1b), dash-dotted 
curve: constant load (Fig. 8.1c); b) predictions for different loads and constant displacement 


conditions. 


8.3.2 Results at 625 and 700°C 


Residual deformations are also available at increased temperatures [8.4]. In addition to 
the results at 550°C (circles), Fig. 8.4 represents the data for 625°C (squares) and 700°C 
(triangles), all measured at 493 MPa. The diffusivities for computations are about D(0)=3 
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10°5m?/s at 625°C (blue square in Fig. 8.2b) and D(0)=7x 10°’ m?/s at 700°C (red square 
in Fig. 8.2b). The hydroxyl concentration at 625°C is about 7% higher than at 650°C and 
about 15% smaller at 700°C. For D(O) and Sp see Table 8.1. 


The solid lines in Fig. 8.4a represent the additional computations. Whereas for 550 a 
tenable agreement between measurement and computations is visible, at higher tempera- 
tures of 625 and 700°C the measurements are clearly underestimated by the compu- 
tations. At 625°C the predictions are 80% and at 700°C about 55% of the measured 
values. Under the assumption of the diffusivity D(0) as a free parameter, least-squares fits 
enable to determine the best temperature-dependent values. The resulting diffusivities are 
entered in Fig. 8.4b by the squares. The related straight-line in Fig. 8.4b is steeper than 
that for the results by Davis and Tomozawa [8.7]. 


Whereas for the data of [8.7] the activation is 81 kJ/mol, the value for the squares in Fig. 
8.4.b is obtained to be 122 kJ/mol. However, there is also the possibility that the glasses 
from references [8.4] and [8.7] may behave differently. 
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Fig. 8.4 a) Residual curvatures at different temperatures (NRA-data used). Symbols: Measurements by 
Tomozawa et al. [8.4]. Solid lines: predictions using the diffusivities D(0) from Fig. 8.2b, (Table 
8.1), dashed line: Least-squares fit to the measurements. In the case of 550°C the solid and dashed 


lines are identical; b) diffusivities obtained for Capp=493 MPa (squares) via data fit in a). 
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8.4 Predictions including the effect of swelling 
stresses on diffusivity 


The predictions made in Section 8.3 were performed with eq.(8.2.3) that implicitly 
assumes that the stresses acting in a volume element would be externally applied stresses 


exclusively. In presence of swelling stresses, this equation must read 


AV, 
D= D, Ex (Arm + | (8.4.1) 


The hydrostatic swelling stress term for anisotropic swelling, reads according to Section 
6.3 
KE 


A aa (l+a)S=-41GPaxS (8.4.2) 
is 6(1-v) 


—.l 
On sw = rile +o 


and for isotropic swelling with a=1/3 to be applied in the computation of Dy from D(0) 


for the unloaded case Oappi=0, 


Gru oF ge 18-7 Ghats (8.4.3) 
9(1-v) 


Equations (8.4.1) and (8.4.3) imply for the unknown diffusivity Də in the absence of any 


stress: 


D, = D(0)exp 18.7 (GPa) x S, £% (8.4.4) 

RT 
We varied the volume A Vy and computed the residual curvatures for the three stresses at 
550°C and compared the results with the measurements by using the Mathematica 
Subroutine NonlinearFit [8.6]. The best parameter set was found by matching the compu- 


ted curves to the measured data by Tomozawa et al. [8.4] as A V,=27 cm?/mol. 


The related set of curves is shown in Fig. 8.5a. Application to the higher temperatures 
again underestimates the measurements as is visible from Fig. 8.5b. It therefore remains 
to be stated that a moderate variation in AV, has no significant influence on the 


predictions. 
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Fig. 8.5 a) Curve fitting according to (8.4.1) for 9=550°C, resulting in A V,=27 cm*/mol, b) predictions for 
temperatures of 625 °C and 700°C. 


8.5 Recovery of curvature after unloading 


The hydroxyl concentration at the end of a bending test (Otena=493 MPa, t=36 h) is 
shown in Fig. 8.6a as the red curve. The S-distribution is non-symmetric with respect to 
the centreline. On the tensile side, the concentration S is enhanced due to the tensile 
stresses and on the compression side reduced. The hydroxyl] profile in the absence of 
externally applied stresses is represented by the black curve. This curve is of course 


symmetric. 


If the bent fiber is unloaded and cooled down to room temperature, the S' -distribution 
remains unchanged, and the S-profile caused by bending is “frozen in”. No time-depen- 


dent “recovery” of the residual curvature is possible. 


When the fiber is unloaded at the soaking temperature, the equilibrium constant must 
decrease because it is no longer stress-enhanced. The consequence is that the equilibrium 
of reaction (1.1.1) tends more to the left side generating molecular water. This is indi- 


cated by the downward arrow in the former tensile region. The time necessary for the 
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8 Prediction of residual deformation 


reaction can be rather short since diffusion processes are not necessary. In the compres- 
sion region, the equilibrium constant must increase and, consequently, the ‚S-concen- 
tration (indicated by an arrow upwards). This reaction is expected to be slow since the 


additional water must diffuse from the surface. 


The changes in S-concentration cause changes in volumetric strain as can be directly 
concluded from eq.(2.1.11). Changes in swelling strains give rise for a change in swelling 
stresses. In the former tensile region, the swelling stresses decrease and in the former 


compression region, swelling stresses increase. 


The mechanical response is a time—dependent reduction of the curvature (Fig. 8.6b). This 
type of “relaxation” is of course only a consequence of variation in S' with time as has to 


be concluded from the sequence of proportionalities 


S=f()>0,,0>M,,(0>1/R) (8.5.1) 
550°C, 493 MPa 1/R 
S 
L = (1/cm) 
nea t=36h 
0.4) 
0.3} 
Obend=493 MPa 
0.2} 
01! 
Obend=493 MPa l Obend=0O 
=! 1 1 L 1 L L 1 
20 40 60 80 100 120 2 4 6 8 10 12 
R-x (um) Vt m’) 


Fig. 8.6 a) Hydroxyl concentration along the diameter of a bent fiber (red curve) compared with the S'profile 
in the absence of bending, b) expected mechanical response of curvature for changes in the 


equilibrium constant and hydroxyl concentration. 
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9 Comparison between prediction and 
experiment 


9.1 Comparison of experimental results with 
predictions from swelling 


In this Section the experimental and predicted results that scattered over the previous 
subject areas are summarized and compared. These results are compiled in Table 9.1 and 
Fig. 9.1 sequential numbers in the first column of the table are intended to facilitate the 


assignment of the data in the tables to those in the text and in the figure. 


No. Property predicted experimental ratio 


(1) | -M (90°C) (Nm/m) [9.1] | 0.0047 0.00413 0.878 
[0.00363, 0.00463] | [0.772, 0.985] 


(2) | -M, (90°C) (Nm/m) [9.1] | 0.0040 1.032 
[0.908, 1.16] 
(3) | -M, (200°C) (Nm/m) [9.1] | 0.0145 0.0173 1.19 


[0.0128, 0.0217] | [0.882, 1.496] 


Table 9.1 Comparison of predicted properties by using swelling effect and measured data. 


Under (1) the bending moment for water-soaked silica disks at 90°C by Wiederhorn et al. 
[9.1] is given together with the prediction based on water content via IR-measurement 
(see Fig. 2.3). The experimental numbers in brackets represent +1 standard deviation 
(S.D.). 


Under (2) the same experimental data are compared with the predictions using water con- 


centrations at the surface as reported by Zouine et al. [9.2]. 
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9 Comparison between prediction and experiment 


Finally, (3) shows the results for the tests at 200°C (see Fig. 2.5). In Fig. 9.1 all these data 
are plotted as the ratio of measured results and predictions. It can be seen that these ratios 
do not significantly deviate from the value “1”, i.e. predictions and measurements are in 
rather good agreement. In Sections 7 and 8, the predictions of strength increase by 
swelling stresses [9.3] and residual deformation in 2-point bending tests were made based 


on the value Yoe™,,,= 14.4 cm?/mole as was concluded in [9.3]. 


The strength predictions are compiled in Fig. 9.2a by the solid curves (based on NRA- 
and IR-measurements) together with the data from [9.4]. In Table 9.2, the measured and 
predicted strength increases are compiled. Numbers in brackets represent +1 Standard 


Deviation. In addition, these results are introduced in Fig. 9.1. 
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Fig. 9.1 Comparison of predicted properties by using swelling effect and measured data. The serial 


numbers of the tests are given in parentheses, see Tables 9.1-9.3. 


For the second data set in Fig. 9.1, where strength increases are considered, the fact that 
about 10 % of the 15 test specimens did already fail under the 3 GPa pre-stressing has an 
influence on the strength of the surviving specimens as has been outlined in Section 7.1. 


The slightly reduced strength at 3 GPa is shown in Fig. 9.3 and introduced as test (7) in 
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9.1 Comparison of experimental results with predictions from swelling 


Fig. 9.1 and Table 9.2. In both cases, the predictions via the NRA-concentrations are 
introduced as red circles, bars, and numbers. However, if the predictions are based on the 
IR data, there are larger deviations (Table 9.2 and Fig. 9.2). The predictions based on the 
IR-data are given in blue colour. 

The residual curvatures after 2-point bending at 550°C [9.5] are shown in Fig. 9.2b. In 
Table 9.3, the measured and predicted residual curvatures are compiled. The numbers in 


brackets represent the 90% Confidence Intervals. 


Strength increase Ao, 


measured predicted measured/predicted 


No. | o (GPa) 
(4) 1 0.51 [0.297, 0.754] | 0.45 0.25 | 1.13 [0.66, 1.67] 2.04 [1.19, 3.0] 
2 1.21 [0.92, 1.49] 1.11 0.77 | 1.095 [0.83, 1.35] 1.57 [1.19, 1.94] 

(6) 3 2.205 [1.89, 2.54] | 1.89 1.48 | 1.17 [1.00, 1.34] 1.49 [1.28, 1.72] 
3 2.11 [1.80, 2.44] 1.89 1.48 | 1.12 [0.95, 1.29] 1.43 [1.21, 1.65] 


Table 9.2 Measured [9.4] and predicted strength increase Ao: [9.3], red: via NRA, blue: via IR. 
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Fig. 9.2 Predictions compared with measurements (red: So from NRA, blue: from IR-analysis), a) strength 


predictions, b) residual curvature for 550°C. 
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9 Comparison between prediction and experiment 


Residual curvature (1/R)/ Vt 
No. | o (MPa) measured predicted | measured/predicted 
(8) 309 0.056 [0.050, 0.063] 0.062 | 0.903 [0.806, 1.02] 
(9) 361 0.060 [0.055, 0.065] 0.074 0.81 [0.74, 0.88] 
(10) 493 0.101 [0.08197, 0.122] | 0.108 0.935 [0.76, 1.13] 


Table 9.3 Measured and predicted residual curvature (curve slopes: R in cm, fin h). 
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Fig. 9.3 Strength predictions compared with measurements (data point at 3 GPa corrected according to Fig. 
7.2b). 
Concluding remark 


In our opinion, the mechanical effect of the water/silica reaction is a volume increase, i.e. 
a swelling of the hydroxyl containing regions. This effect was experimentally proved via 
density measurements by Briickner [9.6, 9.7], Shackelford [9.8] and Shelby [9.9]. As 
illustrated in Fig. 9.1 there is good agreement between predictions on swelling and 
experimental results. This holds also for the strength and residual curvatures in Fig. 9.2. 
In this context, it has to be taken into account that the strength predictions only included 
the global swelling effect and ignored strength-increasing local swelling around the 


individual surface cracks. 


130 


9.2 Comparison with other explanatory approach 


Our description is based on equation (2.1.11) relating the hydroxyl concentration Sas the 
product of the silica-water-reaction (1.1.1) to a volume increase ¢,. The proportionality 


factor has to be understood as a material-specific property. 


Time dependent effects must of course occur during experiments. This is the case when 


the hydroxyl concentration depends on time: 


1. In soaking experiments there is first an increase of the molecular water at the 


surface (see Section 1) with the consequence QÀ, KÒ. 


2; For a fixed given volume element below the surface there is in addition a con- 


tinuously increase in water due to the diffusion process. 


3. Finally, the equilibrium is not spontaneously reached but may be approached with 


time delay due to finite reaction rates especially at lower temperatures. 


In all these cases the reason for mechanical responses is volume swelling and any time 


effect is caused by varying S-concentrations 
E(t) =Kx S(t) (9.1.1) 


For an application see Section 8.5, where recovery effects are discussed which have to be 


expected after bending tests when the temperature is not reduced. 


9.2 Comparison with other explanatory approach 


Some studies focus on stress relaxation as a possible explanation for the effects of the 
silica/water reaction. But effects occurring in the absence of pre-existing stresses as for 
instance the measured generation of bending moments in completely stress-free silica 


disks cannot be a consequence of any type of relaxation. 
Our conclusions from the point of view of the disk results [9.1] are: 


° Before the tests, the silica disks were etched in 5% volume fraction HF in order to 
remove 5-10um of glass, which might have contained a diffusion layer of water at 
the surface or might have residual stresses due to polishing. Consequently, stress 
relaxation could not be the reason for curvature since no stresses to relax were 


present. 
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9 Comparison between prediction and experiment 


° We have to conclude that the mechanical effect of the reaction (1.1.1) is an active 
process that generates strains and stresses depending on the mechanical boundary 
conditions. It is not a passive effect as is stress relaxation that needs pre-existing 


stresses. 


e On the other hand, Fig. 9.1 clearly shows that the predictions and measurements of 
curvature (or bending moments) are within 0.85-1.20 of the ideal agreement (dash- 
dotted line =1) for the case that the NRA results are used. We have to state that there 
is no free space for any further stress-generating effect. This confirms the assump- 
tion that the water species causing swelling strains is exclusively the hydroxyl S 


even at rather low temperatures. 


e Moreover that it can be concluded that no relaxation is present at 200°C. If 
relaxation would be present at 200°C, the generated stresses causing the observed 
curvature must have relaxed simultaneously to their creation and of course not any 


remaining curvature could result in contradiction to the observation. 


Conclusions with respect to time-dependent diffusion profiles and diffusivity, Sections 1 
and 3: 


e It could be shown that the change of the diffusion profile shapes with increasing 
time is a trivial consequence of the penetration limiting reaction parameter or mass 


transfer coefficient in combination with swelling stresses. 


e The same holds for the apparent time-dependence of the diffusivity. An inter- 


pretation in terms of relaxation effects is to our knowledge not available. 


It should be emphasized again that the effect of swelling is already proved experimentally 
by the density reduction due to hydroxyl generation [9.6-9.9]. Direct measurements of a 


relaxation are not known after knowledge of the authors up to now. 
Conclusion with respect to the strength and residual curvature results of silica fibers: 


° The explanation of the results in Fig. 9.2, given in [9.4,9.5], is based on a complete 
relaxation of pre-existing bending and tensile stresses in the thin surface layer con- 
taining water. After unloading, these layers should show an inverse stress. Unfortu- 
nately, the authors do not give a critical water concentration, at which the total stress 


relaxation should take place. 
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9.2 Comparison with other explanatory approach 


e Our competing explanation of the mechanical effect of reaction (1.1.1), suggested in 
[9.10,9.11,9.12], is local swelling of the glass volume due to hydroxyl generation 
as was experimentally proved in literature [9.6-9.9]. By mechanical boundary condi- 
tions, the swelling strains result in swelling stresses proportional to the amount of 
hydroxyl] concentration. This attempt is appropriate to predict quantitatively the me- 
chanical effect of swelling on strength and residual curvature [9.4, 9.5] and allows a 
qualitatively understanding of effects by simplified computations as were presented 


here (see Sections 6-8). 


e If relaxation exists, this effect is most likely to manifest itself in the tensile strength 
measurements. Table 9.2 makes this clear under the assumption that the IR data 
must be used for prediction. This could explain the remaining discrepancies between 


the prediction via the swelling effect and the measurement results. 


From the standpoint of our current investigations, we can interpret many of the experi- 
mentally available results in terms of volume swelling. Nevertheless, we consider relaxa- 
tion and also an interaction of swelling and relaxation to be quite conceivable. Swelling 
occurs under tension and pressure (to a lesser extent) and, in our opinion, a catalytic 
effect of water molecules could occur only in the tensile region but not under compressive 


stress. 


A possible reason for the existence of an apparent relaxation effect exclusively under 
tensile stress is addressed in Appendix A4. From modelling the breaking of SiO2 bonds 
by water molecules catalytically without involving a chemical reaction according to eq. 
(1.1.1) the possibility of a stress reduction appears to be simply possible. 

This was suggested very early on by the authors in a purely mechanical model analysis 
[A4.13, A4.14, A4.15]. 


In addition to the power law of subcritical crack growth, v oc A”, this model provided two 


main results: 
e The occurrence of a threshold value for permanent bond fracture by molecular water. 


e — On the other hand, it was implicitly found that only tensile stresses in the glass show 


this effect. 


e _ Furthermore, a threshold of the applied load should occur. 
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9 Comparison between prediction and experiment 


These results can be transferred to the strength behavior under the influence of molecular 
water by replacing the stress intensity factors at the crack with stresses in the glass. 
Especially a stress threshold value should be of importance for the application of this 


model on strength behavior of silica. 


The question that has to remain unanswered at the moment is whether the bond splitting 
will completely or partially decrease again after unloading. This question is irrelevant in 
the case of subcritical crack growth, since after each break of a bond the tip of the crack 
advances to the next bond and very large crack openings occur at their original location, 


which ruled out re-bonding. 
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10 Dissolution effect for liquid-water 
soaking 


10.1 Dissolution data from literature 


Although fused silica is a very corrosion resistant material, it shows at higher temperature 
dissolution effects in liquid water and, consequently, a surface removal. Diffusion layer 
thicknesses will be reduced compared to surfaces undergoing water vapour soaking. 


Consequently, swelling layers will be partially removed. 


An estimation of the dissolution effect needs the dissolution rate and the saturation con- 
centration. The thickness of the dissolved surface layer depends on the dissolution rate 
Dais and the concentration of dissolved SiO2 in the water, Csio2, described by the classical 


power law relation 


Dis = Das f _ Cso2 | (10.1.1) 


eq 


with the dissolution rate Daiso for Csio2=0 and the equilibrium concentration Cio2,.. = Ga. 


The exponent is often given as n>2. 


Dissolution rates Daiso for eq.(10.1.1), obtained from literature are plotted in Fig. 10.1a. 
The circles correspond to a soda-lime glass that should have the same dissolution behav- 
ior as silica [10.1]. The square shows the rate for silica at 150°C in a 0.0167 M CaCl: 
solution [10.2]. From these data we expect e.g. for 250°-260°C, where the specimens for 
the strength measurements in [10.3, 10.4] were prepared, a dissolution rate of Duis ~ 1 
nm/s (triangle in Fig. 10.1a) in buffered water environment. This high rate implies that 
long-time soaking in water must cause a strong surface removal. This would of course 


only be the case if always fresh water is available at the silica surface. 
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10 Dissolution effect for liquid-water soaking 


For 192h at 250°C it would result a surface removal of about 550 um. The pure water in 
our experiment and the results for the 0.0167 M CaCl: solution may differ somewhat, 


nevertheless the dissolution effect would be dramatic. 


10.2 Effective diffusion zone size 


In order to estimate the maximum effect of dissolution, any redeposition of dissolved 
material may be ignored in the following computation (i.e. a specimen in a large auto- 


clave is assumed). 
The swelling zone size b given as 


b=,/D,t (10.2.1) 


increases with respect to time tby a rate of 


db d [D D 
-L rg Ss tem 12v (10.2.2) 
dt d`” “Vt ?°b 


t 


If the dissolution rate in water of 250°C is abbreviated as Das=d W/dr, with half specimen 


thickness W, the total zone thickness results as 


A A een, (10.2.3) 


A final steady-state zone thickness is approached for dber/d FO instead of a continuously 
increasing thickness as suggested by eq.(10.2.1). This simply results in an equilibrium 


thickness reached after the time 


ee (10.2.4) 


Equation (10.2.3) reads in integrated form 


p= „Dt = Dist for t< teq (10.2.5) 
1 be, tt, 
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10.3 Saturation of the water environment 


For equilibrium: 


D, 
bg 45 (10.2.6) 


dis 


From Diis = 1 nm/s and Dy = 6.2x10°'® m?/s at 260°C it results a nearly negligible zone 


depth of only ber 0.16 um reached after a few minutes. 


10.3 Saturation of the water environment 


The rather dramatic rate effect is fortunately limited by the saturation uptake of SiO2 in 
water. The dissolution effect in a limited water environment is counteracted by the in- 
crease of the silica concentration in the liquid as results from eq.(10.1.1) and the increase 
of redeposition. The saturation concentration Cq increases with temperature, too. Figure 
10.1b shows an average curve that represents literature data for the equilibrium concen- 


tration Gq, compiled by Gunnarsson and Arnorsson [10.5]. 


In our autoclave the water volume was less the glass volume. The solubility at a soaking 
temperature of for instance 250°C results from Fig. 10.1b as 1250 ppm [10.5], i.e. Ca = 
0.00125. 


In the case of narrow packed specimens with a water volume smaller than the glass vol- 
ume as was the case for our soaking procedure in [10.3, 10.4], a simplified estimation of 
the dissolution can be based on the equilibrium concentration exclusively, resulting in the 


maximum possible effect. 


The maximum dissolved silica in the water containment of volume V, and density pw= 1 
g/cm? is mais = Vw Ce. The amount of silica removed in a surface layer of area Syiass, thick- 


ness dand density p =2.2 g/cm? is Mglass = Sglass XPgtass Xd. Since Mais=Mgtass, it follows 


daw Aæ c (10.3.1) 
S G 


glass P glass 


For our tests we used a very small autoclave of 30 cm? total volume. After placing a 
number of N specimens of each Vetass=0.3x0.4x4.5cm=0.54 cm? the autoclave was filled 


with water to maximum 2/3 of its height. This corresponds to the water volume 
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V,=20cm? — Nx Vetass. The total surface of the Nsilica specimen is Sgiass=2 N (0.3+0.4) 4.4 
cm? + 2Nx0.3x0.4 cm?=N6.4 cm’. 


By introducing the data estimated before into eq.(10.3.1) with the number of N=6-7 as 
was used for most of the soaking treatments, it results for the maximum possible layer 


removal: d=2.5um. 


The diffusion layer generated during 24h water soaking is with D=4.6x 10"! m?/s: b= 6.3 
um. Comparison of d and b makes clear that dissolution is not negligible especially for 


rather short soaking times. 
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Fig. 10.1 a) Rate for dissolution of silica in water (circles: data for soda-lime glass R6 with about equal 
dissolution behaviour as silica, reported by Perera and Doremus [10.1], square: result from Dove et 
al. [10.2], triangle: expected for 260°C), b) solubility of silica in water based on a literature review 


by Gunnarsson and Arnorsson [10.5]. 
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APPENDIX 


Al Crack problems 


Al.1 Single edge crack 


In the following considerations it is assumed for clearness that the equilibrium of reaction 
(1.1.1) is already reached. A single crack embedded in a swelling zone of height œ and 
constant volume swelling is shown in Fig. Al.1. The crack length, a, and the zone height 


are negligibly small compared to the specimen dimensions. 


Fig. Al.1 Single edge crack in a semi-infinite specimen embedded in a swelling zone. 


Zones of different heights were modelled in a finite element (FE) study. The case of a 
half-space was realized by a plate of width 50m and height 100. Solid continuum 
elements (8-node bi-quadratic) were chosen and the computations carried out with 
ABAQUS Version 6.8. The volume strain was replaced by the equivalent thermal 
problem by heating the swelling zone by A7 keeping 7=0 in the rest of the structure. 


This results in the strain 
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(a=thermal expansion coefficient). 


The shielding stress intensity factors from FE-computations normalized on the value 
K3(0), i.e. for the semi-infinite crack with w/aO [A1.1], are plotted in Fig. Al.2 as a 


function of the ratio of @/a. The FE-results given by the circles can be expressed by the 


relation 


with 


and different zone shape 
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Fig. Al.2 Shielding stress intensity factor as a function of the relative zone height w/a. 
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A1.2 Multiple edge cracks 


From the solution by eqs.(A1.1.2-A1.1.4), the shielding stress intensity factor for any 
swelling distribution can be determined using the Green’s function as outlined in [A1.2]. 


The shielding contribution of a zone increment of thickness dw’, (Fig. A1.3), is 


dK, = -y SE ep 1.18] 04237- 2a (A1.1.5) 
a a 


(l-v)vo! 


defining the Green’s function 
E o' o' 
h(o') = -y ———=1.18 exp] —1.18— || 0.4237- — (A1.1.6) 
OV yar a| al z 


This function allows to determine the shielding stress intensity factor for an arbitrarily 


given strain distribution &(®’) 


© 


Ku = fe) h(o')da' (A1.1.7) 


Fig. A1.3 Shape of a differential layer of thickness dœ’ in a zone with constant hydroxyl concentration. 


A1.2 Multiple edge cracks 


A1.2.1 Applied stress intensity factor 
Fig. Al.4a shows an array of identical edge cracks. The applied stress intensity factors 


are known from literature [A1.3]. We can represent these data by the simple relation 
[A1.4] 
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11/5 
K appi = VFR un 1.122433 ©) (A1.2.1) 
c 


with the constant solution for long cracks, a>c/3 


K c/2 (A1.2.2) 


appl = O appl 


independent of the crack depth a. This fact is very important. Equation (A1.2.1) is plotted 
in Fig. A1.4b. For half-penny shaped surface cracks, the solution deviates for a/c>0.5 less 
than 10% as can be seen from Isida and Tsuru [A1.5]. In the following computations we 


therefore consider the surface cracks as edge cracks. 


Whereas the single-crack solution increases continuously with increasing crack length, 


the stress intensity factor for the crack array remains constant for larger than a/c~ 1/3. 


See t Beek dedi s Cee ho nk 


© =const. 


Fig. Al.4 a) Array of edge cracks under remote load o, b) stress intensity factor solution. 
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A1.2 Multiple edge cracks 


A1.2.2 Shielding stress intensity factor 


Swelling behaviour of multiple cracks is the same as for a single crack. This is illustrated 
in Fig. Al.5a for a subcritically grown crack with heart-shaped zone front contour and in 
Fig. Al.5b for a crack array showing semi-circular zone fronts in the absence of external 
load. 


The crack array with the semi-circular zone front shape and the same swelling zone of 


height œ for each crack yields the shielding stress intensity factors of Fig. A1.6. 


EOP S A E E A E O cele S T 


7/ 


Fig. A1.5 a) Swelling zones of growing multiple cracks, b) modelling by a single crack introducing symmetry 
planes. 
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Fig. A1.6 Shielding stress intensity factors for multiple cracks (semi-circular zone ends). 


A1.3 Stress-enhanced swelling zones at short surface 
cracks 


In the case of a Jong crack of length a with a swelling zone height œ small compared to 
the crack length, w<<a, the heart-shaped zone at a crack tip yields a disappearing 
shielding stress intensity factor Ka=0 [Al.1]. The height of the contour for the hydro- 


static stress is in this case 


ey) (£) (A1.3.1) 
437 O, 


For the computation of the swelling zones ahead of very short surface cracks, œ~a, a 
Finite Element analysis was performed. Fig. A1.7 shows such a crack together with 
contours for constant hydrostatic stresses On, which under “weak swelling” conditions 


describe the shape of the swelling zones (dashed contours). This treatment is similar to 


150 


A1.3 Stress-enhanced swelling zones at short surface cracks 


the case of “weak phase transformation” assuming that the stress field caused by a phase 
transformation, e.g. in zirconia ceramics, does not affect the stress field triggering the 


transformation [A1.1, A1.6]. 


a 


Fig. Al.7 Contours of hydrostatic stress fields ahead of a short crack with depth a in a semi-infinite plate 
loaded by a tensile stress Oappı (solid curves). The dashed contours show the related long-crack 


solution with w/a50. 


The solid contours in Fig. A1.7 show the contours of hydrostatic stresses op for a tensile 


stress Capp. For a crack of depth a, the “applied stress intensity factor” results as 
K =1.1220,,,./74 (A1.3.2) 


In Fig. A1.8 the zone heights from FE are compared with the near-tip stresses as resulting 


from the long-crack solution 


cos(p/2), r= — (A1.3.3) 


K 
„2rr sin 7/3 


For w/a<'A, differences are negligible. 


o, =4(1+v) 
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Fig. Al.8 Hydrostatic stress at a crack under tension, comparison of short-crack results obtained via FE with 
the solution for a long crack according to McMeeking and Evans [A1.1]; contour for constant 


hydrostatic stress as insert. 


A1.4 Diffusion-controlled heart-shaped swelling zones 


A1.4.1 Zone shape 


The diffusion coefficient, D, which determines the rate of water diffusion into a solid 


depends on the stresses in the material and can be expressed by 
D=D, a re (A1.4.1) 


where ozis the hydrostatic stress, T the absolute temperature, R the universal gas constant, 
and AV, the activation volume for the diffusion of water in silica. Stress-enhanced 
diffusion occurs especially in the singular stress field of mechanically loaded cracks where 


o-o, 
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A1.4 Diffusion-controlled heart-shaped swelling zones 


If K denotes the stress intensity factor at the tip of a crack, then the singular near-tip 


hydrostatic stress field is given by the following equation: 


6, = idr vool?) (A1.4.2) 


„j2rr 


where rand @ are the the polar coordinates with origin at the crack tip. This results in 


heart-shaped diffusion and swelling zones around the crack tip, Fig. Al.9. 


On=const. 
=> gy=const. 


Fig. A1.9 Contour line for constant hydrostatic stress on. 


A1.4.2 Shielding stress intensity factor 


In order to compute the shielding stress intensity factor caused by an arbitrary radial 
distribution ¢,(z), the case of a constant swelling strain may be used, Fig. A1.10a. For this 
zone, McMeeking and Evans [A1.1] showed that 


K,, =0 (A1.4.3) 


Since this result does not depend on the zone height, œ, it results for variable swelling 


strains [A1.2] 


K, -feof Se) da'=0 (A1.4.4) 
0 E, do €=const 


i.e. eq.(A1.4.3) holds also in the case of any radial distribution of swelling, Fig. Al.10b. 


153 


Al Crack problems 


Fig. A1.10 a) Swelling zone with constant strain £v, b) zone increment with local strain &v(®). 


A1.5 Shielding stress intensity factor for cracks in a 
half-space 


A1.5.1 Edge crack 


Finite Element computations were carried out for a small edge crack in a rectangular 
bending bar, a/ W=0.005, under plane strain conditions, Fig. Al.11a. The plate height was 
chosen as H/ W=2 representing sufficiently the infinitely long bar. Swelling was assumed 
taking place from both specimen side-surfaces. Figure Al.11b represents the shielding 
stress intensity factor as a function of b/a. The FE results are introduced as the circles. 

In order to obtain an appropriate fitting relation the analytically known limit values have 
to be included. 

For b/a—0 it was derived in [A1.7]: 


E&E D E€, E 1.465b (A1.5.1) 
l-v 31-v) Ja 


K,, = -yY 


with y=0.255. The first term reflects the swelling zone along the crack surface. Its value 
was given by McMeeking and Evans [A1.1] and determined in [A1.7] also by FE-com- 


putations. 
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A1.5 Shielding stress intensity factor for cracks in a half-space 


of b) 


SswoVa | Eq.(A1.5.6) 


Fig. Al.11 a) Specimen dimensions, b) Finite Element results (circles) compared with the interpolation 


relation of Eq.(A1.5.6), dash-dotted curve: contribution of the side-surface swelling zone. 


Under the assumption of an erfc-shaped swelling strain distribution extending from the 


specimen surface, the swelling stresses in x-direction are 


T, =O, erta Š | (A1.5.2) 
; 2b 


with the compressive surface stress ©Osw,o <0. The contribution of this layer to the shielding 


stress intensity factor can be computed as 


K 


sh, surface 


= [o,n(x,a)a (A1.5.3) 


with the normal stress o, on the crack and the weight function A as given e.g. in [A1.4]. 
Figure Al.12 shows A for a crack in the half-space. For very thin swelling zones, b/a<<1, 
it holds 4=1.465/Va. 


The shielding stress intensity factor Ken surface for b<<a then reads 
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b/a>0 
r 


K sn surface = Os, 0 Va Ja 


which is about 12% larger than the solution for a step-shaped swelling distribution. 


For the limit b/a>>1 it holds 


Ka =o oF (74 (A1.5.5) 


with 7=1.1215. The solid curve in Fig. Al.11b represents an interpolation relation 
including the limit cases that fits best to the FE-results within 0<6/a<1: 


K,, =1.122,/za0,, toes + an (A1.5.6) 
É a a 


The dash-dotted curve in Fig. Al.11b represents the contribution caused by the side- 
surface swelling zone exclusively. For b/a>1 this contribution counts for more than 90% 


of the total shielding stress intensity factor. 


hva | 


mm 
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Fig. Al.12 Weight function A for an edge-cracked half-space. 
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A1.5 Shielding stress intensity factor for cracks in a half-space 


A1.5.2 Semi-circular crack 


The geometry and shielding stress intensity factors from FE-computations on a semi- 


circular crack in the semi-infinite body are shown in Fig. A1.13. 


The asymptotic limits are well known (upper limit for point B: 1.29, for point A: 1.17). 


The bold interpolation curves are given as 


Kan 2129. a0. le + 00642] (A1.5.7) 
a a 


Kaa IAT ty an 008 E +03 v] (A1.5.8) 
a a 


point (A) 


02 04 06 08 1 1.2 14 


b/a 


Fig. A1.13 Shielding stress intensity factors at the deepest point (A) and the surface points (B) of a semi- 
circular crack in a semi-infinite body as a function of the ratio of diffusion length bto the crack 


depth a. 
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A1.5.3 Semi-elliptical crack 


Two contributions make up the shielding stress intensity factor: one coming from the 
diffusion zone originating from the crack faces, Koni, the second originating from the 
external surfaces of the specimen, Ks. The total shielding stress intensity factor is the 


sum of both of these: 
Kona = Koni + Kon2 (A1.5.9) 


Crack-face zone: The diffusion zone developing from the crack faces is given by the 


following equation [A1.1], which holds both at point A and point B: 


rey ~ vb (A1.5.10) 
=p 


The shielding stress intensity factors for points (A) and (B) of the semi-elliptic crack, 
K2, are plotted in Fig. A1.14 versus a/c and b/a. It is obvious that the shielding stress 
intensity factor at the surface points (B) is roughly independent of a/c for about a/c>0.6. 


For small b/a, the stress intensity factor Ks,a2 is much less than Ays.. This solution 


1.13 000,° wnt (£) | 
c c (A1.5.11) 


a 1.65 
1+1 a“) 
C 


plotted in Fig. A1.14 can be written as 


Kap =Oovy7a Pas 


where 


s B, A1.5.12 
>» ( ) 


and 


a, ja B, (A1.5.13) 
F; -$1 - ; 
n=l a a a a 
1+[01186+ 01282 0.130%,+00288; )n 
c c c 


2b 
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Al.6 Deformed steady-state swelling zones 


0.42 
A, = 0.434 0.266 “| (A1.5.14) 
Cc 
0.42 
A, = 0.264 ©) (A1.5.15) 
c 
Bi=-0.37664, By=2.7314, By=—1.35521 
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Fig. A1.14 Influence of aspect ratio a/c and relative depth of the swelling zone b/a on the shielding stress 


intensity factors Ksn,a at the deepest and Ks, at the surface points of a semi-elliptical surface crack. 


A1.6 Deformed steady-state swelling zones 


The computations of the swelling zones are mostly carried out for diffusion in a motion- 
less coordinate system with a fixed crack tip, i.e. for an arrested crack. Whereas for a 
growing crack the transversal diffusion normal to the crack plane is hardly affected by the 
moving crack tip, the diffusion and the crack propagation compete in crack direction 
[A1.8, A1.9]. The consequence for strongly stress enhanced diffusion zones with cordial 
shape is a reduced ratio of the zone size in length direction m to the zone height œ (Fig. 
A1.15), which for an arrested crack is [A1.10] 
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r(0)_ 8 (A1.6.1) 


o(0) 27 


(n(0), @(0) stands for m(v=0) and &( v-0), respectively). For the stress intensity factor of 


the deformed zones one can use some single results from literature. 


arrested crack crack growing under 
steady state conditions 


Fig. Al.15 a) Swelling zone ahead the tip of an arrested crack, b) expected reduced zone length for a crack 


growing at a constant crack-growth rate v (schematic). 


8/\(27) 
| 
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Fig. Al.16 a) Swelling zones with different shape at the zone end, b) shielding stress intensity factors by 
eq.(A1.7.1) for the zones given in a). 
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A1.7 Edge-crack passing and leaving a swelling zone 


In [Al.11] it was found for crack rates v>10° m/s under the assumption that the 


diffusion normal to the crack-plane direction is not affected by the crack rate (v,=0): 


8 nM), 8 (A1.6.2) 


Tr 


Finally, Fig. A1.16 shows the coefficient wy for shielding stress intensity factors, 
eq.(A1.5.10). 


A1.7 Edge-crack passing and leaving a swelling zone 


During the soaking treatment under tension the crack-tip region is surrounded by a 
swelling zone. The shielding stress intensity factor Ky, caused by this zone was given in 
[A1.1] according to eq.(A1.5.10) 


K,, =-Vv =" Jo (A1.7.1) 


1 


with the coefficient y depending on the zone shape in the crack-tip region. 

When crack propagation has taken place during soaking, the swelling zone extends with 
crack length as illustrated in Fig. A1.17b. For the case of weak swelling as discussed by 
McMeeking and Evans [A1.1] the heart-shaped contour remains ahead of the tip and the 
coefficient is y=0.22. After high-temperature soaking, these swelling zones are main- 


tained at room temperature where the strength tests are performed. 


When the crack undergoes continuously increasing loading, the crack must extend and 
will escape of the initial swelling zone. Concequently, the coeffient y in eq.(A1.7.1) must 
become dependent on the crack extension Aa in the strength test. The red curve in Fig. 
Al.17c represents analytical results for a heart-shaped swelling zone we got by 
application of the McMeeking and Evans [A1.1] procedure. A maximum of w is reached 
at Aa/ =8/V27 =1.54 with wx0.3. In the case of a grown crack with extended swelling 


zone (black curve) it results w=0.33. 


At higher temperatures and after long time, also unloaded cracks may develop a swelling 
zone of thickness b by water diffusion into the crack faces. Figure Al.18a represents a 


crack with a circular zone front in the crack tip region. 
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If such a crack grows by an amount of Aa, it will escape from the initial swelling zone, as 
illustrated in Fig. Al.18a-Fig. Al.d for different zone ends. The shielding stress intensity 
factors for two limit cases of the zone shape near the crack tip are plotted in Fig. Al.18e. 
For the circular zone end (Fig. A1.18b), the shielding stress intensity factor is given by 
the squares [A1.12]. A straight zone end with /5=0 yields the dash-dotted curve. The 
two dependencies can be described by eq.(A1.7.1) where now y is a function of Aa/b. 


LLL 


0.15} 
0.1 


0.05 


Fig. A1.17 a) Swelling zone in the absence of crack propagation, b) grown crack for weak swelling, c) 
shielding stress intensity factor for crack extension through the initially swelling zone, red curve: 
crack passing the initial zone developing for an arrested crack, obtained by application of the 
procedure by McMeeking and Evans [A1.1]; black: grown crack after additional extension Aa 
obtained from FE-analysis, [A1.12]. 


The stress intensity factor solution for the case ™/b=¥2 (Fig. Al.18d) can be obtained by 
linear interpolation of the two limit cases resulting in the dashed curve of Fig. Al.18e. 


The coefficient y may be approximated by 


re 4 
PR 0.305 + 4Aa/b for Aa/bs (A1.7.2) 
0.34exp[-2,/Aa/b-4]  Aa/b>4 
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A1.7 Edge-crack passing and leaving a swelling zone 


It should be mentioned that the curves in Fig. Al.17c and Al.18e change abruptly caused 
by the assumption of a sharp layer boundary. For an arbitrarily varying strain distribution, 


the weight function or Green’s function method from [A1.2] can be used. 
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Fig. A1.18 a) Crack with a thick circular swelling layer of constant swelling strain, b) crack after additional 
extension Aa for a circular zone, c) for a rectangular zone, d) for a zone with m/5=1/2, e) variation 
of the shielding stress intensity factor for a crack growing in and outside the initial layer according 
to [A1.10] (dash-dotted curve) and [A1.12] (squares); dashed curve: interpolated from the 
continuous and dash-dotted curves for m/b=1/2. 


The stress intensity factor caused by a swelling layer of thickness db’ in the surface 


distance 203’ (Fig. A1.19a) undergoing the strain (b) results as 
d(K,,) = K,,(b'+db') — K,,(b') (A1.7.3) 


This gives rise for a Green’s- or “weight function” A defined as 
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en (A1.7.4) 


db' 


For any swelling distribution &(z)=e(b’), the related shielding stress intensity factor is 
then obtained by summing over all the zone increments dd’ multiplied with the local 


swelling strain 
Ky =+ Jeb ho)yab'. £ = 6(b'=0) (AL75) 
Eo 0 


In the absence of stress-enhanced diffusion, the swelling profiles are represented by the 


complementary error function. The volumetric strain at the location z=’ is in this case 


E = £, erfol wa (A1.7.6) 
D,t 


where & is the strain at the free surface, D the diffusion coefficient and ¢ the time after 


the first water contact. 
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Fig. A1.19 a) Green’s function as the shielding stress intensity factor of a thin swelling layer of thickness b’, b) 
stress intensity factors according to the dashed curve in (d) for an erfc-shaped strain distribution 


(insert); curves according to Fig. Al.18e. 
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A1.7 Edge-crack passing and leaving a swelling zone 


Equation (A1.7.5) has been evaluated for the case of m/b=”2. The resulting shielding 
stress intensity factors related to the erfc-shaped swelling profiles are shown in Fig. 


A1.19b. The numerical result for @/5=2 can be approximated by equation (A1.7.1) with 


0.305 (A1.7.7) 


2 
b 
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A2 Stresses at pores 


In order to allow application of continuum mechanics relations in terms of stresses and 
strains, the complicated SiO» micro structure is simplified as shown in Fig. A2.1. Figure 
A2.la shows a volume element of the silica network near the free surface containing SiO2 
rings of different size. The largest ones are indicated by the dashed circles. Figure A2.1b 


illustrates the simplification by spherical pores. 


Fig. A2.1 a) Silica network, b) large SiOz rings replaced by nano-pores. 


A2.1 Two-dimensional array of half-spherical surface 
pores 


Fig. A2.2a shows a symmetry element of a two-dimensional array of semi-spherical pores 


at a flat surface. The radius of the half-sphere is R and the distance to the next one 2H. 
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A2 Stresses at pores 


The applied stress is denoted as ©». The stresses oy along the equator between points (A) 
and (B) are represented in Fig. A2.2b for several values of H/R and Poisson’s ratio 
v=0.17. Fig. A2.2c shows the reciprocal stress at the surface and the deepest point as a 
function of R/H. 
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Fig. A2.2 a) Semi-spherical surface pore, b) stress distribution oy along the equator, c) reciprocal stresses, 


1/oy at the surface (B) and the deepest point (A), v=0.17. 
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A2.2 Broken bonds as tablet-shaped nano-pores 


The stresses oy at the surface points, (B), and the deepest point (A) can be approximated 
for R/H<0.8 by 


o| _ 1 (A2.1.1) 
2 4 
Fake) osi-o2s4 2) -025 £) 
H H 
o| _ 1 (A2.1.2) 
[07 


2; 
‘) 0.503-0.1 25 
H 


A2.2 Broken bonds as tablet-shaped nano-pores 


In the preceding considerations, it was implicitly assumed that no interaction of the bond- 
breaking events is present and consequently the hydroxyl concentration was considered 
to be constant over the cross section. For very high stresses, this cannot be guaranteed. 
On one hand side, a stress concentration occurs around a single broken bond since the 
surrounding bonds must transfer the load that was originally carried by this bond, Fig. 
A2.3. The increased load in the neighboured bonds causes an increased failure probability 
by stress-enhanced hydroxyl generation. The result would be an uneven distribution of 


the hydroxide content on the molecular level. 


A SiO, ring of the silica network may be described by a spherical (nano)-void in a 
continuous material (left part of Fig. A2.3). This allows a fracture mechanics treatment 
by an FE-study. According to the finite crack-tip curvature for cracks, the radius of such 


a void has to be expected in the order of p=0.5 nm [A2.1]. 


Under a load in zdirection, the red Si-O-bond may be cracked by reaction with mole- 
cular water (blue/red). Since the open bond cannot transfer load, the void now acts as a 
disk-like defect of length 2(p+Z) oriented perpendicular to the load direction. The stress 
increase Omax/O., at the notch root is shown in Fig. A2.4a as a function of the length L 


(circles). A rough representation of these data may be given by the simple relation 


Ti jags lisiri (A2.2.1) 
On P 
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Fig. A2.3 Left: Bond breaking under uni-axial tension by reaction with molecular water (attacked SiO-bond 


in red); right: open bond represented by a disk-like notch of diameter 2(p+Z). 
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Fig. A2.4 a) Stress concentration at the notch root for several values of notch length Z (circles: FE-results, 


curve: approximation by eq.(A2.2.1), b) displacements in the axis direction. 
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A2.2 Broken bonds as tablet-shaped nano-pores 


The maximum displacement in the centre is given in Fig. A2.4b. The straight-line depen- 


dency is represented by 


Z 


§ E =24+1275Ł (A2.2.2) 
p 


o 
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A3 Effect of surface cracks on the 
diffusion profile 


A3.1 Crack modelling 


In ground silica surfaces cracks are present even in the case of finishing by mechanical 
polishing or chemical etching. This has been outlined by Suratwala et al. [A3.1]. Such a 
surface state is schematically shown in Fig. A3.1a in side- and front view. The individual 
distributions of the depths a and the widths Z are approximated in [A3.1] by exponential 


distribution functions. 


For a simple explanation of the effect of such cracks on diffusion, let us replace the 
distributions of a and L by periodically arranged cracks, Fig. A3.1b, of average crack 


surface 
S =2(ax Ly (A3.1.1) 


in the surface element of length A and width B, where the factor 2 counts for two crack 


faces of each crack. The angle brackets (...) indicate averaging. 
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Fig. A3.1 a) Surface cracks due to surface treatment, front and side view, b) cracks of a) replaced by a 


periodical element with surface area AxB containing a single crack of mean crack surface S. 


A3.2 Diffusion profiles 


When the “unit cell” in Fig. A3.1b is exposed to a humid environment, water can diffuse 
into the whole surface given by the original surface area AxB and the crack surface S, 
Ax B+S, shown in Fig. A3.3a. The concentration of the water below the original surface, 


Q2), is given by the complementary error function 


C(z)=C, erid Z| (A3.2.1) 
2b 


with the concentration G at the surface. The depth b, the diffusion length, at which the 


concentration decreased to 50%, depends on diffusivity D and time /via 
b=./Dt (A3.2.2) 


The water profile normal to the crack faces is the same with yinstead of zin eq.(A3.2.1). 
Deviations occur in the crack-tip region. In a FE-study [A3.2], the diffusion problem for 


water-soaked specimens was solved. 


The water concentration contours are shown in Fig. A3.2a. It can be seen that the zone 
size ahead of the tip is reduced. In contrast to the profiles by eq.(A3.2.1), there is also a 


height reduction visible near the tip. 


174 


A3.2 Diffusion profiles 


Figure A3.2b gives the concentration profiles for the z-direction ahead of the crack tip 
and the profile by eq.(A3.2.1) for larger distances behind the tip region. For reasons of 
simplicity we approximate the water distribution in the diffusion zone by a layer of 
thickness b and a constant water concentration C as indicated by Fig. A3.2b. Near the 
crack tip, the missing water in y-direction is balanced by the water ahead the tip, result- 


ing for C/G=0.5 in the red step-shaped distribution in Fig. A3.2a. Figure A3.3a illus- 
trates the water distribution in the unit cell after a short time with b<< a. 
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Fig. A3.2 Water diffusion through crack surfaces of a semi-infinite crack, a) contour lines for constant water 


concentration, b) concentration profiles in length and height directions. 


175 


A3 Effect of surface cracks on the diffusion profile 


Fig. A3.3 Water diffusion into the crack surface containing a crack of depth a; a) water distribution for short 
times with diffusion zone b<<a, b) water diffusion zone according to eq.(A3.2.1) replaced by a 
zone of thickness band constant concentration showing the same amount of water, c) interacting 


diffusion profiles near the crack mouth (schematic), d) separation of overlapping diffusion profiles. 


The amount of water, penetrated into the glass is then approximated by a constant con- 
centration of C=C in a layer of thickness b and C=0 outside, Fig. A3.3b. For very short 
times, the profiles from the surface and from the crack faces interact only at the crack 
mouth, Fig. A3.3c. As long as this region is small, ie. for b<<a, it holds for the total 


amount of water in the unit cell 
m. = C,b[AB + 2La] (A3.2.3) 


In terms of an apparent layer thickness bap for very short time £>0, eq.(A3.2.3) reads 


<S> 
bap, = Disa dot = if ae AB (A3.2.4) 


For not too large layer thicknesses, let say b< 1/3 Min(a, B/2), the moderate diffusion 
interactions may be separated in diffusion layers from the specimen surface over the 
whole width B and from the crack surface over the reduced distance of a-b as illustrated 


in Fig. A3.3d. In this case we obtain 
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A3.3 Apparent diffusivity 


2L <S> 
Dis =o(1+ 26 (a-)]-of14 AB a-b/a)| (A3.2.5) 


This dependency is shown in Fig. A3.4. 
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Fig. A3.4 Apparent diffusion layer thickness bpp as a function of time. 


A3.3 Apparent diffusivity 


Equations (A3.2.4) and (A3.2.5) allow the determination of an apparent diffusivity Dapp, 


defined via 


= /D_t (A3.3.1) 


D, sofish) (A3.3.2) 
pp AB 
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Two limit cases can be identified. For very short times, it results from eq.(A3.3.2) 


2 
<S> 
Do > Dizisi = oft a AB ) (A3.3.3) 


For long times fand thick zones b, the limit diffusivity is trivially 


D =D =D (A3.3.4) 


app,© ~ ~~ app,t—>% 


Equation (A3.3.2) makes clear that the apparent diffusivity resulting from measurements 


depends on time as is plotted in the representation of Fig. A3.5. 
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Fig. A3.5 Apparent diffusivity Dap as a function of time for several ratios of crack surface <S> to specimen 


surface AxB, arrow indicates the asymptotically reached diffusivity Dapp >o. 


A3.4 Interpolation of limit case solutions 


Having in mind that the crack sizes and distances are continuously distributed [A3.1], the 
results from Section A3.2 may be “smoothed” and expressed by interpolation relations. 


Usual descriptions by exponential functions result in 
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<S> 
Dapp sois = exp) (A3.4.1) 
my b= fi + = (1- anha) ) (A3.4.2) 


Eq.(A3.4.2) 
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a) Interpolation of limit solutions according to eqs.(A3.4.1) and (A3.4.2), b) interpolations for 


different crack parameters, c) related apparent diffusivities. 
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Figure A3.6a shows these interpolations for the case of <S>/(AB)=2. The straight lines 
are again the limit solutions for £>0 and +o. The dashed curve represents eq.(A3.4.1) 
and the solid curve eq.(A3.4.2). Use of (A3.4.2) may be recommended. Figure A3.6b 
shows the interpolations for different parameters <S>/(AB). Finally, Fig. A3.6c gives the 
apparent diffusivity from eq.(A3.4.2) as 


D,,, = of! i —- (1-tanh(JDt ia] (A3.4.3) 


A3.5 Comparison with experimental results 


Results of the diffusivity as a function of time are available in the literature. Profiles of 
time-dependent residual water content by stepwise surface removal were observed by 
Davis and Tomozawa [A3.3] under water vapour pressure of 355mm Hg at 350°C. The 
results given in Fig. A3.7a, showed reduced normalized diffusion depths with increasing 
soaking time. This behaviour is plotted again in Fig. A3.7b in terms of an effective diffu- 


sivity Des (red circles). 
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Fig. A3.7 a) Profiles of residual water content during stepwise surface removal by Davis and Tomozawa 
[A3.3], b) effective diffusivities by Davis and Tomozawa [A3.3] (red circles) and by Oehler and 


Tomozawa [A3.4] at 250°C and 39 bar water vapour pressure (blue circles). 


Measurements of effective diffusivities by Oehler and Tomozawa [A3.4] at 250°C and 
39 bar water vapour pressure are introduced by the blue circles. In both cases a clear 
decrease with time is evident. This effect was discussed in [A3.5] in terms of swelling 
stresses. In case of the results by Davis and Tomozawa [A3.3] it has to be considered that 
under the rather low vapour pressure swelling stresses are not very high. As Fig. A3.6c 
shows, the curves in Fig. A3.7a may also be interpreted as the consequence of surface 


cracks. That will partially hold for the results in [A3.4], too. 
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A4 Stress reduction and subcritical crack 
growth 


A4.1 Analogy to a mechanical model for subcritical 


crack growth 


The possibility of the breaking of SiO bonds by water molecules without a chemical 
reaction according to eq.(1.1.1) appears to be quite trivial. The H>O molecule is retained 
and only reacts analytically. This was suggested very early on by the authors in a purely 
mechanical model analysis [A4.1, A4.2, A4.3]. A single bond in front of the crack tip was 
considered, which had to be broken. The effect of the neighbouring bonds was described 
as a continuum. The consideration in [A4.1-A4.3] made it possible to explain the 
subcritical crack growth including the threshold, which, however, also arises freely for 
the separation process via the chemical reaction [A4.4, A4.5]. The main statements of 


this analysis will be briefly addressed here. 


Figure A4.la illustrates the isolated bond which is prestressed by the springs (1), 
reflecting the linear-elastic continuum, Fig. A4.1b. The spring affected by a water 
molecule is shown in Fig. A4.1c. Figure A4.1d represents two load-displacement curves 
for bond interactions without (dashed curve) and with water (continuous curve). They are 
replaced by a non-linear spring (2), i.e. by a non-linear relation between the force Fand 
the bond displacement 5. For the following considerations a Morse relation is used as is 


common for intermolecular interactions. If F% denotes the maximum force, then 
F = 4F,[exp(-56) — exp(—2s6)], (A4.1.1) 


where sis a measure of the range of bond forces. The load-displacement characteristic of 
a single bond under a thermal activation event is shown by the 778 diagram in Fig. A4.1d 
and in Fig. A4.2 for different external loading (increasing load from a) to d)). In the case 


of a thermal event splitting the bond, the external load becomes reduced with increasing 
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distance of the two bond partners, because the situation shown in Fig. A4.1 is displace- 


ment controlled. 


Figure A4.2d describes a critical load situation. The intersections 6“ and 6 are identical 
and the straight line is the tangent to the bond force curve. This situation corresponds to 
catastrophic failure of the bond whereas the tangent in Fig. A4.2b represents a threshold 
in bond breaking. In Fig. A4.2e the characteristic displacements are related to the 
macroscopic values of the threshold and the critical stress intensity factor values Kj, and 


Kic, and also for the stresses on and o., respectively. This suggests 


6 Ky, a Fn Fn (A4.1.2) 
ô, K,, F; ©. 


Fig. A4.1 a) Model of the loading situation of a stretched bond in the glass network (schematically), b) silica 
network replaced by a continuum, linear springs (1): elastic surrounding, non-linear spring (2): 
interaction forces between two bond partners, c) elastic response of silica network replaced by 
linear springs, reacting bond modelled by a nonlinear spring, d) Morse-type force-displacement 


curves. 
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A4.1 Analogy to a mechanical model for subcritical crack growth 


The critical stress in eq.(A4.1.2) will be equal to 0.=00 in the absence of water where oo 
is the ideal or theoretical strength. For ©. see for instance [A4.6, A4.7] and for oo [A4.8, 
A4.9]. For a small external load (Fig. A4.2a) the bond is displaced by the amount 5°”, 
which represents the state of equilibrium. Even when thermal activation events with high 
energies occur, no additional state of equilibrium will be possible, and the bond will be 
re-established without a new thermal event. In this case the bond cannot be broken per- 


manently. 


sag 


ah ös Ô 


Fig. A4.2 Load-displacement behaviour of an individual bond: a)-d) increasing load, e) definition of critical 


and threshold forces, stress intensity factors and stresses. 


At an increased external load (Fig. A4.2b) there exists a second displacement where the 
straight line and the Morse curve have a common tangent. At a higher load (Fig. A4.2c) 


we find three intersections of the two curves. If in this case a thermal activation event 
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with energy of at least AUoccurs, the bond will find a new position of equilibrium at 5°, 
and the crack can extend by a lattice distance. But this new position with its negligible 
bond interactions is not necessarily a permanent one. If further thermal events with 
energies of at least AW appear at the considered location, which lead to a decrease in the 


distance of the bond partners, the bond can be re-established. 


A4.2 Bond-fracture frequency 


The probability of occurrence of thermal activation events with energies high enough for 


bond breaking is described by the Boltzmann distribution, namely 


= exp(-AU /KT)» (A4.2.1) 


2|2I 


In this expression N / N, (No=number of bonds) is the fraction of events with energies 
greater than the energy barrier AU and & is Boltzmann's constant. For describing the 
frequency of bond-fracture events knowledge of the potential barriers AU and A Wis ne- 
cessary. Numerical computations are reported in [A4.2]. 


During a time span of Arthe number of cracked bonds AN is in a statistical sense 
AN / N, = exp(-AU /KT)[1-exp(-AW /kT)] (A4.2.2) 


From (A4.2.1) it can be seen that as the stress on a crack decreases, the expression in 
brackets [ ] becomes zero, and the share of open bonds vanishes. The related stress value 


is the threshold stress ou. 


A decrease in the Young’s modulus is to be expected from the fact that broken bonds 
cannot transfer any or at least greatly reduced loads. The consequence is a decrease in the 


tensile stresses in the glass. 
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FINAL REMARKS AND OUTLOOK 


Extremely high stresses at crack tips and the water entrance in the case of humid envi- 
ronments cause high hydroxyl concentrations (Sections 5, 6) accompanied by large vol- 
ume swelling and strong damage. In this Section, a first insight into the problematic na- 
ture may be tried for cracks and slender notches applying purely elastic material behav- 


iour. In actual studies, the authors investigate the theme by: 


a) Use of a micro-structurally motivated approach, where the crack tip regions are con- 
sidered as slender notches with root radii in the order of the average radius of the SiO2 


rings, 


b) Alternatively by application of well-established models borrowed from Elastic-Plastic 
Fracture Mechanics (EPFM). 


c) In our considerations on fiber strength, we made our predictions for the high- 
temperature range of >450°C. In this region, the equilibrium constant is very large, so 
that molecular water could hardly contribute to the total amount of water. For predic- 
tions in the low-temperature region, we had assumed so far that the molecular water in 
the glass is located in molecular size holes in the silica glass network where the water 
does not contribute to the volume of the glass. Consequently, molecular water should not 
contribute stress effects, with the consequence of &c=0. To our knowledge, no trustwor- 
thy measurements of volume changes by molecular water &c#0 are actually available. 
Therefore, we assumed disappearing volume changes by molecular water. In terms of 
molar units for Sand C, this relation would read in the simplest case of a mixing rule in 
the form of eq.(2.1.5b). 


In a recent paper by Garofalini et al. [1] some first results on this problem were obtained 
by MD-computations. From this paper the volume changes for a water concentration of 
Cv=5 mol-% were taken and re-plotted in Fig. la. On the abscissa the hydroxyl concen- 
trations are normalized on their maximum value of S"®=2C, and the ordinates on the 


value & obtained for the volume strain at S=2 Cy. An approximate description reads 
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ee 0) 
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For S/(2 Cx») 2 0.25, the data may be represented by a straight line as introduced as the 
dash-dotted line. The related description is given by the first term in eq.(1), i.e. by 


oe, Soc) (2) 
2C 


Finally, Fig. 1b represents the same data by using the ratio S/Cas the abscissa. The inter- 
section of the curve with the abscissa is found at S/C=0.22. According to eq.(1.1.3) this 
corresponds to an estimated temperature of 0=-14°C. The straight line approximation for 
S/(2Cy)20.25 indicates that the temperature should be 0>45°C for the simplified descrip- 
tion. Summarising these results via molecular-dynamics calculations we can conclude, 
that an effect of molecular water is only of interest for temperatures near room tempera- 
ture and can be ignored in all cases considered in this booklet. On the other hand, it 
should be noted that MD calculations can never replace missing experimental results. 


However, they allow at least qualitative statements about trends. 
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Fig. 1 a) Relative volume strain vs. relative hydroxyl concentration (curve: eq.(1)), dash-dotted line 


according to eq.(2), b) relative volume strain vs. equilibrium constant S/C. 
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d) In some experiments, the effects of swelling on samples soaked in liquid water were 
significantly stronger than on samples in water vapour at the same elevated temperature 
and saturation pressure. In Fig. 2a, the predictions on the basis of hydroxyl concentra- 
tions in liquid water, derived from measurements of total water concentration by Zouine 
et al. [2], are shown as the blue circles. The red squares indicate the results from the disk 
tests in Section 2. From measurements of surface swelling stresses from crack- 
terminating angles in DCDC-specimens [3, 4], the blue upright triangles were obtained. 
Inert strength measurements on water- and vapour-soaked bending bars, [5], yielded the 
black downward facing triangles. In all of the results, the values under vapour pressure 
are below those in water. It is somewhat astonishing, that at 90°C the two results are 
identical. This is the reason why the data points were slightly shifted from one another 


with respect to the temperature axis. 
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Fig. 2 a) Predicted (blue circles) and indirectly measured swelling stress results from disk curvature (red 


squares), terminating angles (blue triangles), and inert strengths (black downward triangles). The 
squares for 90°C are slightly displaced, b) Hydroxyl water concentration ŞS normalized on satura- 


tion value ‚So vs. time. 
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Figure 2b represents again results from Section 1. The hydroxyl] concentrations for meas- 
urements in water vapour are plotted in normalized form with the ordinate normalized 


on the individual saturation value So obtained by extrapolation to >”. 


As a clear trend it can be seen that the normalized surface concentrations S/So for the 
same time decrease with increasing temperature. However, it should be noted that in 
addition to the temperatures, different pressures and glasses with possibly different prop- 
erties were used. Whereas the results by Helmich and Rauch [6] (triangles) show satura- 
tion after about 50-100 h, saturation for the results by Oehler and Tomozawa [7] at 
250°C/39.5 bar and by Davis and Tomozawa [8] at 350°C/355 Torr is significantly not 
reached below 1000 h vapour soaking. The reason for these differences is assumed to be 
affected by a variation of mass-transfer conditions at the surface (see Sections 1.2 and 


1.3). This topic has to be deepened in future work. 
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SCHRIFTENREIHE DES INSTITUTS FUR ANGEWANDTE MATERIALIEN 


Water diffusing into silica surfaces gives rise for several effects on diffusion behaviour 
and mechanical properties. Water added to silica glass increases its specific volume 
so that the silica expands near the surface. Mechanical boundary conditions give rise 
for compressive “swelling stresses”. This fact provides a tool for the interpretation 
of many experimental observations from literature. Some of the most interesting 
observations discussed in this booklet are: 


- Density reduction with increasing water concentration. 
- Time-dependent water concentration at silica surfaces. 
- Decrease of measured water diffusivity with time. 

- Strength and toughness increase due to swelling. 

- Residual deformation after bending. 
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